Changeset 1928 for LMDZ4/branches

Timestamp:

Jan 16, 2014, 1:43:34 PM (11 years ago)

Author:

idelkadi

Message:

• Correction du bug dans l'interface avec la physique : On ne veut pas que la distinction entre les nuages convectifs et stratiformes soit prise en compte dans les calculs Cosp. Dans l'interface avec la physique, sont passees en entree pour Cosp, les quantites totales (stratiforme + convective) du contenus en eau et d'autres variables. La fraction nuageuse convective calculee dans la physique est passee en entree pour Cosp dans la version buggee. Elle est remise a 0 dans cette version corrigee.

• Mise a jour pour ISCCP : la fraction d'ensoleillement calculee par LMDZ est passee en argument pour Isccp

• Rajouts d'autre variables diagnostiques : sunlit (=0 si nuit et =1 si jour) parasol_crefl (reflectance integree)

Location:

LMDZ4/branches/LMDZ4_AR5/libf

Files:

• cosp/ini_histdayCOSP.h (modified) (4 diffs)
• cosp/ini_histhfCOSP.h (modified) (4 diffs)
• cosp/ini_histmthCOSP.h (modified) (3 diffs)
• cosp/phys_cosp.F90 (modified) (6 diffs)
• cosp/write_histdayCOSP.h (modified) (2 diffs)
• cosp/write_histhfCOSP.h (modified) (2 diffs)
• cosp/write_histmthCOSP.h (modified) (2 diffs)
• phylmd/physiq.F (modified) (1 diff)

LMDZ4/branches/LMDZ4_AR5/libf/cosp/ini_histdayCOSP.h

@@ -21 +21 @@
         CALL histvert(nid_day_cosp,"height","height","m",Nlevout,vgrid%z,nvert)
        else
-        CALL histvert(nid_day_cosp,"presnivs","Vertical levels","mb",Nlevout,presnivs,nvert)
+        CALL histvert(nid_day_cosp,"presnivs","Vertical levels","mb",Nlevout,presnivs,nvert,"down")
        endif
 
@@ -28 +28 @@
         CALL histvert(nid_day_cosp,"sza","solar_zenith_angle","degrees",PARASOL_NREFL,PARASOL_SZA,nvertp)
 
-        CALL histvert(nid_day_cosp,"pressure2","pressure","mb",7,ISCCP_PC,nvertisccp)
+        CALL histvert(nid_day_cosp,"pressure2","pressure","mb",7,ISCCP_PC,nvertisccp,"down")
 
         CALL histvert(nid_day_cosp,"column","column","count",Ncolumns,column_ax(1:Ncolumns),nvertcol)
@@ -77 +77 @@
                         iim,jj_nb,nhori, PARASOL_NREFL,1, PARASOL_NREFL, nvertp,32, &
                         "ave(X)", zout,zstoday)   
+            CALL histdef(nid_day_cosp, "parasol_crefl", &
+                        "PARASOL-like mono-directional reflectance","1", &
+                        iim,jj_nb,nhori, PARASOL_NREFL,1, PARASOL_NREFL, nvertp,32, &
+                        "ave(X)", zout,zstoday)
+            CALL histdef(nid_day_cosp, "Ncrefl", &
+                        "0 ou 1 si reflectance <= ou > 0","1", &
+                        iim,jj_nb,nhori, PARASOL_NREFL,1, PARASOL_NREFL, nvertp,32, &
+                        "ave(X)", zout,zstoday)
            endif
            if (cfg%Latb532) then
@@ -134 +142 @@
 ! Sorties ISCCP
         if (cfg%Lisccp_sim) then
+          CALL histdef(nid_day_cosp, "sunlit", &
+                     "1 si jour 0 si nuit", "1", &
+                     iim, jj_nb,nhori,1,1,1,-99,32, &
+                     "ave(X)", zout,zstomth)
          if (cfg%Lclisccp2) then
             do ii=1,7

LMDZ4/branches/LMDZ4_AR5/libf/cosp/ini_histhfCOSP.h

@@ -21 +21 @@
         CALL histvert(nid_hf_cosp,"height","height","m",Nlevout,vgrid%z,nvert)
        else
-        CALL histvert(nid_hf_cosp,"presnivs","Vertical levels","mb",Nlevout,presnivs,nvert)
+        CALL histvert(nid_hf_cosp,"presnivs","Vertical levels","mb",Nlevout,presnivs,nvert,"down")
        endif
 
@@ -28 +28 @@
         CALL histvert(nid_hf_cosp,"sza","solar_zenith_angle","degrees",PARASOL_NREFL,PARASOL_SZA,nvertp)
 
-        CALL histvert(nid_hf_cosp,"pressure2","pressure","mb",7,ISCCP_PC,nvertisccp)
+        CALL histvert(nid_hf_cosp,"pressure2","pressure","mb",7,ISCCP_PC,nvertisccp,"down")
 
         CALL histvert(nid_hf_cosp,"column","column","count",Ncolumns,column_ax(1:Ncolumns),nvertcol)
