Changeset 5977 for LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90

Timestamp:

Dec 19, 2025, 2:22:27 PM (12 days ago)

Author:

yann meurdesoif

Message:

For better GPUmorphosis parsing use "PUBLIC :: var" instead of "PUBLIC var"
merge of commit r5923
YM

File:

• LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90 (modified) (9 diffs)

LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -1034 +1034 @@
     REAL, DIMENSION(klon, klev, nbsrf) :: flux_u_x, flux_v_x, flux_u_w, flux_v_w
     REAL, DIMENSION(klon, nbsrf)       :: fluxlat_x, fluxlat_w
-
+    
 
 
@@ -1042 +1042 @@
     INTEGER :: knon
     INTEGER :: ni(klon)
-
+    REAL,  DIMENSION(:, :), ALLOCATABLE ::  ytemp00, ytemp01, ytemp02, ytemp03, ytemp04, ytemp05, ytemp06, ytemp07, ytemp08, ytemp09
+    REAL,  DIMENSION(:, :), ALLOCATABLE ::  ytemp10, ytemp11, ytemp12, ytemp13, ytemp14, ytemp15, ytemp16, ytemp17, ytemp18, ytemp19
+
+    
       IF (is_master) WRITE(lunout,*) "****************** CHECKSUM IN ****************************"
       CALL checksum("dtime", dtime)
@@ -1268 +1271 @@
       ENDIF
     ENDDO
+
+    ALLOCATE(ytemp00(knon,klev+1), ytemp01(knon,klev+1), ytemp02(knon,klev+1), ytemp03(knon,klev+1), ytemp04(knon,klev+1), ytemp05(knon,klev+1), ytemp06(knon,klev+1), ytemp07(knon,klev+1), ytemp08(knon,klev+1), ytemp09(knon,klev+1))
+    ALLOCATE(ytemp10(knon,klev+1), ytemp11(knon,klev+1), ytemp12(knon,klev+1), ytemp13(knon,klev+1), ytemp14(knon,klev+1), ytemp15(knon,klev+1), ytemp16(knon,klev+1), ytemp17(knon,klev+1), ytemp18(knon,klev+1), ytemp19(knon,klev+1))
+
     
     CALL pbl_surface_subsrf( nsrf, knon, ni(1:knon),  &
@@ -1321 +1328 @@
        pblt_w, therm_w, trmb1_w, trmb2_w, trmb3_w, pblh, plcl, capCL, oliqCL, cteiCL, pblT, &
        therm, trmb1, trmb2, trmb3, alb, snowerosion, iflag_split_ref, delp, d_t_diss_x, d_t_diss_w, flux_t_x, flux_q_x, flux_t_w, flux_q_w,&
-       flux_u_x, flux_v_x, flux_u_w, flux_v_w, fluxlat_x, fluxlat_w)
+       flux_u_x, flux_v_x, flux_u_w, flux_v_w, fluxlat_x, fluxlat_w, &
+       ytemp00, ytemp01, ytemp02, ytemp03, ytemp04, ytemp05, ytemp06, ytemp07, ytemp08, ytemp09,& 
+       ytemp10, ytemp11, ytemp12, ytemp13, ytemp14, ytemp15, ytemp16, ytemp17, ytemp18, ytemp19)
+
   
   ENDDO
