Changeset 3827 for trunk/LMDZ.VENUS/libf/phyvenus/conf_phys.F90

Timestamp:

Jul 3, 2025, 3:28:27 PM (4 weeks ago)

Author:

ikovalenko

Message:

 

File:

• trunk/LMDZ.VENUS/libf/phyvenus/conf_phys.F90 (modified) (51 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/conf_phys.F90

@@ -9 +9 @@
    use init_print_control_mod, only: init_print_control
    use print_control_mod, only: lunout
-   use IOIPSL, only: getin
+   use ioipsl_getin_p_mod, only: getin_p
    use age_of_air_mod, only: ok_aoa, reinit_aoa, lev_aoa
 
@@ -44 +44 @@
        ! do we read a startphy.nc file? (default: .true.)
        startphy_file=.true.
-       CALL getin("startphy_file",startphy_file)
+       CALL getin_p("startphy_file",startphy_file)
 
 !Config  Key  = cycle_diurne
@@ -52 +52 @@
 !Config         Peut etre util pour accelerer le code !
        cycle_diurne = .TRUE.       
-       call getin('cycle_diurne',cycle_diurne)
+       call getin_p('cycle_diurne',cycle_diurne)
 
 !Config  Key  = soil_model
@@ -61 +61 @@
 !Config         plus de choix ! Ou meme une liste d'options !
        soil_model = .true.
-       call getin('soil_model',soil_model)
+       call getin_p('soil_model',soil_model)
 
 !Config  Key  = ok_orodr
@@ -69 +69 @@
 !Config         
        ok_orodr = .false.
-       call getin('ok_orodr',ok_orodr)
+       call getin_p('ok_orodr',ok_orodr)
 
 !Config  Key  =  ok_orolf
@@ -76 +76 @@
 !Config  Help = GW lift orographie (pas utilise)
        ok_orolf = .false.
-       call getin('ok_orolf', ok_orolf)
+       call getin_p('ok_orolf', ok_orolf)
 
 !Config  Key  = ok_gw_nonoro
@@ -83 +83 @@
 !Config  Help = GW drag non-orographique
        ok_gw_nonoro = .false.
-       call getin('ok_gw_nonoro',ok_gw_nonoro)
+       call getin_p('ok_gw_nonoro',ok_gw_nonoro)
 
 !Config  Key  = nbapp_rad
@@ -91 +91 @@
 !Config         par jour.
        nbapp_rad = 12
-       call getin('nbapp_rad',nbapp_rad)
+       call getin_p('nbapp_rad',nbapp_rad)
        print*,"nbapp_rad",nbapp_rad
 !Config  Key  = nbapp_chem
@@ -99 +99 @@
 !Config         par jour.
        nbapp_chem = 24000
-       call getin('nbapp_chem',nbapp_chem)
+       call getin_p('nbapp_chem',nbapp_chem)
 
 !Config  Key  = iflag_con
@@ -110 +110 @@
 !Config         3 pour CCM(NCAR)  
        iflag_con = 0
-       call getin('iflag_con',iflag_con)
+       call getin_p('iflag_con',iflag_con)
 
 !******************* fin parametres anciennement lus dans gcm.def
@@ -121 +121 @@
 !
   ok_journe = .false.
-  call getin('OK_journe', ok_journe)
+  call getin_p('OK_journe', ok_journe)
 !
 !Config Key  = OK_mensuel
@@ -130 +130 @@
 !
   ok_mensuel = .false.
-  call getin('OK_mensuel', ok_mensuel)
+  call getin_p('OK_mensuel', ok_mensuel)
 !
 !Config Key  = OK_instan
@@ -139 +139 @@
 !
   ok_instan = .false.
-  call getin('OK_instan', ok_instan)
+  call getin_p('OK_instan', ok_instan)
 !
 !Config  Key  = ecritphy
@@ -148 +148 @@
 !
        ecriphy = 1.
-       call getin('ecritphy', ecriphy)
+       call getin_p('ecritphy', ecriphy)
 !
 !
@@ -158 +158 @@
 !
   if_ebil = 0
-  call getin('if_ebil', if_ebil)
+  call getin_p('if_ebil', if_ebil)
 !!
 !! Parametres orbitaux
@@ -170 +170 @@
 ! R_ecc = 0.006787
   R_ecc   = 0.0
-  call getin('R_ecc', R_ecc)
+  call getin_p('R_ecc', R_ecc)
 !!
 !Config Key  = R_peri
@@ -179 +179 @@
 ! VENUS
   R_peri = 0.
-  call getin('R_peri', R_peri)
+  call getin_p('R_peri', R_peri)
 !!
 !Config Key  = R_incl
@@ -188 +188 @@
 ! VENUS
   R_incl = 0.0
-  call getin('R_incl', R_incl)
+  call getin_p('R_incl', R_incl)
 !
 !
@@ -206 +206 @@
 ! 6-9 = schema des thermiques Fred
   iflag_pbl = 1
-  call getin('iflag_pbl',iflag_pbl)
+  call getin_p('iflag_pbl',iflag_pbl)
 
 !
@@ -215 +215 @@
 !
   ksta = 1.0e-7
-  call getin('ksta',ksta)
+  call getin_p('ksta',ksta)
 
 !
@@ -224 +224 @@
 !
   z0 = 1.0e-2
-  call getin('z0',z0)
+  call getin_p('z0',z0)
 
 !
@@ -233 +233 @@
 !
   lmixmin = 35.
-  call getin('lmixmin',lmixmin)
+  call getin_p('lmixmin',lmixmin)
 
 !
@@ -242 +242 @@
 !
   ok_kzmin = .false.
-  call getin('ok_kzmin',ok_kzmin)
+  call getin_p('ok_kzmin',ok_kzmin)
 
   ok_clmain = .true.
-  call getin('ok_clmain',ok_clmain)
+  call getin_p('ok_clmain',ok_clmain)
 
   physideal = .false.
-  call getin('physideal',physideal)
+  call getin_p('physideal',physideal)
 
 !Config Key  = iflag_ajs
@@ -256 +256 @@
 !
   iflag_ajs = 0
-  call getin('iflag_ajs',iflag_ajs)
+  call getin_p('iflag_ajs',iflag_ajs)
 
 !
@@ -265 +265 @@
 !
   inertie = 2000.
-  call getin('inertie',inertie)
+  call getin_p('inertie',inertie)
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -277 +277 @@
 !
   lev_histins = 0
-  call getin('lev_histins',lev_histins)
+  call getin_p('lev_histins',lev_histins)
 
 !
@@ -286 +286 @@
 !
   lev_histday = 1
-  call getin('lev_histday',lev_histday)
+  call getin_p('lev_histday',lev_histday)
 
 !
@@ -295 +295 @@
 !
   lev_histmth = 2
-  call getin('lev_histmth',lev_histmth)
+  call getin_p('lev_histmth',lev_histmth)
 
 
@@ -317 +317 @@
 !       Need ok_chem or ok_cloud
   tr_scheme = 0
-  call getin('tr_scheme',tr_scheme)
+  call getin_p('tr_scheme',tr_scheme)
 
 !Config Key = ok_aoa
@@ -323 +323 @@
 !Config Def = .false.
   ok_aoa=.false.
-  call getin('ok_aoa',ok_aoa)
+  call getin_p('ok_aoa',ok_aoa)
 
 !Config Key = reinit_aoa
@@ -329 +329 @@
 !Config Def = .true.
   reinit_aoa=.true.
-  call getin('reinit_aoa',reinit_aoa)
+  call getin_p('reinit_aoa',reinit_aoa)
   
 !Config Key = lev_aoa
@@ -335 +335 @@
 !Config Def = 1 (surface)
   lev_aoa=1
-  call getin('lev_aoa',lev_aoa)
+  call getin_p('lev_aoa',lev_aoa)
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -348 +348 @@
 !
   reinit_trac = .FALSE.
-  call getin('reinit_trac',reinit_trac)
+  call getin_p('reinit_trac',reinit_trac)
   
 !
@@ -357 +357 @@
 !
   ok_cloud = .false.
-  call getin('ok_cloud',ok_cloud)
+  call getin_p('ok_cloud',ok_cloud)
 
 !
@@ -369 +369 @@
 
   cl_scheme = 2
-  call getin('cl_scheme',cl_scheme)
+  call getin_p('cl_scheme',cl_scheme)
 
 !
@@ -378 +378 @@
 !
   ok_chem = .false.
-  call getin('ok_chem',ok_chem)
+  call getin_p('ok_chem',ok_chem)
 
   if (((tr_scheme.ne.3).and.(ok_chem.or.ok_cloud)).or. &
@@ -406 +406 @@
 !
   ok_sedim = .false.
-  call getin('ok_sedim',ok_sedim)
+  call getin_p('ok_sedim',ok_sedim)
 
 !
@@ -415 +415 @@
 !
   nb_mode = 0
-  call getin('nb_mode',nb_mode)
+  call getin_p('nb_mode',nb_mode)
   
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -430 +430 @@
 ! 2 = Generic module
   solarchoice = 1
-  call getin('solarchoice',solarchoice)
+  call getin_p('solarchoice',solarchoice)
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -443 +443 @@
 !
   callnlte = .false.
-  call getin('callnlte',callnlte)
+  call getin_p('callnlte',callnlte)
 
 !
@@ -451 +451 @@
 !
   callnirco2 = .false.
-  call getin('callnirco2',callnirco2)
+  call getin_p('callnirco2',callnirco2)
 
 !
@@ -460 +460 @@
 !
   nircorr = 0
-  call getin('nircorr',nircorr)
+  call getin_p('nircorr',nircorr)
 
 !
@@ -469 +469 @@
 !
   callthermos = .false.
-  call getin('callthermos',callthermos)
+  call getin_p('callthermos',callthermos)
 
 !
@@ -478 +478 @@
 !
   nltemodel = 0
-  call getin('nltemodel',nltemodel)
+  call getin_p('nltemodel',nltemodel)
 
 !
@@ -487 +487 @@
 !
   solvarmod = 1
-  call getin('solvarmod',solvarmod)
+  call getin_p('solvarmod',solvarmod)
 
 !
@@ -496 +496 @@
 !
   fixed_euv_value =140. 
-  call getin('fixed_euv_value',fixed_euv_value)
+  call getin_p('fixed_euv_value',fixed_euv_value)
 
 !
@@ -505 +505 @@
 !
   euveff = 0.22
-  call getin('euveff',euveff)
+  call getin_p('euveff',euveff)
 
 !Config Key  = tuneupperatm
@@ -513 +513 @@
 !
   tuneupperatm = .false.
-  call getin('tuneupperatm',tuneupperatm)
+  call getin_p('tuneupperatm',tuneupperatm)
 
 !
@@ -522 +522 @@
 !
   ok_jonline = .false.
-  call getin('ok_jonline',ok_jonline)
+  call getin_p('ok_jonline',ok_jonline)
 
 !
@@ -531 +531 @@
 !
   ok_ionchem = .false.
-  call getin('ok_ionchem',ok_ionchem)
+  call getin_p('ok_ionchem',ok_ionchem)
 
   if ((ok_jonline.eqv..false.).and.(.true..eqv.ok_ionchem)) then
@@ -549 +549 @@
 !
   ok_iondiff = .false.
-  call getin('ok_iondiff',ok_iondiff)
+  call getin_p('ok_iondiff',ok_iondiff)
 
   if ((ok_ionchem.eqv..false.).and.(.true..eqv.ok_iondiff)) then