@@ -77 +77 @@
                         iim,jj_nb,nhori, PARASOL_NREFL,1, PARASOL_NREFL, nvertp,32, &
                         "ave(X)", zout,zstohf)   
+            CALL histdef(nid_hf_cosp, "parasol_crefl", &
+                        "PARASOL-like mono-directional reflectance","1", &
+                        iim,jj_nb,nhori, PARASOL_NREFL,1, PARASOL_NREFL, nvertp,32, &
+                        "ave(X)", zout,zstohf)
+            CALL histdef(nid_hf_cosp, "Ncrefl", &
+                        "1 ou 0 si cldt > o <= 0 ","1", &
+                        iim,jj_nb,nhori, PARASOL_NREFL,1, PARASOL_NREFL, nvertp,32, &
+                        "ave(X)", zout,zstohf)
            endif
            if (cfg%Latb532) then
@@ -134 +142 @@
 ! Sorties ISCCP
         if (cfg%Lisccp_sim) then
+
+          CALL histdef(nid_hf_cosp, "sunlit", &
+                     "0 si nuit 1 si jour", "1", &
+                     iim, jj_nb,nhori,1,1,1,-99,32, &
+                     "ave(X)", zout,zstomth)
          if (cfg%Lclisccp2) then
             do ii=1,7

LMDZ4/branches/LMDZ4_AR5/libf/cosp/ini_histmthCOSP.h

@@ -27 +27 @@
         CALL histvert(nid_mth_cosp,"sza","solar_zenith_angle","degrees",PARASOL_NREFL,PARASOL_SZA,nvertp)
 
-        CALL histvert(nid_mth_cosp,"pressure2","pressure","mb",7,ISCCP_PC,nvertisccp)
+        CALL histvert(nid_mth_cosp,"pressure2","pressure","mb",7,ISCCP_PC,nvertisccp,"down")
 
         CALL histvert(nid_mth_cosp,"column","column","count",Ncolumns,column_ax(1:Ncolumns),nvertcol)
@@ -84 +84 @@
                         iim,jj_nb,nhori, PARASOL_NREFL,1, PARASOL_NREFL, nvertp,32, &
                         "ave(X)", zout,zstomth)   
+            CALL histdef(nid_mth_cosp, "parasol_crefl", &
+                        "PARASOL-like mono-directional reflectance","1", &
+                        iim,jj_nb,nhori, PARASOL_NREFL,1, PARASOL_NREFL, nvertp,32, &
+                        "ave(X)", zout,zstomth)
+            CALL histdef(nid_mth_cosp, "Ncrefl", &
+                        "0 ou 1 si reflectance <= ou > 0","1", &
+                        iim,jj_nb,nhori, PARASOL_NREFL,1, PARASOL_NREFL, nvertp,32, &
+                        "ave(X)", zout,zstomth)
            endif
            if (cfg%Latb532) then
@@ -141 +149 @@
 ! Sorties ISCCP
         if (cfg%Lisccp_sim) then
+
+             CALL histdef(nid_mth_cosp, "sunlit", &
+                     "1 si jour 0 si nuit", "1", &
+                     iim, jj_nb,nhori,1,1,1,-99,32, &
+                     "ave(X)", zout,zstomth)
          if (cfg%Lclisccp2) then
             do ii=1,7

LMDZ4/branches/LMDZ4_AR5/libf/cosp/phys_cosp.F90

@@ -1 +1 @@
 ! Simulateur COSP : Cfmip Observation Simulator Package
+
 ! ISCCP, Radar (QuickBeam), Lidar et Parasol (ACTSIM), MISR, RTTOVS
 !Idelkadi Abderrahmane Aout-Septembre 2009