+  
+  DEALLOCATE(ytemp00, ytemp01, ytemp02, ytemp03, ytemp04, ytemp05, ytemp06, ytemp07, ytemp08, ytemp09)
+  DEALLOCATE(ytemp10, ytemp11, ytemp12, ytemp13, ytemp14, ytemp15, ytemp16, ytemp17, ytemp18, ytemp19)
 
   CALL pbl_surface_uncompressed_post( &
@@ -2356 +2369 @@
        therm, trmb1, trmb2, trmb3, alb, snowerosion, iflag_split_ref, &
        delp, d_t_diss_x, d_t_diss_w, flux_t_x, flux_q_x, flux_t_w, flux_q_w,&
-       flux_u_x, flux_v_x, flux_u_w, flux_v_w, fluxlat_x, fluxlat_w)
-!$gpum horizontal knon
+       flux_u_x, flux_v_x, flux_u_w, flux_v_w, fluxlat_x, fluxlat_w, &
+       ytemp00, ytemp01, ytemp02, ytemp03, ytemp04, ytemp05, ytemp06, ytemp07, ytemp08, ytemp09,&
+       ytemp10, ytemp11, ytemp12, ytemp13, ytemp14, ytemp15, ytemp16, ytemp17, ytemp18, ytemp19)
+!$gpum horizontal knon klon
 
 
@@ -2954 +2969 @@
     REAL, DIMENSION(knon) :: yzxtsol     ! temperature at surface
     REAL, DIMENSION(knon)                   :: ypblh_tmp ! temporaire pblh compressed 
+    REAL, DIMENSION(knon) :: zgeo1_over_RG
 #ifdef ISO
     REAL, DIMENSION(knon)       :: h1
     INTEGER                     :: ixt
 #endif
+
+    REAL, INTENT(OUT), DIMENSION(knon, klev+1) ::  ytemp00, ytemp01, ytemp02, ytemp03, ytemp04, ytemp05, ytemp06, ytemp07, ytemp08, ytemp09
+    REAL, INTENT(OUT), DIMENSION(knon, klev+1) ::  ytemp10, ytemp11, ytemp12, ytemp13, ytemp14, ytemp15, ytemp16, ytemp17, ytemp18, ytemp19
+   
+   ytemp00=0 ; ytemp01=0 ; ytemp02=0 ; ytemp03=0 ; ytemp04=0 ; ytemp05=0 ; ytemp06=0 ; ytemp07=0 ; ytemp08=0 ; ytemp09=0
+   ytemp10=0 ; ytemp11=0 ; ytemp12=0 ; ytemp13=0 ; ytemp14=0 ; ytemp15=0 ; ytemp16=0 ; ytemp17=0 ; ytemp18=0 ; ytemp19=0
+
+!$gpum nocall checksum
+
       IF (using_xios) THEN
         missing_val=missing_val_xios
@@ -3312 +3337 @@
             ENDDO
             !
-
-            CALL eff_surf_param(knon, nbtersrf, yz0m_tersrf, yfrac_tersrf, 'CDN', yz0m, zgeo1/RG)
-            CALL eff_surf_param(knon, nbtersrf, yz0h_tersrf, yfrac_tersrf, 'CDN', yz0h, zgeo1/RG)
+            zgeo1_over_RG(1:knon) = zgeo1(1:knon)/RG
+            CALL eff_surf_param(knon, nbtersrf, yz0m_tersrf, yfrac_tersrf, 'CDN', yz0m, zgeo1_over_RG)
+            CALL eff_surf_param(knon, nbtersrf, yz0h_tersrf, yfrac_tersrf, 'CDN', yz0h, zgeo1_over_RG)
             !
           ENDIF
@@ -3990 +4015 @@
 
         end if
-
+           zgeo1_over_RG(1:knon) = zgeo1(1:knon)/RG
            CALL surf_ocean(rlon, rlat, ysolsw, ysollw, yalb_vis, &
                ywindsp, yrmu0, yfder, yts, &
                itap, dtime, jour, knon, ni, &
-               ypplay(:,1), zgeo1(1:knon)/RG, ycdragh, ycdragm, yrain_f, ysnow_f, ybs_f, yt(:,1), yq(:,1),&    ! ym missing init
+               ypplay(:,1), zgeo1_over_RG, ycdragh, ycdragm, yrain_f, ysnow_f, ybs_f, yt(:,1), yq(:,1),&    ! ym missing init
                AcoefH, AcoefQ, BcoefH, BcoefQ, &
                AcoefU, AcoefV, BcoefU, BcoefV, &
@@ -4011 +4036 @@
 #endif               
          &      )
-!$gpum nocall
+
            CALL checksum("yalb_dir_new_ocean",yalb_dir_new(1:knon,:))
       IF (prt_level >=10) THEN