Changeset 5512 for LMDZ6/trunk/libf/phylmd/lmdz_thermcell_flux2.f90

Timestamp:

Jan 28, 2025, 7:07:51 PM (2 days ago)

Author:

yann meurdesoif

Message:

Implement GPU automatic port for :

• Thermics
• acama_gwd_rando
• flott_gwd_rando

YM

File:

• LMDZ6/trunk/libf/phylmd/lmdz_thermcell_flux2.f90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_thermcell_flux2.f90

@@ -16 +16 @@
 !---------------------------------------------------------------------------
 
-      USE lmdz_thermcell_ini, ONLY : prt_level,iflag_thermals_optflux
+      USE lmdz_thermcell_ini, ONLY : prt_level,iflag_thermals_optflux, thermals_fomass_max, thermals_alphamax 
       IMPLICIT NONE
       
@@ -48 +48 @@
       REAL f_old,ddd0,eee0,ddd,eee,zzz
 
-      REAL,SAVE :: fomass_max=0.5
-      REAL,SAVE :: alphamax=0.7
-!$OMP THREADPRIVATE(fomass_max,alphamax)
-
       logical check_debug,labort_physic
 
-      character (len=20) :: modname='thermcell_flux2'
+      character (len=20), PARAMETER :: modname='thermcell_flux2'
       character (len=80) :: abort_message
 
@@ -391 +387 @@
         do ig=1,ngrid
            if (zw2(ig,l+1).gt.1.e-10) then
-           zfm=rhobarz(ig,l+1)*zw2(ig,l+1)*alphamax
+           zfm=rhobarz(ig,l+1)*zw2(ig,l+1)*thermals_alphamax
            if ( fm(ig,l+1) .gt. zfm) then
               f_old=fm(ig,l+1)
@@ -430 +426 @@
             eee0=entr(ig,l)
             ddd0=detr(ig,l)
-            eee=entr(ig,l)-masse(ig,l)*fomass_max/ptimestep
+            eee=entr(ig,l)-masse(ig,l)*thermals_fomass_max/ptimestep
             ddd=detr(ig,l)-eee
             if (eee.gt.0.) then
@@ -470 +466 @@
                          print*,'detr',detr(ig,l)
                          print*,'masse',masse(ig,l)
-                         print*,'fomass_max',fomass_max
-                         print*,'masse(ig,l)*fomass_max/ptimestep',masse(ig,l)*fomass_max/ptimestep
+                         print*,'thermal_fomass_max',thermals_fomass_max
+                         print*,'masse(ig,l)*fomass_max/ptimestep',masse(ig,l)*thermals_fomass_max/ptimestep
                          print*,'ptimestep',ptimestep
                          print*,'lmax(ig)',lmax(ig)