Changeset 4147 for trunk/LMDZ.COMMON/libf/evolution/output.F90

Timestamp:

Mar 23, 2026, 2:55:10 PM (3 weeks ago)

Author:

jbclement

Message:

PEM:

• Simplification of subroutines to convert data between the physical and the dynamical/lon-lat grids + making them more robust.
• Correction for air mass to give back to the PCM. The variable is extensive so poles must be treated specifically.
• Making the PEM able to do 0 year.
• Explicit information about the frost values computed by the PEM + enforcing positivity of yearly minima.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/output.F90 (modified) (4 diffs)

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

@@ -383 +383 @@
 call check_nc(nf90_enddef(ncid_diagevo),'leaving define mode')
 allocate(var_dyn(nlon_loc,nlat))
-call vect2dyngrd(nlon + 1,nlat,ngrid,longitudes,var_dyn)
+call vect2dyngrd(longitudes,var_dyn)
 call put_var_nc('longitude',var_dyn)
-call vect2dyngrd(nlon + 1,nlat,ngrid,latitudes,var_dyn)
+call vect2dyngrd(latitudes,var_dyn)
 call put_var_nc('latitude',var_dyn)
 call put_var_nc('ap',ap)
 call put_var_nc('bp',bp)
 call put_var_nc('soildepth',mlayer)
-call vect2dyngrd(nlon + 1,nlat,ngrid,cell_area,var_dyn)
+call vect2dyngrd(cell_area,var_dyn)
 call put_var_nc('cell_area',var_dyn)
 deallocate(var_dyn)
@@ -499 +499 @@
                 rank(1)
                     allocate(var_dyn(nlon_loc,nlat))
-                    call vect2dyngrd(nlon + 1,nlat,ngrid,var,var_dyn)
+                    call vect2dyngrd(var,var_dyn)
                     call put_var_nc(var_name,var_dyn,itime)
                     deallocate(var_dyn)
@@ -505 +505 @@
                     allocate(var1d_dyn(nlon_loc,nlat,size(var,2)))
                     do i = 1,size(var,2)
-                        call vect2dyngrd(nlon + 1,nlat,ngrid,var(:,i),var1d_dyn(:,:,i))
+                        call vect2dyngrd(var(:,i),var1d_dyn(:,:,i))
                     end do
                     call put_var_nc(var_name,var1d_dyn,itime)
@@ -513 +513 @@
                     do i = 1,size(var,2)
                         do j = 1,size(var,3)
-                            call vect2dyngrd(nlon + 1,nlat,ngrid,var(:,i,j),var2d_dyn(:,:,i,j))
+                            call vect2dyngrd(var(:,i,j),var2d_dyn(:,:,i,j))
                         end do
                     end do