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

Timestamp:

Jan 31, 2008, 6:20:42 PM (17 years ago)

Author:

Laurent Fairhead

Message:

modifications pour faire de l'aquaplanète FH
LF

File:

• LMDZ4/trunk/libf/phylmd/conf_phys.F90 (modified) (6 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, &
+ &                     solarlong0,qsol0,seuil_inversion, &
  &                     fact_cldcon, facttemps,ok_newmicro,iflag_radia,&
  &                     iflag_cldcon, &
@@ -26 +26 @@
 include "compbl.h"
 include "control.h"
+include "comsoil.h"
 !
 ! Configuration de la "physique" de LMDZ a l'aide de la fonction
@@ -93 +94 @@
   INTEGER,SAVE :: iflag_pdf_omp
   REAL,SAVE :: rad_froid_omp, rad_chau1_omp, rad_chau2_omp
+  REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_ice_omp
+  REAL      :: qsol0
+  REAL,SAVE :: qsol0_omp
+  REAL      :: solarlong0
+  REAL,SAVE :: solarlong0_omp
   INTEGER,SAVE :: top_height_omp,overlap_omp
   REAL,SAVE :: cdmmax_omp,cdhmax_omp,ksta_omp,ksta_ter_omp
@@ -536 +542 @@
   ratqshaut_omp = 0.3
   call getin('ratqshaut',ratqshaut_omp)
+
+!
+!-----------------------------------------------------------------------
+! Longitude solaire pour le calcul de l'ensoleillement en degre
+! si on veut imposer la saison. Sinon, solarlong0=-999.999
+!Config Key  = solarlong0
+!Config Desc =  
+!Config Def  = -999.999 
+!Config Help = 
+!
+  solarlong0_omp = -999.999
+  call getin('solarlong0',solarlong0_omp)
+!
+!-----------------------------------------------------------------------
+!  Valeur imposee de l'humidite du sol pour le modele bucket.
+!Config Key  = qsol0
+!Config Desc =  
+!Config Def  = -1.
+!Config Help = 
+!
+  qsol0_omp = -1.
+  call getin('qsol0',qsol0_omp)
+!
+!-----------------------------------------------------------------------
+!
+!Config Key  = inertie_ice
+!Config Desc =  
+!Config Def  = 2000.
+!Config Help = 
+!
+  inertie_ice_omp = 2000.
+  call getin('inertie_ice',inertie_ice_omp)
+!
+!Config Key  = inertie_sno
+!Config Desc =  
+!Config Def  = 2000.
+!Config Help = 
+!
+  inertie_sno_omp = 2000.
+  call getin('inertie_sno',inertie_sno_omp)
+!
+!Config Key  = inertie_sol
+!Config Desc =  
+!Config Def  = 2000.
+!Config Help = 
+!
+  inertie_sol_omp = 2000.
+  call getin('inertie_sol',inertie_sol_omp)
 
 !
@@ -872 +926 @@
     reevap_ice = reevap_ice_omp
     iflag_pdf = iflag_pdf_omp
+    solarlong0 = solarlong0_omp
+    qsol0 = qsol0_omp
+    inertie_sol = inertie_sol_omp
+    inertie_ice = inertie_ice_omp
+    inertie_sno = inertie_sno_omp
     rad_froid = rad_froid_omp
     rad_chau1 = rad_chau1_omp
@@ -996 +1055 @@
   write(numout,*)' type_run = ',type_run 
   write(numout,*)' ok_isccp = ',ok_isccp 
+  WRITE(numout,*)' solarlong0 = ', solarlong0
+  write(numout,*)' qsol0 = ', qsol0
+  write(numout,*)' inertie_sol = ', inertie_sol
+  write(numout,*)' inertie_ice = ', inertie_ice
+  write(numout,*)' inertie_sno = ', inertie_sno
+
   write(numout,*)' lonmin lonmax latmin latmax bilKP_ins =',&
  & lonmin_ins, lonmax_ins, latmin_ins, latmax_ins