Changeset 878 for LMDZ4/trunk/libf/phylmd/conf_phys.F90

Timestamp:

Jan 14, 2008, 1:03:39 PM (16 years ago)

Author:

Laurent Fairhead

Message:

Bascule de la physique du LMD vers la physique avec thermiques
LF

File:

• LMDZ4/trunk/libf/phylmd/conf_phys.F90 (modified) (9 diffs)

LMDZ4/trunk/libf/phylmd/conf_phys.F90

@@ -6 +6 @@
 
   subroutine conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, ok_instan, ok_hf, &
+ &                     seuil_inversion, &
  &                     fact_cldcon, facttemps,ok_newmicro,iflag_cldcon, &
-!IM&                   ratqsbas,ratqshaut,ip_ebil_phy, &
- &                     ratqsbas,ratqshaut, &
+ &                     iflag_ratqs,ratqsbas,ratqshaut, &
  &                     ok_ade, ok_aie, &
  &                     bl95_b0, bl95_b1,&
@@ -51 +51 @@
   real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut
   integer              :: iflag_cldcon
+  integer              :: iflag_ratqs
 
   character (len = 6),SAVE  :: ocean_omp
@@ -61 +62 @@
   real,SAVE           :: ratqshaut_omp
   integer,SAVE        :: iflag_cldcon_omp, ip_ebil_phy_omp
+  integer,SAVE        :: iflag_ratqs_omp
 
 ! Local
   integer              :: numout = 6
   real                 :: zzz
+
+  real :: seuil_inversion
+  real,save :: seuil_inversion_omp
 
   integer :: iflag_thermals,nsplit_thermals
@@ -82 +87 @@
   REAL,SAVE :: cdmmax_omp,cdhmax_omp,ksta_omp,ksta_ter_omp
   LOGICAL,SAVE :: ok_kzmin_omp
-  REAL, SAVE ::  fmagic_omp
+  REAL, SAVE :: fmagic_omp
   INTEGER,SAVE :: iflag_pbl_omp,lev_histhf_omp,lev_histday_omp,lev_histmth_omp
   CHARACTER*4, SAVE :: type_run_omp
@@ -207 +212 @@
 !Config Help = 
 !               
-!
   ip_ebil_phy_omp = 0
   call getin('ip_ebil_phy', ip_ebil_phy_omp)
+!
+!Config Key  = seuil_inversion
+!Config Desc = Seuil ur dTh pour le choix entre les schemas de CL
+!Config Def  = -0.1
+!Config Help = 
+!               
+  seuil_inversion_omp = -0.1
+  call getin('seuil_inversion', seuil_inversion_omp)
+
 !!
 !! Constante solaire & Parametres orbitaux & taux gaz effet de serre BEG
@@ -438 +451 @@
   reevap_ice_omp = .false.
   call getin('reevap_ice',reevap_ice_omp)
+
+!Config Key  = iflag_ratqs
+!Config Desc =
+!Config Def  = 1
+!Config Help =
+!
+  iflag_ratqs_omp = 1
+  call getin('iflag_ratqs',iflag_ratqs_omp)
+
 !
 !Config Key  = iflag_cldcon 
@@ -835 +857 @@
     ratqshaut = ratqshaut_omp
     iflag_cldcon = iflag_cldcon_omp
+    iflag_ratqs = iflag_ratqs_omp
     ip_ebil_phy = ip_ebil_phy_omp
     iflag_thermals = iflag_thermals_omp
@@ -840 +863 @@
     type_run = type_run_omp
     ok_isccp = ok_isccp_omp
+    seuil_inversion=seuil_inversion_omp
     lonmin_ins = lonmin_ins_omp
     lonmax_ins = lonmax_ins_omp
@@ -891 +915 @@
   write(numout,*)' iflag_pdf = ', iflag_pdf
   write(numout,*)' iflag_cldcon = ', iflag_cldcon
+  write(numout,*)' iflag_ratqs = ', iflag_ratqs
+  write(numout,*)' seuil_inversion = ', seuil_inversion
   write(numout,*)' fact_cldcon = ', fact_cldcon
   write(numout,*)' facttemps = ', facttemps