Changeset 5274 for LMDZ6/trunk/libf/phylmd/dyn1d

Timestamp:

Oct 25, 2024, 3:41:23 PM (3 months ago)

Author:

abarral

Message:

Replace yomcst.h by existing module

Location:

LMDZ6/trunk/libf/phylmd/dyn1d

Files:

• 1DUTILS.h (modified) (7 diffs)
• 1Dconv.h (modified) (2 diffs)
• mod_1D_cases_read_std.f90 (modified) (1 diff)
• old_1DUTILS_read_interp.h (modified) (32 diffs)
• old_lmdz1d.f90 (modified) (1 diff)
• scm.f90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/dyn1d/1DUTILS.h

@@ -1352 +1352 @@
 !   sans WTG rajouter une advection horizontale
 !----------------------------------------------------------------------
-        implicit none
-INCLUDE "YOMCST.h"
+        USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
+          , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
+          , R_ecc, R_peri, R_incl                                      &
+          , RA, RG, R1SA                                         &
+          , RSIGMA                                                     &
+          , R, RMD, RMV, RD, RV, RCPD                    &
+          , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12        &
+          , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w                    &
+          , RCW, RCS                                                 &
+          , RLVTT, RLSTT, RLMLT, RTT, RATM                           &
+          , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS            &
+          , RALPD, RBETD, RGAMD
+implicit none
+
 !        argument
         integer llm
@@ -1426 +1438 @@
 !   sans WTG rajouter une advection horizontale
 !----------------------------------------------------------------------
-        implicit none
-INCLUDE "YOMCST.h"
+        USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
+          , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
+          , R_ecc, R_peri, R_incl                                      &
+          , RA, RG, R1SA                                         &
+          , RSIGMA                                                     &
+          , R, RMD, RMV, RD, RV, RCPD                    &
+          , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12        &
+          , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w                    &
+          , RCW, RCS                                                 &
+          , RLVTT, RLSTT, RLMLT, RTT, RATM                           &
+          , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS            &
+          , RALPD, RBETD, RGAMD
+implicit none
+
 !        argument
         integer llm,nqtot
@@ -1507 +1531 @@
   USE dimphy
 
-  implicit none
+  USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
+          , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
+          , R_ecc, R_peri, R_incl                                      &
+          , RA, RG, R1SA                                         &
+          , RSIGMA                                                     &
+          , R, RMD, RMV, RD, RV, RCPD                    &
+          , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12        &
+          , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w                    &
+          , RCW, RCS                                                 &
+          , RLVTT, RLSTT, RLMLT, RTT, RATM                           &
+          , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS            &
+          , RALPD, RBETD, RGAMD
+implicit none
 
 ! ========================================================
@@ -1526 +1562 @@
 ! Declaration des constantes et des fonctions thermodynamiques
 !
-include "YOMCST.h"
+
 include "YOETHF.h"
 !
@@ -1590 +1626 @@
   USE dimphy
 
-  implicit none
+  USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
+          , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
+          , R_ecc, R_peri, R_incl                                      &
+          , RA, RG, R1SA                                         &
+          , RSIGMA                                                     &
+          , R, RMD, RMV, RD, RV, RCPD                    &
+          , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12        &
+          , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w                    &
+          , RCW, RCS                                                 &
+          , RLVTT, RLSTT, RLMLT, RTT, RATM                           &
+          , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS            &
+          , RALPD, RBETD, RGAMD
+implicit none
 
 ! ========================================================
@@ -1620 +1668 @@
 ! Declaration des constantes et des fonctions thermodynamiques
 !
-include "YOMCST.h"
+
 include "YOETHF.h"
 !
@@ -1747 +1795 @@
 
        USE dimensions_mod, ONLY: iim, jjm, llm, ndm
+USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
+          , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
+          , R_ecc, R_peri, R_incl                                      &
+          , RA, RG, R1SA                                         &
+          , RSIGMA                                                     &
+          , R, RMD, RMV, RD, RV, RCPD                    &
+          , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12        &
+          , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w                    &
+          , RCW, RCS                                                 &
+          , RLVTT, RLSTT, RLMLT, RTT, RATM                           &
+          , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS            &
+          , RALPD, RBETD, RGAMD
 implicit none
 
