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

Timestamp:

Mar 23, 2026, 2:55:10 PM (2 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/surf_temp.F90 (modified) (2 diffs)

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

@@ -169 +169 @@
 call print_msg("> Adapting surface temperature where ice disappeared",LVL_NFO)
 ! Convert from reduced grid to lon-lat grid
-call vect2lonlat(nlon,nlat,ngrid,latitudes,latitude_ll)
+call vect2lonlat(latitudes,latitude_ll)
 do islope = 1,nslope
-    call vect2lonlat(nlon,nlat,ngrid,tsurf_avg(:,islope),tsurf_ll(:,:,islope))
-    call vect2lonlat(nlon,nlat,ngrid,surfice(:,islope),co2ice_ll(:,:,islope))
-    call vect2lonlat(nlon,nlat,ngrid,merge(1._dp,0._dp,is_ice_ini(:,islope)),mask_co2ice_ini(:,:,islope))
-    call vect2lonlat(nlon,nlat,ngrid,merge(1._dp,0._dp,ice_disappeared(:,islope)),mask_co2ice_disappeared(:,:,islope))
+    call vect2lonlat(tsurf_avg(:,islope),tsurf_ll(:,:,islope))
+    call vect2lonlat(surfice(:,islope),co2ice_ll(:,:,islope))
+    call vect2lonlat(merge(1._dp,0._dp,is_ice_ini(:,islope)),mask_co2ice_ini(:,:,islope))
+    call vect2lonlat(merge(1._dp,0._dp,ice_disappeared(:,islope)),mask_co2ice_disappeared(:,:,islope))
 end do
 
@@ -233 +233 @@
 ! Convert back from lon-lat grid to reduced grid
 do islope = 1,nslope
-    call lonlat2vect(nlon,nlat,ngrid,tsurf_ll(:,:,islope),tsurf_avg(:,islope))
-    call lonlat2vect(nlon,nlat,ngrid,mask_co2ice_disappeared(:,:,islope),tmp)
+    call lonlat2vect(tsurf_ll(:,:,islope),tsurf_avg(:,islope))
+    call lonlat2vect(mask_co2ice_disappeared(:,:,islope),tmp)
     where (tmp > 0.5_dp) ice_disappeared(:,islope) = .true.
 end do