Changeset 1244 for LMDZ4

Timestamp:

Sep 17, 2009, 4:59:40 PM (15 years ago)

Author:

Laurent Fairhead

Message:

Le calcul de ratqs faisait osciller le modele et donnait une solution
asymptotique bidon et dépendant fortement du pas de temps.

Location:

LMDZ4/branches/LMDZ4-dev/libf/phylmd

Files:

• conf_phys.F90 (modified) (6 diffs)
• physiq.F (modified) (3 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/conf_phys.F90

@@ -11 +11 @@
                        fact_cldcon, facttemps,ok_newmicro,iflag_radia,&
                        iflag_cldcon, &
-                       iflag_ratqs,ratqsbas,ratqshaut, &
+                       iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
                        ok_ade, ok_aie, aerosol_couple, &
                        flag_aerosol, new_aod, &
@@ -65 +65 @@
   LOGICAL              :: new_aod
   REAL                 :: bl95_b0, bl95_b1
-  real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut
+  real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut,tau_ratqs
   integer              :: iflag_cldcon
   integer              :: iflag_ratqs
@@ -81 +81 @@
   real,SAVE           :: fact_cldcon_omp, facttemps_omp,ratqsbas_omp
   real,SAVE           :: ratqshaut_omp
+  real,SAVE           :: tau_ratqs_omp
   integer,SAVE        :: iflag_radia_omp
   integer,SAVE        :: iflag_rrtm_omp
@@ -729 +730 @@
   ratqshaut_omp = 0.3
   call getin('ratqshaut',ratqshaut_omp)
+
+!Config Key  = tau_ratqs
+!Config Desc =  
+!Config Def  = 1800.
+!Config Help = 
+!
+  tau_ratqs_omp = 1800.
+  call getin('tau_ratqs',tau_ratqs_omp)
 
 !
@@ -1375 +1384 @@
     ratqsbas = ratqsbas_omp
     ratqshaut = ratqshaut_omp
+    tau_ratqs = tau_ratqs_omp
+
     iflag_radia = iflag_radia_omp
     iflag_rrtm = iflag_rrtm_omp
@@ -1500 +1511 @@
   write(numout,*)' ratqsbas = ',ratqsbas 
   write(numout,*)' ratqshaut = ',ratqshaut 
+  write(numout,*)' tau_ratqs = ',tau_ratqs 
   write(numout,*)' top_height = ',top_height 
   write(numout,*)' overlap = ',overlap 

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

@@ -878 +878 @@
 c
       REAL ratqss(klon,klev),ratqsc(klon,klev)
-      real ratqsbas,ratqshaut
-      save ratqsbas,ratqshaut
-c$OMP THREADPRIVATE(ratqsbas,ratqshaut)
+      real ratqsbas,ratqshaut,tau_ratqs
+      save ratqsbas,ratqshaut,tau_ratqs
+c$OMP THREADPRIVATE(ratqsbas,ratqshaut,tau_ratqs)
       real zpt_conv(klon,klev)
 
@@ -1248 +1248 @@
      .                  solarlong0,seuil_inversion,
      .                  fact_cldcon, facttemps,ok_newmicro,iflag_radia,
-     .                  iflag_cldcon,iflag_ratqs,ratqsbas,ratqshaut,
+     .         iflag_cldcon,iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs,
      .                  ok_ade, ok_aie, aerosol_couple, 
      .                  flag_aerosol, new_aod,
@@ -2530 +2530 @@
          enddo
 
-!   les ratqs sont une conbinaison de ratqss et ratqsc
-!   1800s, en dur pour le moment, est le temps de
-!   relaxation des ratqs
-
-         facteur=exp(-pdtphys/1800.)
-
-         print*,'WARNING ratqs a revoir '
+!   les ratqs sont une combinaison de ratqss et ratqsc
+       print*,'PHYLMD NOUVEAU TAU_RATQS ',tau_ratqs
+
+         if (tau_ratqs>1.e-10) then
+            facteur=exp(-pdtphys/tau_ratqs)
+         else
+            facteur=0.
+         endif
          ratqs(:,:)=ratqsc(:,:)*(1.-facteur)+ratqs(:,:)*facteur
-         ratqs(:,:)=sqrt(ratqs(:,:)**2+ratqss(:,:)**2)
-
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! FH 22/09/2009
+! La ligne ci-dessous faisait osciller le modele et donnait une solution
+! assymptotique bidon et dépendant fortement du pas de temps.
+!        ratqs(:,:)=sqrt(ratqs(:,:)**2+ratqss(:,:)**2)
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+         ratqs(:,:)=max(ratqs(:,:),ratqss(:,:))
       else
 !   on ne prend que le ratqs stable pour fisrtilp