Changeset 4350 for LMDZ6/trunk/libf/phylmd

Timestamp:

Nov 18, 2022, 2:41:55 PM (2 years ago)

Author:

fhourdin

Message:

Details pour la fonctionalie Replay

Location:

LMDZ6/trunk/libf/phylmd

Files:

• dyn1d/replay1d.F90 (modified) (4 diffs)
• iophys.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/dyn1d/replay1d.F90

@@ -1 +1 @@
 PROGRAM rejouer
 
-  USE mod_const_mpi, ONLY: comm_lmdz
-  USE inigeomphy_mod, ONLY: inigeomphy
- USE comvert_mod, ONLY: presnivs
+USE mod_const_mpi, ONLY: comm_lmdz
+USE inigeomphy_mod, ONLY: inigeomphy
+USE comvert_mod, ONLY: presnivs
 USE comvert_mod, only :  preff, pa
+USE ioipsl, only: getin
+
 
 
@@ -17 +19 @@
 
 integer ntime
-integer jour0,mois0,an0
-integer it
+integer jour0,mois0,an0,day_step,anneeref,dayref
 integer klev,klon
+CHARACTER (len=10) :: calend
+CHARACTER(len=20) :: calendrier
+
 
 !---------------------------------------------------------------------
@@ -31 +35 @@
 
 preff=101325.
+!preff=100000.
 pa=50000.
-open(82,file='dump_param.bin',form='unformatted',status='old')
   CALL disvert()
   CALL inigeomphy(1,1,llm, &
@@ -43 +47 @@
 
 CALL suphel
-ntime=100
-jour0=1
-mois0=1
-an0=2000
+!ntime=4320
+ntime=10000000
+dayref=1
+anneeref=2000
+call getin('dayref',dayref)
+call getin('anneeref',anneeref)
+call getin('calend',calend)
+call getin('day_step',day_step)
+calendrier=calend
+if ( calendrier == "earth_360d" ) calendrier="360d"
+
+
+jour0=dayref
+mois0=(jour0-1)/30+1
+jour0=jour0-30*((jour0-1)/30)
+an0=anneeref
+
+!print*,"REPLAY1D jour0,mois0,an0",jour0,mois0,an0
+
+
 klon=1
 klev=llm
-call iotd_ini('phys.nc',1,1,klev,0.,0.,presnivs,jour0,mois0,an0,0.,900.,'360d')
+call iotd_ini('phys.nc',1,1,klev,0.,0.,presnivs,jour0,mois0,an0,0.,86400./day_step,calendrier)
+! Consistent with ... CALL iophys_ini(600.)
 
 !---------------------------------------------------------------------
 ! Initialisation de la parametrisation
 !---------------------------------------------------------------------
-      call get_ini_module
+call get_ini_module
 
 !---------------------------------------------------------------------
 ! Boucle en temps sur l'appel à la parametrisation
 !---------------------------------------------------------------------
-      DO it=1,ntime
-         print*,'Pas de temps ',it,klon,klev
-         call call_param_replay(klon,klev)
-      ENDDO
-      END
+call call_param_replay(klon,klev)
 
+end
 
 !---------------------------------------------------------------------

LMDZ6/trunk/libf/phylmd/iophys.F90

@@ -170 +170 @@
 mois0=day_ref/30+1
 an0=annee_ref
-t0=(day_ini-1)*RDAY
+!FH BIZARE QUAND 1D ...  t0=(day_ini-1)*RDAY
+t0=0.
 calendrier=calend
-
 if ( calendrier == "earth_360d" ) calendrier="360d"
+
+print*,'iophys_ini annee_ref day_ref',annee_ref,day_ref,day_ini,calend,t0
+
 
 call iotd_ini('phys.nc', &