-INCLUDE "YOMCST.h"
+
 
 

LMDZ6/trunk/libf/phylmd/dyn1d/1Dconv.h

@@ -6 +6 @@
      &       Ts,imp_fcg,ts_fcg,Tp_fcg,Turb_fcg)                                 
 !
-        implicit none
- 
+        USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
+          , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
+          , R_ecc, R_peri, R_incl                                      &
+          , RA, RG, R1SA                                         &
+          , RSIGMA                                                     &
+          , R, RMD, RMV, RD, RV, RCPD                    &
+          , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12        &
+          , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w                    &
+          , RCW, RCS                                                 &
+          , RLVTT, RLSTT, RLMLT, RTT, RATM                           &
+          , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS            &
+          , RALPD, RBETD, RGAMD
+implicit none
+
 !ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 ! cette routine permet d obtenir u_convg,v_convg,ht,hq et ainsi de
@@ -13 +25 @@
 !cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 
-      INCLUDE "YOMCST.h"
+
 
       INTEGER klev

LMDZ6/trunk/libf/phylmd/dyn1d/mod_1D_cases_read_std.f90

@@ -932 +932 @@
 
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
+USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
+          , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
+          , R_ecc, R_peri, R_incl                                      &
+          , RA, RG, R1SA                                         &
+          , RSIGMA                                                     &
+          , R, RMD, RMV, RD, RV, RCPD                    &
+          , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12        &
+          , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w                    &
+          , RCW, RCS                                                 &
+          , RLVTT, RLSTT, RLMLT, RTT, RATM                           &
+          , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS            &
+          , RALPD, RBETD, RGAMD
 implicit none
 
-    INCLUDE "YOMCST.h"
+
 
 

LMDZ6/trunk/libf/phylmd/dyn1d/old_1DUTILS_read_interp.h

@@ -1463 +1463 @@
      &             ,ufa_astex,vfa_astex,div_prof,ts_prof,ug_prof,vg_prof   &
      &             ,ufa_prof,vfa_prof)
-        implicit none
+        USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
+          , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
+          , R_ecc, R_peri, R_incl                                      &
+          , RA, RG, R1SA                                         &
+          , RSIGMA                                                     &
+          , R, RMD, RMV, RD, RV, RCPD                    &
+          , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12        &
+          , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w                    &
+          , RCW, RCS                                                 &
+          , RLVTT, RLSTT, RLMLT, RTT, RATM                           &
+          , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS            &
+          , RALPD, RBETD, RGAMD
+      implicit none
 
 !---------------------------------------------------------------------------------------
@@ -1543 +1555 @@
      &'day,annee_ref,day_ini_astex,timeit,it_astex1,it_astex2,SST:',       &
      &day,annee_ref,day_ini_astex,timeit/86400.,it_astex1,                 &
-     &it_astex2,div_prof,ts_prof,ug_prof,vg_prof,ufa_prof,vfa_prof 
+     &it_astex2,div_prof,ts_prof,ug_prof,vg_prof,ufa_prof,vfa_prof
 
         return
@@ -1815 +1827 @@
 ! Time interpolation of a 2D field to the timestep corresponding to day
 !
-! day: current julian day 
+! day: current julian day
 ! day1: first day of the simulation
 ! nt_gabls4: total nb of data in the forcing (e.g. 37 for gabls4)
@@ -2148 +2160 @@
         close(ilesfile)
 
-! Vertical coordinates half levels for eta-coordinates (plev = alpha + beta * psurf) 
+! Vertical coordinates half levels for eta-coordinates (plev = alpha + beta * psurf)
        IF(nlev_max.EQ.19) THEN
        open (ifile,file='proh.inp.19',status='old',iostat=ierr)
@@ -2250 +2262 @@
            stop 'vg'
          endif
-      
+
       ierr=nf90_inq_varid(nid,"wls",var3didin(9))
          if(ierr/=nf90_noerr) then
@@ -2282 +2294 @@
 !dimensions lecture
 !      call catchaxis(nid,ntime,nlevel,time,z,ierr)
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(1),zz)
          if(ierr/=nf90_noerr) then
