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

Timestamp:

Sep 21, 2015, 10:41:04 PM (9 years ago)

Author:

oboucher

Message:

Second batch of changes for diagnosing SW radiative
diagnostics at the physics timestep resolution
radsol is now computed in physiq.F90

File:

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

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -590 +590 @@
   !
   REAL dist, rmu0(klon), fract(klon)
-  REAL zrmu0(klon), zfract(klon), swradcorr(klon)
+  REAL zrmu0(klon), zfract(klon)
   REAL zdtime, zdtime1, zdtime2, zlongi
   !
@@ -1846 +1846 @@
      !  Comme cette routine est appele a tous les pas de temps de la physique
      !  meme si le rayonnement n'est pas appele je remonte en arriere les
-     !  radpas-1 pas de temps suivant. Petite ruse pour prendre en compte le
-     !  premier pas de temps la physique ou itaprad=0 
+     !  radpas-1 pas de temps suivant. Petite ruse avec MOD pour prendre en
+     !  compte le premier pas de temps de la physique pendant lequel itaprad=0 
         zdtime1=dtime*REAL(-MOD(itaprad,4)-1)      
         zdtime2=dtime*REAL(radpas-MOD(itaprad,4)-1) 
@@ -3586 +3586 @@
              topsw0,toplw0,solsw0,sollw0, &
              lwdn0, lwdn, lwup0, lwup,  &
-             zswdn0, zswdn, zswup0, zswup, &
+             swdn0, swdn, swup0, swup, &
              ok_ade, ok_aie, &
              tau_aero, piz_aero, cg_aero, &
@@ -3627 +3627 @@
              zqsat, flwc, fiwc, &
              ref_liq, ref_ice, ref_liq_pi, ref_ice_pi, &
-             heat,heat0,cool,cool0,radsol,albpla, &
+             heat,heat0,cool,cool0,albpla, &
              topsw,toplw,solsw,sollw, &
              sollwdown, &
              topsw0,toplw0,solsw0,sollw0, &
              lwdn0, lwdn, lwup0, lwup,  &
-             zswdn0, zswdn, zswup0, zswup, &
+             swdn0, swdn, swup0, swup, &
              topswad_aero, solswad_aero, &
              topswai_aero, solswai_aero, &
@@ -3683 +3683 @@
                    zqsat, flwc, fiwc, &
                    ref_liq, ref_ice, ref_liq_pi, ref_ice_pi, &
-                   heatp,heat0p,coolp,cool0p,radsolp,albplap, &
+                   heatp,heat0p,coolp,cool0p,albplap, &
                    topswp,toplwp,solswp,sollwp, &
                    sollwdownp, &
                    topsw0p,toplw0p,solsw0p,sollw0p, &
                    lwdn0p, lwdnp, lwup0p, lwupp,  &
-                   zswdn0p, zswdnp, zswup0p, zswupp, &
+                   swdn0p, swdnp, swup0p, swupp, &
                    topswad_aerop, solswad_aerop, &
                    topswai_aerop, solswai_aerop, &
@@ -3730 +3730 @@
      swup=0.    ! MPL 27102011 pour les fichiers AMMA_profiles et AMMA_scalars
      swup0=0.
-     zswdn=0.
-     zswdn0=0.
      lwup=0.
      lwup0=0.
@@ -3739 +3737 @@
 
   !
-  ! Corriger les flux SW pour le cycle diurne ameliore
+  ! Calculer radsol a l'exterieur de radlwsw
+  ! pour prendre en compte le cycle diurne
   ! recode par Olivier Boucher en sept 2015
   !
-
-  DO k=1, klev+1
-    swdn0(:,k)=swradcorr(:)*zswdn0(:,k)
-    swdn(:,k) =swradcorr(:)*zswdn(:,k)
-    swup0(:,k)=swradcorr(:)*zswup0(:,k)
-    swup(:,k) =swradcorr(:)*zswup(:,k)
-  ENDDO
+  radsol=solsw*swradcorr+sollw
   if (ok_4xCO2atm) then
-  DO k=1, klev+1
-    swdn0p(:,k)=swradcorr(:)*zswdn0p(:,k)
-    swdnp(:,k) =swradcorr(:)*zswdnp(:,k)
-    swup0p(:,k)=swradcorr(:)*zswup0p(:,k)
-    swupp(:,k) =swradcorr(:)*zswupp(:,k)
-  ENDDO
+    radsolp=solswp*swradcorr+sollwp
   endif
 
   !
   ! Ajouter la tendance des rayonnements (tous les pas)
+  ! avec une correction pour le cycle diurne dans le SW
   !