Changeset 3907 for trunk/LMDZ.COMMON/libf/evolution/pem.F90

Timestamp:

Sep 1, 2025, 6:45:25 PM (4 months ago)

Author:

jbclement

Message:

PEM:
Correction for the update of surface temperature when ice is disappearing: now it is really set as the surface temperature of the nearest bare soil point because we check the lon-lat grid and not the reduced grid.
JBC

File:

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

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

@@ -75 +75 @@
 use job_timelimit_mod,          only: timelimit, antetime, timewall
 use paleoclimate_mod,           only: h2o_ice_depth, zdqsdif_ssi_tot
+use surf_temp,                  only: update_tsurf_nearest_baresoil
 
 #ifndef CPP_STD
@@ -834 +835 @@
 !------------------------
 ! II_d.1 Update Tsurf
-    write(*,*) "> Updating surface temperature"
-    do ig = 1,ngrid
-        do islope = 1,nslope
-            ! CO2 ice disappeared so we look for the closest point without CO2 ice
-            if (is_co2ice_ini(ig,islope) .and. co2_ice(ig,islope) < 1.e-10 .and. .not. co2ice_disappeared(ig,islope)) then
-                co2ice_disappeared(ig,islope) = .true.
-                if (latitude_deg(ig) > 0.) then ! North hemisphere
-                    outer1: do ig_loop = ig,ngrid ! Go towards equator
-                        do islope_loop = islope - 1,1,-1 ! Go over the slopes (backward numbering - should be equator-ward)
-                            if (.not. is_co2ice_ini(ig_loop,islope_loop) .and. co2_ice(ig_loop,islope_loop) < 1.e-10) then
-                                tsurf_avg(ig,islope) = tsurf_avg(ig_loop,islope_loop)
-                                exit outer1
-                            endif
-                        enddo
-                        do islope_loop = islope + 1,nslope ! Go over the slopes (forward numbering - should be pole-ward)
-                            if (.not. is_co2ice_ini(ig_loop,islope_loop) .and. co2_ice(ig_loop,islope_loop) < 1.e-10) then
-                                tsurf_avg(ig,islope) = tsurf_avg(ig_loop,islope_loop)
-                                exit outer1
-                            endif
-                        enddo
-                    enddo outer1
-                else ! South hemisphere
-                    outer2: do ig_loop = ig,1,-1 ! Go towards equator
-                        do islope_loop = islope + 1,nslope ! Go over the slopes (forward numbering - should be equator-ward)
-                            if (.not. is_co2ice_ini(ig_loop,islope_loop) .and. co2_ice(ig_loop,islope_loop) < 1.e-10) then
-                                tsurf_avg(ig,islope) = tsurf_avg(ig_loop,islope_loop)
-                                exit outer2
-                            endif
-                        enddo
-                        do islope_loop = islope - 1,1,-1 ! Go over the slopes (backward numbering - should be pole-ward)
-                            if (.not. is_co2ice_ini(ig_loop,islope_loop) .and. co2_ice(ig_loop,islope_loop) < 1.e-10) then
-                                tsurf_avg(ig,islope) = tsurf_avg(ig_loop,islope_loop)
-                                exit outer2
-                            endif
-                        enddo
-                    enddo outer2
-                endif
-            else if (co2_ice(ig,islope) > 1.e-10 .and. d_co2ice(ig,islope) > 1.e-10) then ! Put tsurf as Tcond CO2
-                call computeTcondCO2(timelen,ngrid,nslope,vmr_co2_PEM_phys,ps_timeseries,ps_avg_global_ini,ps_avg_global_new,tsurf_avg)
-            endif
-        enddo
-    enddo
+    call update_tsurf_nearest_baresoil(ngrid,nslope,iim,jjm,tsurf_avg,co2_ice,is_co2ice_ini,co2ice_disappeared)
 
     if (soil_pem) then