Changeset 674

Timestamp:

Nov 16, 2005, 3:47:19 PM (19 years ago)

Author:

lmdzadmin

Message:

Initialisations - AC
MAF

Location:

LMDZ4/trunk/libf/phylmd

Files:

• clmain.F (modified) (2 diffs)
• hbtm.F (modified) (1 diff)
• physiq.F (modified) (4 diffs)

LMDZ4/trunk/libf/phylmd/clmain.F

@@ -316 +316 @@
       LOGICAL check
       PARAMETER (check=.false.)
+
+
+c initialisation Anne
+      ytherm(:) = 0.
 C
       if (check) THEN
@@ -500 +504 @@
         yts(j) = ts(i,nsrf)
 cIM "slab" ocean
+c        PRINT *, 'tslab = ', i, tslab(i)
         ytslab(i) = tslab(i)
 c

LMDZ4/trunk/libf/phylmd/hbtm.F

@@ -168 +168 @@
 #include "YOETHF.h"
 #include "FCTTRE.h"
-       b212=sqrt(b1*b2)
-       b2sr=sqrt(b2)
+
+
+
+! initialisations (Anne)
+      th_th(:) = 0.
+      q_star = 0
+      t_star = 0
+
+
+      b212=sqrt(b1*b2)
+      b2sr=sqrt(b2)
 c
 C ============================================================

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -125 +125 @@
 cIM "slab" ocean
       REAL tslab(klon)    !Temperature du slab-ocean
+      SAVE tslab
       REAL seaice(klon)   !glace de mer (kg/m2) 
+      SAVE seaice
       REAL fluxo(klon)    !flux turbulents ocean-glace de mer 
       REAL fluxg(klon)    !flux turbulents ocean-atmosphere
@@ -1222 +1224 @@
       SAVE ok_ade, ok_aie, bl95_b0, bl95_b1
 cym
+c Anne
+         SAVE u10m
+         SAVE v10m
+         SAVE t2m
+         SAVE q2m
+         SAVE ffonte
+         SAVE fqcalving
+         SAVE piz_ae
+         SAVE tau_ae
+         SAVE cg_ae
+         SAVE rain_con
+         SAVE snow_con
+         SAVE bl95_b0
+         SAVE bl95_b1
+         SAVE topswai
+         SAVE topswad
+         SAVE solswai
+         SAVE solswad
+         SAVE d_u_con
+         SAVE d_v_con
+         SAVE rnebcon0
+         SAVE clwcon0
+         SAVE rnebcon
+         SAVE clwcon
+         SAVE paire_ter
+         SAVE nhistoW
+         SAVE histoW
+c SAVE anne 20/09/2005
+         SAVE pblh
+         SAVE plcl
+         SAVE capCL
+         SAVE oliqCL
+         SAVE cteiCL
+         SAVE pblt
+         SAVE therm
+         SAVE trmb1
+         SAVE trmb2
+         SAVE trmb3
+
+c fin Anne
 cjq-end
 c
@@ -1284 +1326 @@
          histoW(:,:,:,:) = 0.0
 ! fin anne
+! Anne 12/09/2005
+
+         pblh(:,:)   =0.        ! Hauteur de couche limite
+         plcl(:,:)   =0.        ! Niveau de condensation de la CLA
+         capCL(:,:)  =0.        ! CAPE de couche limite
+         oliqCL(:,:) =0.        ! eau_liqu integree de couche limite
+         cteiCL(:,:) =0.        ! cloud top instab. crit. couche limite
+         pblt(:,:)   =0.        ! T a la Hauteur de couche limite
+         therm(:,:)  =0.
+         trmb1(:,:)  =0.        ! deep_cape
+         trmb2(:,:)  =0.        ! inhibition 
+         trmb3(:,:)  =0.        ! Point Omega
+! fin Anne
 
 cym
@@ -1787 +1842 @@
 
       fder = dlw
-
 
       CALL clmain(dtime,itap,date0,pctsrf,pctsrf_new,