Changeset 1928 for LMDZ4/branches/LMDZ4_AR5/libf/cosp/ini_histdayCOSP.h

Timestamp:

Jan 16, 2014, 1:43:34 PM (10 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)

File:

• LMDZ4/branches/LMDZ4_AR5/libf/cosp/ini_histdayCOSP.h (modified) (4 diffs)

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