Changeset 1528 for LMDZ4/branches/LMDZ4_AR5/libf/cosp/ini_histmthCOSP.h

Timestamp:

May 26, 2011, 2:06:28 PM (13 years ago)

Author:

idelkadi

Message:

Les modificaions concerne la partie Lidar du simulateur COSP :

1. Correction dans phys_cosp.F90 du bug sur le calcul des altitudes aux milieux et aux interfaces des couches
2. Rajout du teste "Lidar parfait" (iflag_testlidar=.true. dans lidar_simulator.F90)
3. Correction d'une erreur sur la dimension du tableau frac_out dans cosp_lidar.F90
4. Sories sur les niveaux presnivs du modele dans le cas ou l'on desactive l'interpolation sur les 40 niveaux COSP

File:

• LMDZ4/branches/LMDZ4_AR5/libf/cosp/ini_histmthCOSP.h (modified) (3 diffs)

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

@@ -14 +14 @@
         idayref = day_ref
         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
-
         CALL histbeg_phy("histmthCOSP",itau_phy,zjulian,dtime,nhori,nid_mth_cosp) 
 
 ! Definition de l'axe vertical
+        print*,'height_mlev=',vgrid%mz
+       if (use_vgrid) then 
         CALL histvert(nid_mth_cosp,"height","height","m",Nlevout,vgrid%z,nvert)
+       else
+        CALL histvert(nid_mth_cosp,"presnivs","Vertical levels","Pa",Nlevout,presnivs,nvert,"down")
+       endif
         CALL histvert(nid_mth_cosp,"height_mlev","height_mlev","m",Nlevlmdz,vgrid%mz,nvertm)
-!        CALL histvert(nid_mth_cosp,"presnivs","Vertical levels","mb",Nlevout,presnivs,nvert)
 
         CALL histvert(nid_mth_cosp,"sza","solar_zenith_angle","degrees",PARASOL_NREFL,PARASOL_SZA,nvertp)
@@ -26 +29 @@
         CALL histvert(nid_mth_cosp,"pressure2","pressure","mb",7,ISCCP_PC,nvertisccp)
 
-        CALL histvert(nid_mth_cosp,"column","column","count",Ncolumns,column_ax,nvertcol)
+        CALL histvert(nid_mth_cosp,"column","column","count",Ncolumns,column_ax(1:Ncolumns),nvertcol)
 
 ! Sorties LIDAR
@@ -59 +62 @@
                      iim,jj_nb,nhori, Nlevout,1,Nlevout,nvert, 32, &
                      "ave(X)", zout,zstomth)
+!         CALL histdef(nid_mth_cosp, "zlev_m", &
+!                     "Height at Midel model levels", "m", &
+!                     iim,jj_nb,nhori, Nlevlmdz,1,Nlevlmdz,nvertm, 32, &
+!                     "ave(X)", zout,zstomth)
+!         CALL histdef(nid_mth_cosp, "zlev_b", &
+!                     "Height at bottom model levels", "m", &
+!                     iim,jj_nb,nhori, Nlevlmdz,1,Nlevlmdz,nvertm, 32, &
+!                     "ave(X)", zout,zstomth)
          endif
            if (cfg%Lcfad_lidarsr532) then