Changeset 5874 for LMDZ6/branches/PBLSURF_GPUPORT

Timestamp:

Nov 19, 2025, 11:03:49 AM (3 days ago)

Author:

yann meurdesoif

Message:

GPU port of pbl_surface_uncompress_pre

YM

File:

• LMDZ6/branches/PBLSURF_GPUPORT/libf/phylmd/pbl_surface_mod.F90 (modified) (10 diffs)

LMDZ6/branches/PBLSURF_GPUPORT/libf/phylmd/pbl_surface_mod.F90

@@ -324 +324 @@
 
   END SUBROUTINE pbl_surface_init_iso
-
+#endif
 !
 !****************************************************************************************
@@ -1494 +1494 @@
 
 
-
   SUBROUTINE pbl_surface_uncompress_pre( &
        itap,          &
@@ -1541 +1540 @@
 #endif      
      &   )
+!$gpum horizontal klon
+
 !****************************************************************************************
 ! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
@@ -1623 +1624 @@
     use hbtm_mod, only: hbtm
     USE indice_sol_mod
-    USE time_phylmdz_mod,   ONLY : 
     USE mod_grid_phy_lmdz,  ONLY :  grid1dto2d_glo
     USE print_control_mod,  ONLY : prt_level
@@ -1642 +1642 @@
     USE ioipsl_getin_p_mod, ONLY : getin_p
     use phys_state_var_mod, only:  frac_tersrf, albedo_tersrf !AM 
-    use phys_output_var_mod, only: 
-    use lmdz_blowing_snow_ini, only : 
     USE dimsoil_mod_h, ONLY: nsoilmx
     USE surf_param_mod, ONLY: eff_surf_param  !AM
@@ -1870 +1868 @@
     REAL, DIMENSION(klon, nbsrf), INTENT(OUT)       :: trmb3_w
 
-    REAL, PARAMETER                    :: facteur=2./sqrt(3.14)
+    REAL, PARAMETER                    :: facteur = 2. / 1.772  ! ( == 2. / SQRT(3.14))
     REAL, PARAMETER                    :: inertia=2000.
     REAL, DIMENSION(klon,nbsrf), INTENT(OUT)        :: pblh         ! height of the planetary boundary layer
@@ -2257 +2255 @@
 ! >> PC
    IF (carbon_cycle_cpl .AND. carbon_cycle_tr .AND. nbcf_out.GT.0 ) THEN
-       r_co2_ppm(:) = co2_send(:)
+       r_co2_ppm(1:klon) = co2_send(1:klon)
        DO n=1,nbcf_out
-           IF (cfname_out(n) == "atmco2" ) fields_out(:,n) = co2_send(:)
+           IF (cfname_out(n) == "atmco2" ) fields_out(1:klon,n) = co2_send(1:klon)
        ENDDO
    ENDIF
 
    IF ( .NOT. carbon_cycle_tr .AND. nbcf_out.GT.0 ) THEN
-       r_co2_ppm(:) = co2_ppm     ! Constant field
+       r_co2_ppm(1:klon) = co2_ppm     ! Constant field
        DO n=1,nbcf_out
-           IF (cfname_out(n) == "atmco2" ) fields_out(:,n) = co2_ppm
+           IF (cfname_out(n) == "atmco2" ) fields_out(1:klon,n) = co2_ppm
        ENDDO
    ENDIF
@@ -2880 +2878 @@
     REAL, DIMENSION(knon)              :: ydTs_ins, ydqs_ins
 !
-    REAL, PARAMETER                    :: facteur=2./sqrt(3.14)
+    REAL, PARAMETER                    :: facteur = 2. / 1.772  ! ( == 2. / SQRT(3.14))
     REAL, PARAMETER                    :: inertia=2000.
     REAL, DIMENSION(knon)              :: ydtsurf_th
@@ -5430 +5428 @@
     REAL, DIMENSION(klon, nbsrf), INTENT(IN)       :: pblh_w
     REAL, DIMENSION(klon, nbsrf), INTENT(IN)       :: plcl_w
-    REAL, PARAMETER                    :: facteur=2./sqrt(3.14)
+    REAL, PARAMETER                    :: facteur = 2. / 1.772  ! ( == 2. / SQRT(3.14))
     REAL, PARAMETER                    :: inertia=2000.
 
@@ -5612 +5610 @@
     zxtsol(:) = zxtsol(:)+1.5*meansqT(:)/zxtsol(:)
    ENDIF   ! iflag_order2_sollw == 1
-       
+
+!$gpum nocall       
        CALL checksum("n2mout", n2mout) 
        CALL checksum("n2mout_x", n2mout_x)