Changeset 6163


Ignore:
Timestamp:
Apr 10, 2026, 10:24:25 AM (2 hours ago)
Author:
idelkadi
Message:

Rewriting and simplifying the interface between the LMDZ code and the radiative transfer codes (continued)

Location:
LMDZ6/trunk/libf
Files:
2 edited
1 moved

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/lmdz_call_radiatifcodes_m.F90

    r6162 r6163  
    33!
    44
    5 MODULE radlwsw_m
     5MODULE lmdz_call_radiatifcodes_m
    66
    77  IMPLICIT NONE
     
    99contains
    1010
    11   SUBROUTINE radlwsw( &
     11  SUBROUTINE lmdz_call_radiatifcodes( &
    1212                      debut, dist, rmu0, fract, &
    1313                      paprs, pplay, tsol, SFRWL,alb_dir, alb_dif, &
     
    3737    ! A. Idelkadi, F. Hourdin, L. Fairhead, mars 2026 :                                                 
    3838    ! Recriture de l interface entre LMDZ et les codes radiatifs oldrad/rrtm/ecrad =>
    39     ! radlwsw_m.F90     
     39    ! lmdz_call_radiatifcodes_m.F90     
    4040    ! ------------------------------------------------------------------------------------------
    4141    !                   -> - Declarations => in / out / locales
     
    260260   REAL(KIND=8) ZFLDN0(KDLON,KFLEV+1)
    261261   CHARACTER (LEN=80) :: abort_message
    262    CHARACTER (LEN=80) :: modname='radlwsw_m'
     262   CHARACTER (LEN=80) :: modname='lmdz_call_radiatifcodes_m'
    263263
    264264   PRINT*,'size(wo, 1), klon, size(wo, 2), klev', &
    265265           size(wo, 1), klon, size(wo, 2), klev
    266266
    267    CALL assert(size(wo, 1) == klon, size(wo, 2) == klev, "radlwsw wo")
     267   CALL assert(size(wo, 1) == klon, size(wo, 2) == klev, "lmdz_call_radiatifcodes wo")
    268268
    269269   !
     
    360360       !======================================================================
    361361
    362   END SUBROUTINE radlwsw
    363 END MODULE radlwsw_m
     362  END SUBROUTINE lmdz_call_radiatifcodes
     363END MODULE lmdz_call_radiatifcodes_m
  • LMDZ6/trunk/libf/phylmd/physiq_mod.F90

    r6159 r6163  
    119119    USE readaerosol_mod, ONLY : init_aero_fromfile
    120120    USE readaerosolstrato_m, ONLY : init_readaerosolstrato
    121     USE radlwsw_m, ONLY: radlwsw
     121    USE lmdz_call_radiatifcodes_m, ONLY: lmdz_call_radiatifcodes
    122122    USE regr_horiz_time_climoz_m, ONLY: regr_horiz_time_climoz
    123123    USE regr_pr_time_av_m, ONLY: regr_pr_time_av
     
    11061106    ! retros liees aux CRF
    11071107    REAL, dimension(klon, klev) :: cldtaurad   ! epaisseur optique
    1108     ! pour radlwsw pour
     1108    ! pour lmdz_call_radiatifcodes pour
    11091109    ! tester "CRF off"
    11101110    REAL, dimension(klon, klev) :: cldtaupirad ! epaisseur optique
    1111     ! pour radlwsw pour
     1111    ! pour lmdz_call_radiatifcodes pour
    11121112    ! tester "CRF off"
    11131113    REAL, dimension(klon, klev) :: cldemirad   ! emissivite pour
    1114     ! radlwsw pour tester
     1114    ! lmdz_call_radiatifcodes pour tester
    11151115    ! "CRF off"
    11161116    REAL, dimension(klon, klev) :: cldfrarad   ! fraction nuageuse
     
    25702570    ! Calculs de l'orbite.
    25712571    ! Necessaires pour le rayonnement et la surface (calcul de l'albedo).
    2572     ! doit donc etre plac\'e avant radlwsw et pbl_surface
     2572    ! doit donc etre plac\'e avant lmdz_call_radiatifcodes et pbl_surface
    25732573
    25742574    ! !!   jyg 17 Sep 2010 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     
    44284428
    44294429          IF (prt_level .GE.10) THEN
    4430              print *,' ->radlwsw, number 1 '
     4430             print *,' ->lmdz_call_radiatifcodes, number 1 '
    44314431          ENDIF
    44324432          !
    4433           ! AI namelist utilise pour l appel principal de radlwsw (ecrad)
     4433          ! AI namelist utilise pour l appel principal de lmdz_call_radiatifcodes (ecrad)
    44344434          namelist_ecrad_file='namelist_ecrad'
    44354435          !
    4436           CALL radlwsw &
     4436          CALL lmdz_call_radiatifcodes &
    44374437               (debut, dist, rmu0, fract,  &
    44384438                                !albedo SB >>>
     
    45074507             !
    45084508             IF (prt_level .GE.10) THEN
    4509                 print *,' ->radlwsw, number 2 '
     4509                print *,' ->lmdz_call_radiatifcodes, number 2 '
    45104510             ENDIF
    45114511             !
    45124512             namelist_ecrad_file='namelist_ecrad'
    45134513             !
    4514              CALL radlwsw &
     4514             CALL lmdz_call_radiatifcodes &
    45154515                  (debut, dist, rmu0, fract,  &
    45164516                                !albedo SB >>>
     
    45624562             !                print*,'ok_2xcall_ecrad = ',ok_2xcall_ecrad
    45634563             namelist_ecrad_file='namelist_ecrad_s2'
    4564              CALL radlwsw &
     4564             CALL lmdz_call_radiatifcodes &
    45654565                  (debut, dist, rmu0, fract,  &
    45664566                  paprs, pplay,zxtsol,SFRWL,albsol_dir, albsol_dif, &
     
    46344634
    46354635    !
    4636     ! Calculer radsol a l'exterieur de radlwsw
     4636    ! Calculer radsol a l'exterieur de lmdz_call_radiatifcodes
    46374637    ! pour prendre en compte le cycle diurne
    46384638    ! recode par Olivier Boucher en sept 2015
  • LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

    r6159 r6163  
    9292    USE readaerosol_mod, ONLY : init_aero_fromfile
    9393    USE readaerosolstrato_m, ONLY : init_readaerosolstrato
    94     USE radlwsw_m, ONLY: radlwsw
     94    USE lmdz_call_radiatifcodes_m, ONLY: lmdz_call_radiatifcodes
    9595    USE reevap_mod, ONLY: reevap
    9696    USE regr_horiz_time_climoz_m, ONLY: regr_horiz_time_climoz
     
    11801180    ! retros liees aux CRF
    11811181    REAL, DIMENSION(klon, klev) :: cldtaurad   ! epaisseur optique
    1182     ! pour radlwsw pour
     1182    ! pour lmdz_call_radiatifcodes pour
    11831183    ! tester "CRF off"
    11841184    REAL, DIMENSION(klon, klev) :: cldtaupirad ! epaisseur optique
    1185     ! pour radlwsw pour
     1185    ! pour lmdz_call_radiatifcodes pour
    11861186    ! tester "CRF off"
    11871187    REAL, DIMENSION(klon, klev) :: cldemirad   ! emissivite pour
    1188     ! radlwsw pour tester
     1188    ! lmdz_call_radiatifcodes pour tester
    11891189    ! "CRF off"
    11901190    REAL, DIMENSION(klon, klev) :: cldfrarad   ! fraction nuageuse
     
    29672967    ! Calculs de l'orbite.
    29682968    ! Necessaires pour le rayonnement et la surface (calcul de l'albedo).
    2969     ! doit donc etre place avant radlwsw et pbl_surface
     2969    ! doit donc etre place avant lmdz_call_radiatifcodes et pbl_surface
    29702970
    29712971    ! !!   jyg 17 Sep 2010 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     
    59775977
    59785978          IF (prt_level .GE.10) THEN
    5979              print *,' ->radlwsw, number 1 '
     5979             print *,' ->lmdz_call_radiatifcodes, number 1 '
    59805980          ENDIF
    59815981          !
    5982           ! AI namelist utilise pour l appel principal de radlwsw (ecrad)
     5982          ! AI namelist utilise pour l appel principal de lmdz_call_radiatifcodes (ecrad)
    59835983          namelist_ecrad_file='namelist_ecrad'
    59845984          !
    5985           CALL radlwsw &
     5985          CALL lmdz_call_radiatifcodes &
    59865986               (debut, dist, rmu0, fract,  &
    59875987                                !albedo SB >>>
     
    60566056             !
    60576057             IF (prt_level .GE.10) THEN
    6058                 print *,' ->radlwsw, number 2 '
     6058                print *,' ->lmdz_call_radiatifcodes, number 2 '
    60596059             ENDIF
    60606060             !
    60616061             namelist_ecrad_file='namelist_ecrad'
    60626062             !
    6063              CALL radlwsw &
     6063             CALL lmdz_call_radiatifcodes &
    60646064                  (debut, dist, rmu0, fract,  &
    60656065                                !albedo SB >>>
     
    61116111             !                print*,'ok_2xcall_ecrad = ',ok_2xcall_ecrad
    61126112             namelist_ecrad_file='namelist_ecrad_s2'
    6113              CALL radlwsw &
     6113             CALL lmdz_call_radiatifcodes &
    61146114                  (debut, dist, rmu0, fract,  &
    61156115                  paprs, pplay,zxtsol,SFRWL,albsol_dir, albsol_dif, &
     
    61836183
    61846184    !
    6185     ! Calculer radsol a l'exterieur de radlwsw
     6185    ! Calculer radsol a l'exterieur de lmdz_call_radiatifcodes
    61866186    ! pour prendre en compte le cycle diurne
    61876187    ! recode par Olivier Boucher en sept 2015
Note: See TracChangeset for help on using the changeset viewer.