Changeset 3214 for trunk/LMDZ.COMMON/libf/evolution/reshape_XIOS_output.F90

Timestamp:

Feb 12, 2024, 5:37:25 PM (9 months ago)

Author:

jbclement

Message:

PEM:

• It is now possible to set the number of initial PCM calls independently of the number of "inter-PEM" PCM calls. It is useful to get a stable situation for the start of the simulation.
• Correction of a bug: 'reshape_XIOS_output' treated the first two PCM runs instead of the last two. The numbering has been adapted in "launch_pem.sh" to get it right.
• PEM outputs ("diagpem.nc" and "diagsoil_pem.nc") has been improved: 'Time' now evolves according to 'dt_pem' (usually year by year) and 'ecritpem' is a full variable of "time_evol_mod.F90".

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/reshape_XIOS_output.F90 (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/evolution/reshape_XIOS_output.F90

@@ -14 +14 @@
 implicit none
 
-integer                               :: state, ncid, ncid1, ncid2, nDims, nVars, nGlobalAtts, unlimDimID
-integer                               :: i, j, include_parents, cstat
+integer                               :: state, ncid1, ncid2, nDims, nVars, nGlobalAtts, unlimDimID
+integer                               :: i, include_parents, cstat
 integer, dimension(:),    allocatable :: dimids, varids, dimids_2, varids_2, dimid_var
 real, dimension(:),       allocatable :: tempvalues_1d, values_1d
@@ -31 +31 @@
     write(str(1:1),'(i1.1)') numyear
 
-    state = nf90_open(path = "data2reshape"//str//".nc", mode = nf90_nowrite, ncid = ncid1)
+    state = nf90_open(path = "data2reshape_Y"//str//".nc", mode = nf90_nowrite, ncid = ncid1)
     if (state /= nf90_noerr) call handle_err(state)