Changeset 2413 for LMDZ5/trunk/libf/phylmd/radlwsw_m.F90

Timestamp:

Dec 18, 2015, 8:27:39 PM (8 years ago)

Author:

oboucher

Message:

limiting iflag_rrtm=0 to NSW=2 case (NSW=4 and 6 were possible with slightly different averaging procedures)
limiting iflag_rrtm choices to 0 and 1
limiting iflag_albedo choices to 0 and 1
cleaning up surface ocean albedo parametrisations
fmagic and pmagic corrections are now done in surf_ocean_albedo.F90
fmagic and pmagic corrections also possible for iflag_rrtm=1 now

File:

• LMDZ5/trunk/libf/phylmd/radlwsw_m.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/radlwsw_m.F90

@@ -426 +426 @@
 
 !albedo SB >>>
-!      PALBD(i,1) = alb1(iof+i)
-!      PALBD(i,2) = alb2(iof+i)
-!         PALBD_NEW(i,1) = alb1(iof+i)   !!!!! A REVOIR (MPL) PALBD_NEW en
-!         fonction bdes SW
-!         do kk=2,NSW
-!           PALBD_NEW(i,kk) = alb2(iof+i)
-!         enddo
-!      PALBP(i,1) = alb1(iof+i)
-!      PALBP(i,2) = alb2(iof+i)
-!
-!         PALBP_NEW(i,1) = alb1(iof+i)     !!!!! A REVOIR (MPL) PALBP_NEW en
-!         fonction bdes SW
-!         do kk=2,NSW
-!           PALBP_NEW(i,kk) = alb2(iof+i)
-!         enddo
-
-      if(iflag_rrtm==0)then
-        select case(nsw)
-        case(2)
-          PALBD(i,1)=alb_dif(iof+i,1)
-          PALBD(i,2)=alb_dif(iof+i,2)
-          PALBP(i,1)=alb_dir(iof+i,1)
-          PALBP(i,2)=alb_dir(iof+i,2)
-        case(4)
-          PALBD(i,1)=alb_dif(iof+i,1)
-          PALBD(i,2)=(alb_dif(iof+i,2)*SFRWL(2)+alb_dif(iof+i,3)*SFRWL(3) &
-                 +alb_dif(iof+i,4)*SFRWL(4))/(SFRWL(2)+SFRWL(3)+SFRWL(4))
-          PALBP(i,1)=alb_dir(iof+i,1)
-          PALBP(i,2)=(alb_dir(iof+i,2)*SFRWL(2)+alb_dir(iof+i,3)*SFRWL(3) &
-                 +alb_dir(iof+i,4)*SFRWL(4))/(SFRWL(2)+SFRWL(3)+SFRWL(4))
-        case(6)
-          PALBD(i,1)=(alb_dif(iof+i,1)*SFRWL(1)+alb_dif(iof+i,2)*SFRWL(2) &
-                 +alb_dif(iof+i,3)*SFRWL(3))/(SFRWL(1)+SFRWL(2)+SFRWL(3))
-          PALBD(i,2)=(alb_dif(iof+i,4)*SFRWL(4)+alb_dif(iof+i,5)*SFRWL(5) &
-                 +alb_dif(iof+i,6)*SFRWL(6))/(SFRWL(4)+SFRWL(5)+SFRWL(6))
-          PALBP(i,1)=(alb_dir(iof+i,1)*SFRWL(1)+alb_dir(iof+i,2)*SFRWL(2)  &
-                 +alb_dir(iof+i,3)*SFRWL(3))/(SFRWL(1)+SFRWL(2)+SFRWL(3))
-          PALBP(i,2)=(alb_dir(iof+i,4)*SFRWL(4)+alb_dir(iof+i,5)*SFRWL(5)  &
-                 +alb_dir(iof+i,6)*SFRWL(6))/(SFRWL(4)+SFRWL(5)+SFRWL(6))
-        end select
-      elseif(iflag_rrtm==1)then
+!
+      IF (iflag_rrtm==0) THEN
+!
+        PALBD(i,1)=alb_dif(iof+i,1)
+        PALBD(i,2)=alb_dif(iof+i,2)
+        PALBP(i,1)=alb_dir(iof+i,1)
+        PALBP(i,2)=alb_dir(iof+i,2)
+!
+      ELSEIF (iflag_rrtm==1) THEn
+!
         DO kk=1,NSW
-         PALBD_NEW(i,kk)=alb_dif(iof+i,kk)
-         PALBP_NEW(i,kk)=alb_dir(iof+i,kk)
+          PALBD_NEW(i,kk)=alb_dif(iof+i,kk)
+          PALBP_NEW(i,kk)=alb_dir(iof+i,kk)
         ENDDO
-      endif
+!
+      ENDIF
 !albedo SB <<<
-
-
 
 
@@ -666 +634 @@
        ENDIF
 
-              
-          DO i=1,kdlon
-          DO k=1,kflev+1
-         ZSWFT0_i(1:klon,k) = ZFSDN0(1:klon,k)-ZFSUP0(1:klon,k)
-         ZLWFT0_i(1:klon,k)=-ZFLDN0(1:klon,k)-ZFLUP0(1:klon,k)
+       ZSWFT0_i(:,:) = ZFSDN0(:,:)-ZFSUP0(:,:)
+       ZLWFT0_i(:,:) =-ZFLDN0(:,:)-ZFLUP0(:,:)
+
+       DO i=1,kdlon
+       DO k=1,kflev+1
 !        print *,'iof i k klon klev=',iof,i,k,klon,klev
          lwdn0 ( iof+i,k)   = ZFLDN0 ( i,k)
@@ -680 +648 @@
          swup0 ( iof+i,k)   = ZFSUP0 ( i,k)
          swup  ( iof+i,k)   = ZFSUP  ( i,k)
-          ENDDO  
-          ENDDO  
+       ENDDO  
+       ENDDO  
 !          print*,'SW_AR4 ZFSDN0 1 , klev:',ZFSDN0(1:klon,1),ZFSDN0(1:klon,klev)
 !          print*,'SW_AR4 swdn0  1 , klev:',swdn0(1:klon,1),swdn0(1:klon,klev)