@@ -2303 +2315 @@
          endif
 !          write(*,*)'lecture qt ok',qt
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(4),u)
          if(ierr/=nf90_noerr) then
@@ -2374 +2386 @@
 !          write(*,*)'lecture thl_rad ok',thl_rad
 
-         return 
+         return
          end subroutine read_fire
 !=====================================================================
@@ -2387 +2399 @@
             nf90_inq_dimid,nf90_inquire_dimension
 
-      implicit none
-      INCLUDE "YOMCST.h"
+      USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
+          , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
+          , R_ecc, R_peri, R_incl                                      &
+          , RA, RG, R1SA                                         &
+          , RSIGMA                                                     &
+          , R, RMD, RMV, RD, RV, RCPD                    &
+          , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12        &
+          , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w                    &
+          , RCW, RCS                                                 &
+          , RLVTT, RLSTT, RLMLT, RTT, RATM                           &
+          , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS            &
+          , RALPD, RBETD, RGAMD
+implicit none
+
 
       integer ntime,nlevel
@@ -2471 +2495 @@
            stop 'lhf'
          endif
-      
+
       ierr=nf90_inq_varid(nid,"lwup",var3didin(4))
          if(ierr/=nf90_noerr) then
@@ -2551 +2575 @@
 !dimensions lecture
 !      call catchaxis(nid,ntime,nlevel,time,z,ierr)
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(1),zz)
          if(ierr/=nf90_noerr) then
@@ -2558 +2582 @@
          endif
 !          write(*,*)'lecture zz ok',zz
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(11),pres)
          if(ierr/=nf90_noerr) then
@@ -2582 +2606 @@
          endif
 !          write(*,*)'lecture qv ok',qv
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(14),u)
          if(ierr/=nf90_noerr) then
@@ -2709 +2733 @@
 !          write(*,*)'lecture omega ok',omega
 
-         return 
+         return
          end subroutine read_dice
 !=====================================================================
@@ -2717 +2741 @@
 !program reading initial profils and forcings of the Gabls4 case study
 
+      USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
+          , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
+          , R_ecc, R_peri, R_incl                                      &
+          , RA, RG, R1SA                                         &
+          , RSIGMA                                                     &
+          , R, RMD, RMV, RD, RV, RCPD                    &
+          , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12        &
+          , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w                    &
+          , RCW, RCS                                                 &
+          , RLVTT, RLSTT, RLMLT, RTT, RATM                           &
+          , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS            &
+          , RALPD, RBETD, RGAMD
       USE netcdf, ONLY: nf90_open,nf90_nowrite,nf90_noerr,nf90_strerror,nf90_inq_varid,nf90_get_var,&
             nf90_inq_dimid,nf90_inquire_dimension
@@ -2727 +2763 @@
       real*8 time(ntime)
 
-!  ATTENTION: visiblement quand on lit gabls4_driver.nc on recupere les donnees 
+!  ATTENTION: visiblement quand on lit gabls4_driver.nc on recupere les donnees
 ! dans un ordre inverse par rapport a la convention LMDZ
 ! ==> il faut tout inverser  (MPL 20141024)
@@ -2846 +2882 @@
 !dimensions lecture
 !      call catchaxis(nid,ntime,nlevel,time,z,ierr)
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(1),zz_i)
          if(ierr/=nf90_noerr) then
@@ -2852 +2888 @@
             stop "getvarup"
          endif
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(2),depth_sn)
          if(ierr/=nf90_noerr) then
@@ -2858 +2894 @@
             stop "getvarup"
          endif
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(3),ug_i)
          if(ierr/=nf90_noerr) then
@@ -2864 +2900 @@
             stop "getvarup"
          endif
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(4),vg_i)
          if(ierr/=nf90_noerr) then
@@ -2870 +2906 @@
             stop "getvarup"
          endif
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(5),pf_i)
          if(ierr/=nf90_noerr) then
@@ -2894 +2930 @@
             stop "getvarup"
          endif
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(9),u_i)
          if(ierr/=nf90_noerr) then
@@ -2900 +2936 @@
             stop "getvarup"
          endif
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(10),v_i)
          if(ierr/=nf90_noerr) then
@@ -2906 +2942 @@
             stop "getvarup"
          endif
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(11),hadvt_i)
          if(ierr/=nf90_noerr) then
