Changeset 4110 for LMDZ6/trunk/libf/phylmd/dyn1d

Timestamp:

Mar 25, 2022, 10:30:13 PM (3 years ago)

Author:

fhourdin

Message:

Pour la capacité "replay" des paramétrisations

Location:

LMDZ6/trunk/libf/phylmd/dyn1d

Files:

• old_lmdz1d.F90 (modified) (2 diffs)
• replay1d.F90 (modified) (5 diffs)

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

@@ -1033 +1033 @@
 !
 !=====================================================================
+#ifdef OUTPUT_PHYS_SCM
        CALL iophys_ini(timestep)
+#endif
 ! START OF THE TEMPORAL LOOP :
 !=====================================================================
@@ -1286 +1288 @@
      &             +dt_cooling(1:mxcalc))  ! Taux de chauffage ou refroid.
 
-print*,'OLDLMDZ1D IOPH'
+#ifdef OUTPUT_PHYS_SCM
       CALL iophys_ecrit('d_t_adv',klev,'d_t_adv','m/s',d_t_adv)
       CALL iophys_ecrit('d_t_nudge',klev,'d_t_nudge','m/s',d_t_nudge)
+#endif
 
       endif  ! forcing_sandu or forcing_astex

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

@@ -3 +3 @@
   USE mod_const_mpi, ONLY: comm_lmdz
   USE inigeomphy_mod, ONLY: inigeomphy
+ USE comvert_mod, ONLY: presnivs
+USE comvert_mod, only :  preff, pa
+
 
 
@@ -8 +11 @@
 #include "dimensions.h"
 
-real presnivs(llm)
 real :: airefi
 real :: zcufi    = 1.
@@ -27 +29 @@
   rlat_rad(1)=0.
   rlon_rad(1)=0.
+
+preff=101325.
+pa=50000.
+  CALL disvert()
   CALL inigeomphy(1,1,llm, &
                1, comm_lmdz, &
@@ -36 +42 @@
 
 CALL suphel
-
-!---------------------------------------------------------------------
-! Initialisation d'un fichier de sorties
-!---------------------------------------------------------------------
-call system("ncdump -h phys.nc  |grep -i time | head -1 | cut -d\( -f2  | awk ' { print $1 } ' > presnivs")
-call system("ncks -v lev phys.nc | grep 'lev.*.,' | sed -e 's/lev//' -e 's/[=;,]//g' >> presnivs")
-open(10,file="presnivs",form="formatted")
-read(10,*) ntime
-read(10,*) presnivs
-close(10)
+ntime=100
 jour0=1
 mois0=1
@@ -51 +48 @@
 klon=1
 klev=llm
-call iotd_ini('physb.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.,900.,'360d')
 
 !---------------------------------------------------------------------