-
 
   subroutine phys_cosp( itap,dtime,freq_cosp, &
                         ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP, &
                         ecrit_mth,ecrit_day,ecrit_hf, &
-                        Nptslmdz,Nlevlmdz,lon,lat, presnivs,overlaplmdz, &
+                        Nptslmdz,Nlevlmdz,lon,lat, presnivs,overlaplmdz,sunlit, &
                         ref_liq,ref_ice,fracTerLic,u_wind,v_wind,phis,phi,ph,p,skt,t, &
                         sh,rh,tca,cca,mr_lsliq,mr_lsice,fl_lsrainI,fl_lssnowI, &
@@ -144 +144 @@
                                      zlev,zlev_half,mr_ozone,radliq,radice,dtau_s,dem_s,ref_liq,ref_ice
   real,dimension(Nptslmdz,Nlevlmdz) ::  fl_lsrainI,fl_lssnowI,fl_ccrainI,fl_ccsnowI
-  real,dimension(Nptslmdz)        :: lon,lat,skt,fracTerLic,u_wind,v_wind,phis
+  real,dimension(Nptslmdz)        :: lon,lat,skt,fracTerLic,u_wind,v_wind,phis,sunlit    
   real,dimension(Nlevlmdz)        :: presnivs
   integer                         :: itap,k,ip
   real                            :: dtime,freq_cosp
   logical, parameter              :: lCOSP=.FALSE.
+
+  real, dimension(Nptslmdz,PARASOL_NREFL) :: parasolcrefl, Ncref
 
 !
@@ -237 +239 @@
         gbx%q = rh*100.
         gbx%sh = sh
-        gbx%cca = cca !convective_cloud_amount (1)
+        gbx%cca = 0.
         gbx%tca = tca ! total_cloud_amount (1)
         gbx%psfc = ph(:,1) !pression de surface
@@ -253 +255 @@
         gbx%u_wind  = u_wind !eastward_wind (m s-1)
         gbx%v_wind  = v_wind !northward_wind
-! Attention
-        gbx%sunlit  = 1
+
+      do ip = 1, Npoints
+        if (sunlit(ip).le.0.) then
+           gbx%sunlit(ip)=0
+        else
+           gbx%sunlit(ip)=1
+        endif
+      enddo
 
 ! A voir l equivalent LMDZ
@@ -280 +288 @@
      gbx%Reff(:,:,I_CVCLIQ) = ref_liq*1e-6
      gbx%Reff(:,:,I_CVCICE) = ref_ice*1e-6
-
         ! ISCCP simulator
         gbx%dtau_s   = dtau_s
@@ -423 +430 @@
      if(isccp%boxptop(ip,ii).eq.-1.E+30)then
        isccp%boxptop(ip,ii)=0.
+     endif
+    enddo
+   enddo
+
+   do k=1,PARASOL_NREFL
+    do ii=1, Npoints
+     if (stlidar%cldlayer(ii,4).gt.0.01) then
+        parasolcrefl(ii,k)=(stlidar%parasolrefl(ii,k)-0.03*(1.-stlidar%cldlayer(ii,4)))/ &
+                   stlidar%cldlayer(ii,4)
+        Ncref(ii,k) = 1.
+     else
+        parasolcrefl(ii,k)=0.
+        Ncref(ii,k) = 0.
      endif
     enddo

LMDZ4/branches/LMDZ4_AR5/libf/cosp/write_histdayCOSP.h

@@ -38 +38 @@
            CALL histwrite_phy(nid_day_cosp,lCOSP,"parasol_refl", &
      itau_wcosp,stlidar%parasolrefl)
+           CALL histwrite_phy(nid_day_cosp,lCOSP,"parasol_crefl", &
+     itau_wcosp,parasolcrefl)
+           CALL histwrite_phy(nid_day_cosp,lCOSP,"Ncrefl", &
+     itau_wcosp,Ncref)
          endif
          if (cfg%Latb532) then
@@ -91 +95 @@
 ! Sorties ISCCP
         if (cfg%Lisccp_sim) then
+          CALL histwrite_phy(nid_day_cosp,lCOSP,"sunlit", &
+     itau_wcosp,sunlit)
           if (cfg%Lclisccp2) then
             do ii=1,7

LMDZ4/branches/LMDZ4_AR5/libf/cosp/write_histhfCOSP.h

@@ -38 +38 @@
            CALL histwrite_phy(nid_hf_cosp,lCOSP,"parasol_refl", &
      itau_wcosp,stlidar%parasolrefl)
+           CALL histwrite_phy(nid_hf_cosp,lCOSP,"parasol_crefl", &
+     itau_wcosp,parasolcrefl)
+           CALL histwrite_phy(nid_hf_cosp,lCOSP,"Ncrefl", &
+     itau_wcosp,Ncref)
          endif
          if (cfg%Latb532) then
@@ -91 +95 @@
 ! Sorties ISCCP
         if (cfg%Lisccp_sim) then
+
+          CALL histwrite_phy(nid_hf_cosp,lCOSP,"sunlit", &
+                        itau_wcosp,gbx%sunlit)
           if (cfg%Lclisccp2) then
             do ii=1,7

LMDZ4/branches/LMDZ4_AR5/libf/cosp/write_histmthCOSP.h

@@ -38 +38 @@
            CALL histwrite_phy(nid_mth_cosp,lCOSP,"parasol_refl", &
      itau_wcosp,stlidar%parasolrefl)
+           CALL histwrite_phy(nid_mth_cosp,lCOSP,"parasol_crefl", &
+     itau_wcosp,parasolcrefl)
+           CALL histwrite_phy(nid_mth_cosp,lCOSP,"Ncrefl", &
+     itau_wcosp,Ncref)
          endif
          if (cfg%Latb532) then
@@ -91 +95 @@
 ! Sorties ISCCP
         if (cfg%Lisccp_sim) then
+       CALL histwrite_phy(nid_mth_cosp,lCOSP,"sunlit",itau_wcosp,gbx%sunlit)
           if (cfg%Lclisccp2) then
             do ii=1,7

LMDZ4/branches/LMDZ4_AR5/libf/phylmd/physiq.F

@@ -3447 +3447 @@
      $                   ecrit_mth,ecrit_day,ecrit_hf,
      $                   klon,klev,rlon,rlat,presnivs,overlap,
-     $                   ref_liq,ref_ice,
+     $                   fract,ref_liq,ref_ice,
      $                   pctsrf(:,is_ter)+pctsrf(:,is_lic),
      $                   zu10m,zv10m,pphis,