Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/lmdz

Timestamp:

Aug 2, 2024, 9:58:25 PM (6 months ago)

Author:

abarral

Message:

Put dimensions.h and paramet.h into modules

Location:

LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/lmdz

Files:

• aeropt_5wv_ecrad.F90 (modified) (1 diff)
• calcul_cloud_overlap_decorr_len.F90 (modified) (1 diff)
• radiation_scheme_mod.F90 (modified) (3 diffs)
• radiation_setup.F90 (modified) (1 diff)
• readaerosol_optic_ecrad.F90 (modified) (9 diffs)
• readaerosolstrato_ecrad.F90 (modified) (4 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/lmdz/aeropt_5wv_ecrad.F90

@@ -1 +1 @@
 !aeropt_5wv_ecrad.F90 2022-09-20 A. Idelkadi et O. Boucher 
-!
+
      SUBROUTINE AEROPT_5WV_ECRAD(istartcol,iendcol,istartlev,iendlev, &
                                  config,thermodynamics,aerosol)

LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/lmdz/calcul_cloud_overlap_decorr_len.F90

@@ -18 +18 @@
 !      Ld=mid_decorrelation_length si pres entre [440,680]
 !      Ld=high_decorrelation_length si pres < 440hPa
-!
+
 ! TO DO 
 ! Que faire de PDECORR_LEN_WATER_M ?

LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/lmdz/radiation_scheme_mod.F90

@@ -687 +687 @@
 
 ! RADIATION_SCHEME - Interface to modular radiation scheme
-!
+
 ! (C) Copyright 2015- ECMWF.
-!
+
 ! This software is licensed under the terms of the Apache Licence Version 2.0
 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
-!
+
 ! In applying this licence, ECMWF does not waive the privileges and immunities
 ! granted to it by virtue of its status as an intergovernmental organisation
 ! nor does it submit to any jurisdiction.
-!
+
 ! PURPOSE
 ! -------
@@ -704 +704 @@
 !   it to the radiation scheme.  It returns net fluxes and surface
 !   flux components needed by the rest of the model. 
-!
+
 !   Lower case is used for variables and types taken from the
 !   radiation library
-!
+
 ! INTERFACE
 ! ---------
@@ -713 +713 @@
 !    SETUP_RADIATION_SCHEME routine (in the RADIATION_SETUP module)
 !    should have been run first.
-!
+
 ! AUTHOR
 ! ------
 !   Robin Hogan, ECMWF
 !   Original: 2015-09-16
-!
+
 ! MODIFICATIONS
 ! -------------
-!
+
 ! TO DO
 ! -----
-!
+
 !-----------------------------------------------------------------------
 

LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/lmdz/radiation_setup.F90

@@ -2 +2 @@
 
 ! RADIATION_SETUP - Setting up modular radiation scheme
-!
+
 ! AUTHOR
 ! ------
 !   Robin Hogan, ECMWF
 !   Original: 2015-09-16
-!
+
 ! MODIFICATIONS
 ! -------------
 !   Abderrahmane Idelkadi LMD, juillet 2023
-!
+
 !-----------------------------------------------------------------------
 

LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/lmdz/readaerosol_optic_ecrad.F90

@@ -1 +1 @@
 ! $Id: readaerosol_optic_ecrad.F90
-!
+
 SUBROUTINE readaerosol_optic_ecrad(debut, aerosol_couple, ok_alw, ok_volcan, &
      flag_aerosol, flag_bc_internal_mixture, itap, rjourvrai, &
@@ -102 +102 @@
   !    
   !****************************************************************************************
-  !
-  !
+
+
   IF (aerosol_couple) THEN   !--we get aerosols from tr_seri array from INCA
-     !
+
      !--copy fields from INCA tr_seri 
      !--convert to ug m-3 unit for consistency with offline fields
-     !
+
      itr = 0
      DO iq = 1,nqtot
@@ -145 +145 @@
      nitrcoarse(:,:)   =   tr_seri(:,:,id_CSNO3M)                        *zrho(:,:)*1.e9  ! CSNO3M
      nitrinscoarse(:,:)=   tr_seri(:,:,id_CINO3M)                        *zrho(:,:)*1.e9  ! CINO3M
-     !
+
      bcsol_pi(:,:)        =   0.0 ! ASBCM pre-ind
      pomsol_pi(:,:)       =   0.0 ! ASPOMM pre-ind
@@ -159 +159 @@
      nitrcoarse_pi(:,:)   =   0.0 ! CSNO3M pre-ind
      nitrinscoarse_pi(:,:)=   0.0 ! CINO3M
-     !
+
   ELSE !--not aerosol_couple
-     !
+
      ! Read and interpolate sulfate
      IF ( flag_aerosol .EQ. 1 .OR. flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN 
@@ -222 +222 @@
         loaddust=0.
      ENDIF
-     !
+
      ! Read and interpolate asno3m, csno3m, cino3m
      IF (flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN 
@@ -240 +240 @@
         loadno3(:)=0.0
      ENDIF
-     !
+
      ! CSSO4M is set to 0 as not reliable
      sulfcoarse(:,:)      =   0.0 ! CSSO4M (=SO4) + CSMSAM (=MSA) 
@@ -247 +247 @@
   ENDIF !--not aerosol_couple
 
-  !
+
   ! Store all aerosols mixing ratios in one variable for radiation scheme (unit kg/kg for ECRAD)
   ! present-day values
@@ -281 +281 @@
   m_allaer_pi(:,:,id_STRAT_phy)  = 0.0
 
-  !
+
   ! Calculate the total mass of all soluble aersosols (in unit ug /m3)
   ! to be revisited for AR6
@@ -289 +289 @@
   !****************************************************************************************
   ! 2) Calculate optical properties for the aerosols
-  !
+
   !****************************************************************************************
   DO k = 1, klev

LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/lmdz/readaerosolstrato_ecrad.F90

@@ -1 @@
-!
+
 ! $Id: readaerosolstrato_ecrad.F90 tlurton $
-!
+
 SUBROUTINE readaerosolstrato_ecrad(config, debut, ok_volcan)
 
@@ -438 +438 @@
 
     IF (.NOT. ok_volcan) THEN
-!
+
 !--this is the default case 
 !--stratospheric aerosols are added to both index 2 and 1 for double radiation calls
@@ -464 +464 @@
     ENDWHERE
     ENDDO
-!
+
     ELSE
-!
+
 !--this is the VOLMIP case
 !--stratospheric aerosols are only added to index 2 in this case 
@@ -523 +523 @@
       ENDWHERE
     ENDDO
-!
+
     ELSE
-!
+
 ! mod ThL
 !--this is the VOLMIP case