Changeset 1258 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F

Timestamp:

Nov 9, 2009, 4:24:36 PM (15 years ago)

Author:

idelkadi

Message:

-Rajout des champs des precipitations verticales lsc (liq et ice) et conv (liq et ice)
-Appel dans physiq.F de l'interface avec le simulateur COSP sous la cles CPP_COSP avec une frequence=freq_COSP(s)
-Declaration de freq_COSP dans clesphys.h et lecture dans physiq.def via conf_phys.F90
-Sorie dans newmicro.F des rayons effectifs des particules liq et ice des nuages necessaires pour cosp

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F (modified) (4 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F

@@ -888 +888 @@
       logical ok_newmicro
       save ok_newmicro
+      real ref_liq(klon,klev), ref_ice(klon,klev)
 c$OMP THREADPRIVATE(ok_newmicro)
       save fact_cldcon,facttemps
@@ -1115 +1116 @@
 cIM 100106 END : pouvoir sortir les ctes de la physique
 c
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+c Declarations pour Simulateur COSP
+c============================================================
+      real :: mr_ozone(klon,klev)
 c======================================================================
 ! Ecriture eventuelle d'un profil verticale en entree de la physique.
@@ -2875 +2880 @@
      e            mass_solu_aero, mass_solu_aero_pi,
      e            bl95_b0, bl95_b1,
-     s            cldtaupi, re, fl)
+     s            cldtaupi, re, fl, ref_liq, ref_ice)
       else
       CALL nuage (paprs, pplay,
@@ -3168 +3173 @@
 c
 c
+!====================================================================
+! Interface Simulateur COSP (Calipso, ISCCP, MISR, ..)
+!====================================================================
+! Abderrahmane 24.08.09
+
+!      IF (ok_cosp) THEN
+! adeclarer 
+#ifdef CPP_COSP
+       IF (MOD(itap,NINT(freq_cosp/dtime)).EQ.0) THEN
+
+       print*,'freq_cosp',freq_cosp
+          mr_ozone=wo * dobson_u * 1e3 / zmasse
+!       print*,'Dans physiq.F avant appel cosp ref_liq,ref_ice=',
+!     s        ref_liq,ref_ice
+          call phys_cosp(itap,dtime,freq_cosp,
+     $                 ecrit_mth,ecrit_day,ecrit_hf,overlap,
+     $                   klon,klev,rlon,rlat,presnivs,
+     $                   ref_liq,ref_ice,
+     $                   pctsrf(:,is_ter)+pctsrf(:,is_lic),
+     $                   zu10m,zv10m,
+     $                   zphi,paprs(:,1:klev),pplay,zxtsol,t_seri,
+     $                   qx(:,:,ivap),zx_rh,cldfra,rnebcon,flwc,fiwc,
+     $                   prfl(:,1:klev),psfl(:,1:klev),
+     $                   pmflxr(:,1:klev),pmflxs(:,1:klev),
+     $                   mr_ozone,cldtau, cldemi)
+!     L          calipso2D,calipso3D,cfadlidar,parasolrefl,atb,betamol,
+!     L          cfaddbze,clcalipso2,dbze,cltlidarradar,
+!     M          clMISR,
+!     R          clisccp2,boxtauisccp,boxptopisccp,tclisccp,ctpisccp,
+!     I          tauisccp,albisccp,meantbisccp,meantbclrisccp)
+
+        ENDIF
+
+!      ENDIF
+#endif
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 cAA
 cAA Installation de l'interface online-offline pour traceurs