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

Timestamp:

Jan 17, 2025, 5:14:18 PM (2 weeks ago)

Author:

jbclement

Message:

PEM:

• Albedo is now updated only at the end of the PEM (and not at every iteration) + Correct way to set it taking into account CO2/H2O ice and frost.
• Cosmetic cleanings.

JBC

File:

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

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

@@ -42 +42 @@
 
 do numyear = 1,2
-    write(*,*) 'numyear',numyear
-    write(str(1:1),'(i1.1)') numyear
-
-    state = nf90_open(path = "data2reshape_Y"//str//".nc", mode = nf90_nowrite, ncid = ncid1)
+    write(str,'(i1.1)') numyear
+    write(*,*) 'Reshaping of variables from "data2reshape_Y'//str//'.nc"...'
+
+    state = nf90_open(path = "data2reshape_Y"//str//".nc",mode = nf90_nowrite,ncid = ncid1)
     if (state /= nf90_noerr) call handle_err(state)
 
@@ -57 +57 @@
         endif
     endif
-    state = nf90_create(path = "data_PCM_Y"//str//".nc", cmode=or(nf90_noclobber,nf90_64bit_offset), ncid = ncid2)
+    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)
 
@@ -69 +69 @@
     allocate(varids_2(nvars))
 
-    state = nf90_inq_dimids(ncid1, ndims, dimids, include_parents)
-    if (state /= nf90_noerr) call handle_err(state)
-    state = nf90_inq_varids(ncid1, nvars, varids)
+    state = nf90_inq_dimids(ncid1,ndims,dimids,include_parents)
+    if (state /= nf90_noerr) call handle_err(state)
+    state = nf90_inq_varids(ncid1,nvars,varids)
     if (state /= nf90_noerr) call handle_err(state)
 
     do i = 1,ndims
-        state = nf90_inquire_dimension(ncid1, dimids(i), name_, len_)
+        state = nf90_inquire_dimension(ncid1,dimids(i),name_,len_)
         if (state /= nf90_noerr) call handle_err(state)
         if (name_ == "lon" .or. name_ == "longitude") then
@@ -82 +82 @@
             len_ = len_ + 1
         else if (name_ == "lat".or. name_ == "latitude") then
-            dimid_lat=dimids(i)
-            len_lat=len_
+            dimid_lat = dimids(i)
+            len_lat = len_
         else if (name_ == "time_counter".or. name_ ==  "Time") then
-            dimid_time=dimids(i)
-            len_time=len_
+            dimid_time = dimids(i)
+            len_time = len_
         else if (name_ == "soil_layers".or. name_ ==  "subsurface_layers") then
-            dimid_soil=dimids(i)
+            dimid_soil = dimids(i)
             len_soil = len_
         endif
-        state = nf90_def_dim(ncid2, name_,len_,dimid_2)
+        state = nf90_def_dim(ncid2,name_,len_,dimid_2)
         if (state /= nf90_noerr) call handle_err(state)
         dimids_2(i) = dimid_2
@@ -97 +97 @@
 
     do i = 1,nvars
-        state = nf90_inquire_variable(ncid1, varids(i),name = namevar,xtype = xtype_var,ndims = numdims,natts = numatts)
-        write(*,*) "namevar00= ", namevar
+        state = nf90_inquire_variable(ncid1,varids(i),name = namevar,xtype = xtype_var,ndims = numdims,natts = numatts)
+        write(*,*) '> Treatment of '//namevar
         if (state /= nf90_noerr) call handle_err(state)
         allocate(dimid_var(numdims))
@@ -109 +109 @@
                 state = nf90_get_var(ncid1,varids(i),tempvalues_1d)
                 if (state /= nf90_noerr) call handle_err(state)
-                state = nf90_def_var(ncid2,namevar,xtype_var, dimid_var, varids_2(i))
+                state = nf90_def_var(ncid2,namevar,xtype_var,dimid_var,varids_2(i))
                 if (state /= nf90_noerr) call handle_err(state)
                 values_1d(1:len_lon) = tempvalues_1d(:)
@@ -219 +219 @@
 
     deallocate(dimids,varids,dimids_2,varids_2)
+    write(*,*) 'Done!'
 enddo