Changeset 2360 for LMDZ5/trunk/libf/phylmd/physiq.F90

Timestamp:

Sep 5, 2015, 12:06:39 PM (9 years ago)

Author:

oboucher

Message:

Small structural changes to rmu0 calculations
Now the ocean albedo sees the high-resolution rmu0
if iflag_cycle_diurne=2

File:

• LMDZ5/trunk/libf/phylmd/physiq.F90 (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -1820 +1820 @@
   IF (abs(solarlong0-1000.)<1.e-4) then
      call zenang_an(iflag_cycle_diurne.GE.1,jH_cur,rlat,rlon,rmu0,fract)
+     JrNt = 1.0
   ELSE
   ! recode par Olivier Boucher en sept 2015
@@ -1826 +1827 @@
      !  Sans cycle diurne
         CALL angle(zlongi, rlat, fract, rmu0)
+        swradcorr = 1.0
+        JrNt = 1.0
+        zrmu0 = rmu0
      CASE(1)  
      !  Avec cycle diurne sans application des poids
@@ -1832 +1836 @@
         zdtime=dtime*REAL(radpas) ! pas de temps du rayonnement (s)
         CALL zenang(zlongi,jH_cur,0.0,zdtime,rlat,rlon,rmu0,fract)
+        zrmu0 = rmu0
+        swradcorr = 1.0
+     ! Calcul du flag jour-nuit
+        JrNt = 0.0
+        WHERE (fract.GT.0.0) JrNt = 1.0
      CASE(2)  
      !  Avec cycle diurne sans application des poids
@@ -1842 +1851 @@
         zdtime2=dtime*REAL(radpas-MOD(itaprad,4)-1) 
         CALL zenang(zlongi,jH_cur,zdtime1,zdtime2,rlat,rlon,rmu0,fract)
+     !
+     ! Calcul des poids
+     !
+        zdtime1=-dtime !--on corrige le rayonnement pour representer le
+        zdtime2=0.0    !--pas de temps de la physique qui se termine
+        CALL zenang(zlongi,jH_cur,zdtime1,zdtime2,rlat,rlon,zrmu0,zfract)
+        swradcorr = 0.0
+        WHERE (rmu0.GE.1.e-10 .OR. fract.GE.1.e-10) swradcorr=zfract/fract*zrmu0/rmu0
+     ! Calcul du flag jour-nuit
+        JrNt = 0.0
+        WHERE (zfract.GT.0.0) JrNt = 1.0 
      END SELECT 
   ENDIF
@@ -1916 +1936 @@
           dtime,     date0,     itap,    days_elapsed+1, &
           debut,     lafin, &
-          rlon,      rlat,      rugoro,  rmu0,      &
+          rlon,      rlat,      rugoro,  zrmu0,      &
           zsig,      sollwdown, pphi,    cldt,      &
           rain_fall, snow_fall, solsw,   sollw,     &
@@ -3719 +3739 @@
 
   !
-  ! Calculer les poids a appliquer sur le SW 
-  ! sortie JrNt = jour-nuit 
-  ! recode par Olivier Boucher en sept 2015
-  ! 
-
-  SELECT CASE (iflag_cycle_diurne)
-  ! 
-  CASE(0)
-  !  Sans cycle diurne
-     swradcorr=1.0
-     JrNt = 1.0
-  CASE(1)
-  !  Avec cycle diurne sans les poids
-     swradcorr=1.0
-     JrNt=0.0
-     WHERE (fract.GT.0.0) JrNt=1.0
-  CASE(2)
-  !  Avec cycle diurne et les poids
-     zdtime1=-dtime !--on corrige le rayonnement pour representer le
-     zdtime2=0.0    !--pas de temps de la physique qui se termine
-     CALL zenang(zlongi,jH_cur,zdtime1,zdtime2,rlat,rlon,zrmu0,zfract)
-     swradcorr=0.0
-     WHERE (rmu0.GE.1.e-10 .OR. fract.GE.1.e-10) swradcorr=zfract/fract*zrmu0/rmu0
-     JrNt=0.0
-     WHERE (zfract.GT.0.0) JrNt=1.0 
-  END SELECT
-
-  !
   ! Corriger les flux SW pour le cycle diurne ameliore
   ! recode par Olivier Boucher en sept 2015