Changeset 3106 for trunk/LMDZ.COMMON

Timestamp:

Oct 26, 2023, 11:30:35 AM (13 months ago)

Author:

jbclement

Message:

PEM:

• Correction of a bug when running with ifort, XIOS and no sub-slopes. If compiled with ifort, the 'Reshape' program cashed because of sub-slopes variables were written in the XIOS output but not filled. A new file definition without sub-slopes can be chosen in the xml file to avoid this issue.
• Addition in the 'Reshape' program of a check to remove the output files if they exist to prevent a crash.

JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• reshape_XIOS_output.F90 (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/evolution/changelog.txt

@@ -120 +120 @@
 Following r3095, 'iniorbit' is now a subroutine of "planete_h.F90".
 Correction of a bug from r3088: "diagfi.nc" was not changed into "diagpem.nc" in the commited version of "writediagpem.F90"!
+
+== 26/10/2023 == JBC
+    - Correction of a bug when running with ifort, XIOS and no sub-slopes. If compiled with ifort, the 'Reshape' program cashed because of sub-slopes variables were written in the XIOS output but not filled. A new file definition without sub-slopes can be chosen in the xml file to avoid this issue.
+    - Addition in the 'Reshape' program of a check to remove the output files if they exist to prevent a crash.

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

@@ -15 +15 @@
 
 integer                               :: state, ncid, ncid1, ncid2, nDims, nVars, nGlobalAtts, unlimDimID
-integer                               :: i, j, include_parents
+integer                               :: i, j, include_parents, cstat
 integer, dimension(:),    allocatable :: dimids, varids, dimids_2, varids_2, dimid_var
 real, dimension(:),       allocatable :: tempvalues_1d, values_1d
@@ -21 +21 @@
 real, dimension(:,:,:),   allocatable :: tempvalues_3d, values_3d
 real, dimension(:,:,:,:), allocatable :: tempvalues_4d, values_4d
-character(1)                          :: str2
+character(1)                          :: str
 character(30)                         :: name_, namevar
 integer                               :: xtype_var, len_, len_1, len_2, len_lat, len_lon, len_time, len_soil
 integer                               :: dimid_lon, dimid_lat, dimid_time, dimid_soil, dimid_2, numdims, numatts, numyear
+logical                               :: yes
 
 do numyear = 1,2
     write(*,*) 'numyear',numyear
-    write(str2(1:1),'(i1.1)') numyear
+    write(str(1:1),'(i1.1)') numyear
     !nf90_open                 ! open existing netCDF dataset
     !integer :: ncid, state
     !...
-    state = nf90_open(path = "data2reshape"//str2//".nc", mode = nf90_nowrite, ncid = ncid1)
-    if (state /= nf90_noerr) call handle_err(state)
-
-    state = nf90_create(path = "data_PCM_Y"//str2//".nc", cmode=or(nf90_noclobber,nf90_64bit_offset), ncid = ncid2)
+    state = nf90_open(path = "data2reshape"//str//".nc", mode = nf90_nowrite, ncid = ncid1)
+    if (state /= nf90_noerr) call handle_err(state)
+
+    inquire(file = 'data_PCM_Y'//str//'.nc', exist = yes)
+    if (yes) then
+        call execute_command_line('rm data_PCM_Y'//str//'.nc',cmdstat = cstat)
+        if (cstat > 0) then
+            error stop 'Command exection failed!'
+        else if (cstat < 0) then
+            error stop 'Command execution not supported!'
+        endif
+    endif
+    state = nf90_create(path = "data_PCM_Y"//str//".nc", cmode=or(nf90_noclobber,nf90_64bit_offset), ncid = ncid2)
     if (state /= nf90_noerr) call handle_err(state)