@@ -2912 +2948 @@
             stop "getvarup"
          endif
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(12),hadvq_i)
          if(ierr/=nf90_noerr) then
@@ -2918 +2954 @@
             stop "getvarup"
          endif
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(14),tsnow)
          if(ierr/=nf90_noerr) then
@@ -2924 +2960 @@
             stop "getvarup"
          endif
- 
+
          ierr = NF90_GET_VAR(nid,var3didin(15),snow_dens)
          if(ierr/=nf90_noerr) then
@@ -2952 +2988 @@
            hadvq(k,:)=hadvq_i(nlevel+1-k,:)
          enddo
-         return 
+         return
  end subroutine read_gabls4
 !=====================================================================
 
-!     Reads CIRC input files      
+!     Reads CIRC input files
 
       SUBROUTINE read_circ(nlev_circ,cf,lwp,iwp,reliq,reice,t,z,p,pm,h2o,o3,sza)
-      
+
       parameter (ncm_1=49180)
-      INCLUDE "YOMCST.h"
 
       real albsfc(ncm_1), albsfc_w(ncm_1)
@@ -2985 +3020 @@
       open (15, file='cloud_input_case.txt', status='old')
       open (16, file='sfcalbedo_input_case.txt', status='old')
-      
+
 !     Read scalar information
       do iskip=1,5
@@ -3018 +3053 @@
          pm(il)=pm(il)*100.
       enddo
-303   format (i8, 2f9.2, 10(2x,e13.7))      
+303   format (i8, 2f9.2, 10(2x,e13.7))
       close(13)
-      
+
 !     Read aerosol layer information
       do iskip=1,3
@@ -3033 +3068 @@
 304   format (i8, f9.5, 2f8.3)
       close(14)
-      
+
 !     Read cloud information
       do iskip=1,3
@@ -3057 +3092 @@
 306   format(f10.1, 2f12.5, f14.8)
       close(16)
- 
-      return 
+
+      return
       end subroutine read_circ
 !=====================================================================
-!     Reads RTMIP input files      
+!     Reads RTMIP input files
 
       SUBROUTINE read_rtmip(nlev_rtmip,play,plev,t,h2o,o3)
-      
-      INCLUDE "YOMCST.h"
+
 
       real t(nlev_rtmip), pt(nlev_rtmip),pb(nlev_rtmip),h2o(nlev_rtmip), o3(nlev_rtmip)

LMDZ6/trunk/libf/phylmd/dyn1d/old_lmdz1d.f90

@@ -50 +50 @@
    USE dimensions_mod, ONLY: iim, jjm, llm, ndm
    USE dimsoil_mod_h, ONLY: nsoilmx
+USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
+          , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
+          , R_ecc, R_peri, R_incl                                      &
+          , RA, RG, R1SA                                         &
+          , RSIGMA                                                     &
+          , R, RMD, RMV, RD, RV, RCPD                    &
+          , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12        &
+          , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w                    &
+          , RCW, RCS                                                 &
+          , RLVTT, RLSTT, RLMLT, RTT, RATM                           &
+          , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS            &
+          , RALPD, RBETD, RGAMD
 implicit none
 
-      INCLUDE "YOMCST.h"
+
 !!      INCLUDE "control.h"
       INCLUDE "clesphys.h"

LMDZ6/trunk/libf/phylmd/dyn1d/scm.f90

@@ -46 +46 @@
    USE dimensions_mod, ONLY: iim, jjm, llm, ndm
    USE dimsoil_mod_h, ONLY: nsoilmx
+USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
+          , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
+          , R_ecc, R_peri, R_incl                                      &
+          , RA, RG, R1SA                                         &
+          , RSIGMA                                                     &
+          , R, RMD, RMV, RD, RV, RCPD                    &
+          , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12        &
+          , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w                    &
+          , RCW, RCS                                                 &
+          , RLVTT, RLSTT, RLMLT, RTT, RATM                           &
+          , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS            &
+          , RALPD, RBETD, RGAMD
 implicit none
 
-      INCLUDE "YOMCST.h"
+
 !!      INCLUDE "control.h"
       INCLUDE "clesphys.h"