Skip to content

Commit

Permalink
Merge pull request #59 from QuentinRodier/check_mesonh_dev_0.7.0_to_m…
Browse files Browse the repository at this point in the history 
…aster

MésoNH contribution from 5.7.0 --> 5.7.1
  • Loading branch information
@QuentinRodier
QuentinRodier authored Jan 22, 2025
2 parents d575740 + 5031bb9 commit 669c0ef
Show file tree
Hide file tree
Showing 31 changed files with 1,523 additions and 773 deletions.
12 changes: 11 additions & 1 deletion src/arome/ext/aro_adjust_lima.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -86,6 +86,9 @@ SUBROUTINE ARO_ADJUST_LIMA(PHYEX, &
!
USE MODD_PARAM_LIMA
USE MODD_NSV
USE MODD_CH_AEROSOL, ONLY: NSP, NSOA, NCARB, LORILAM
USE MODD_SALT, ONLY: LSALT
USE MODD_DUST, ONLY: LDUST
!
USE MODI_LIMA_ADJUST_SPLIT
USE MODE_SET_CONC_LIMA
Expand Down Expand Up @@ -154,6 +157,9 @@ SUBROUTINE ARO_ADJUST_LIMA(PHYEX, &
REAL, DIMENSION(KLON,1,KLEV), INTENT(INOUT) :: PICEFR! Cloud fraction
REAL, DIMENSION(KLON,1,KLEV), INTENT(INOUT) :: PPRCFR! Cloud fraction
!
CHARACTER(LEN=4) :: HACTCCN ! kind of CCN activation
REAL, DIMENSION(KLON,1,KLEV,10) :: ZSOLORG ![%] solubility fraction of soa
REAL, DIMENSION(KLON,1,KLEV, NSP+NCARB+NSOA) :: ZMI
!
TYPE(TYP_DDH), INTENT(INOUT), TARGET :: YDDDH
TYPE(TLDDH), INTENT(IN), TARGET :: YDLDDH
Expand Down Expand Up @@ -189,7 +195,9 @@ SUBROUTINE ARO_ADJUST_LIMA(PHYEX, &
IF (LHOOK) CALL DR_HOOK('ARO_ADJUST_LIMA',0,ZHOOK_HANDLE)

CALL FILL_DIMPHYEX(YLDIMPHYEX, KLON, 1, KLEV, 0, KFDIA)

HACTCCN=' '
ZMI=0.
ZSOLORG=0.
!
!* 2. TRANSFORMATION INTO PHYSICAL TENDENCIES
! ---------------------------------------
Expand Down Expand Up @@ -334,10 +342,12 @@ SUBROUTINE ARO_ADJUST_LIMA(PHYEX, &

CALL LIMA_ADJUST_SPLIT(D=YLDIMPHYEX, CST=PHYEX%CST, BUCONF=TBUCONF, TBUDGETS=YLBUDGET, KBUDGETS=SIZE(YLBUDGET), &
KRR=KRR, KMI=1, HCONDENS=PHYEX%NEBN%CCONDENS, HLAMBDA3=PHYEX%NEBN%CLAMBDA3, &
KCARB=NCARB, KSOA=NSOA, KSP=NSP, ODUST=LDUST, OSALT=LSALT, OORILAM=LORILAM, &
OSUBG_COND=OSUBG_COND, OSIGMAS=OSIGMAS, PTSTEP=2*PTSTEP, PSIGQSAT=PSIGQSAT, &
PRHODREF=PRHODREF, PRHODJ=PRHODJ, PEXNREF=PEXNREF, PSIGS=PSIGS, LMFCONV=.TRUE., PMFCONV=PMFCONV, &
PPABST=PPABSM, PPABSTT=PPABSM, PZZ=ZZZ, ODTHRAD=.TRUE., PDTHRAD=PDTHRAD, PW_NU=PW_NU, &
PRT=PRT, PRS=PRS, PSVT=PSVT, PSVS=PSVS, &
HACTCCN=HACTCCN, PAERO=PSVT, PSOLORG=ZSOLORG, PMI=ZMI, &
PTHS=PTHS, OCOMPUTE_SRC=.TRUE., PSRCS=PSRCS, PCLDFR=PCLDFR, PICEFR=PICEFR, &
PRC_MF=PRC_MF, PRI_MF=PRI_MF, PCF_MF=PCF_MF )
!
Expand Down
16 changes: 14 additions & 2 deletions src/arome/ext/aro_lima.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -47,6 +47,9 @@ SUBROUTINE ARO_LIMA(PHYEX,KKA,KKU,KKL,KLON,KLEV,KFDIA,KRR, KSV, &
!
USE MODD_PARAM_LIMA
USE MODD_NSV
USE MODD_CH_AEROSOL, ONLY: NSP, NCARB, NSOA, LORILAM
USE MODD_DUST, ONLY: LDUST
USE MODD_SALT, ONLY: LSALT
!
USE MODD_BUDGET
USE MODE_BUDGET_PHY, ONLY: BUDGET_DDH
Expand Down Expand Up @@ -108,6 +111,10 @@ SUBROUTINE ARO_LIMA(PHYEX,KKA,KKU,KKL,KLON,KLEV,KFDIA,KRR, KSV, &
REAL, DIMENSION(KLON,1,KLEV), INTENT(INOUT) :: PICEFR ! ice cloud fraction
REAL, DIMENSION(KLON,1,KLEV), INTENT(INOUT) :: PPRCFR ! precipitation fraction
!
CHARACTER(LEN=4) :: HACTCCN ! kind of CCN activation
REAL, DIMENSION(KLON,1,KLEV,10) :: ZSOLORG ![%] solubility fraction of soa
REAL, DIMENSION(KLON,1,KLEV, NSP+NCARB+NSOA) :: ZMI
!
TYPE(TYP_DDH), INTENT(INOUT), TARGET :: YDDDH
TYPE(TLDDH), INTENT(IN), TARGET :: YDLDDH
TYPE(TMDDH), INTENT(IN), TARGET :: YDMDDH
Expand Down Expand Up @@ -156,6 +163,9 @@ SUBROUTINE ARO_LIMA(PHYEX,KKA,KKU,KKL,KLON,KLEV,KFDIA,KRR, KSV, &
ZDUM3DG=0.
ZDUM3DH=0.
PINPRH=0.
HACTCCN=' '
ZSOLORG=0.
ZMI=0.
IF (PHYEX%MISC%CELEC /='NONE') THEN
CALL ABOR1('ARO_LIMA : CELEC ELECTRICITY SCHEME NOT YET CORRECLY PLUGGED HERE')
! The following value of ZTHVREFZIKB must be removed from the electricity scheme or computed correctly here
Expand Down Expand Up @@ -294,12 +304,14 @@ SUBROUTINE ARO_LIMA(PHYEX,KKA,KKU,KKL,KLON,KLEV,KFDIA,KRR, KSV, &
!
CALL LIMA (D=YLDIMPHYEX, CST=PHYEX%CST, ICED=PHYEX%RAIN_ICE_DESCRN, ICEP=PHYEX%RAIN_ICE_PARAMN, &
ELECD=PHYEX%ELEC_DESCR, ELECP=PHYEX%ELEC_PARAM, &
BUCONF=TBUCONF, TBUDGETS=YLBUDGET, KBUDGETS=SIZE(YLBUDGET), KRR=KRR, &
BUCONF=TBUCONF, TBUDGETS=YLBUDGET,HACTCCN=HACTCCN, KBUDGETS=SIZE(YLBUDGET), KRR=KRR, &
PTSTEP=2*PTSTEP, OELEC=PHYEX%MISC%OELEC, &
PRHODREF=PRHODREF, PEXNREF=PEXNREF, PDZZ=PDZZ, PTHVREFZIKB=ZTHVREFZIKB, &
PRHODJ=PRHODJ, PPABST=PPABSM, &
NCCN=NMOD_CCN, NIFN=NMOD_IFN, NIMM=NMOD_IMM, &
NCCN=NMOD_CCN, NIFN=NMOD_IFN, NIMM=NMOD_IMM, KCARB=NCARB, KSOA=NSOA, KSP=NSP, &
ODUST=LDUST, OSALT=LSALT, OORILAM=LORILAM, &
ODTHRAD=.TRUE., PDTHRAD=PDTHRAD, PTHT=PTHT, PRT=PRT, PSVT=PSVT, PW_NU=PW_NU, &
PAERO=PSVT, PSOLORG=ZSOLORG, PMI=ZMI, &
PTHS=PTHS, PRS=PRS, PSVS=PSVS, &
PINPRC=ZINPRC, PINDEP=ZINDEP, PINPRR=PINPRR, PINPRI=ZINPRI, PINPRS=PINPRS, PINPRG=PINPRG, PINPRH=PINPRH, &
PEVAP3D=PEVAP, PCLDFR=PCLDFR, PICEFR=PICEFR, PPRCFR=PPRCFR, PFPR=PFPR )
Expand Down
82 changes: 82 additions & 0 deletions src/common/aux/ch_aer_activation.F90
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,82 @@
!ORILAM_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!ORILAM_LIC This is part of the ORILAM software governed by the CeCILL-C licence
!ORILAM_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!ORILAM_LIC for details.
!-----------------------------------------------------------------
!! ##############################
MODULE MODI_CH_AER_ACTIVATION
!! ##############################
!!
INTERFACE
!
SUBROUTINE CH_AER_ACTIVATION(PSVT,PTEMP, PWT, PTRAD, &
PRHODREF,PPABST, PNCN, PMCN,&
PSOLORG, PMI, PSMAX)

IMPLICIT NONE

REAL, DIMENSION(:,:), INTENT(INOUT) :: PSVT ! Aerosol concentration
REAL, DIMENSION(:), INTENT(IN) :: PTEMP ! Air temperature (K)
REAL, DIMENSION(:), INTENT(IN) :: PRHODREF ! Air density (kg/m3)
REAL, DIMENSION(:), INTENT(IN) :: PWT ! Activation vertical velocity (m/s)
REAL, DIMENSION(:), INTENT(IN) :: PTRAD ! Activation cooling radiative tendency (K/s)
REAL, DIMENSION(:), INTENT(IN) :: PPABST! Air pressure (Pa)
REAL, DIMENSION(:,:), INTENT(IN) :: PSOLORG ! ratio of SOA in acqueous phase
REAL, DIMENSION(:,:), INTENT(IN) :: PMI ! Molecular weight (g/mol)
REAL, DIMENSION(:), INTENT(OUT) :: PNCN ! Number of activated aerosol (#/m3)
REAL, DIMENSION(:), INTENT(OUT) :: PMCN ! Mass of activated aerosol (ug/m3)
REAL, DIMENSION(:), INTENT(INOUT) :: PSMAX ! Maximum supersaturation


END SUBROUTINE CH_AER_ACTIVATION
!!
END INTERFACE
!!
END MODULE MODI_CH_AER_ACTIVATION
!!
!! #######################################
SUBROUTINE CH_AER_ACTIVATION(PSVT,PTEMP, PWT, PTRAD, &
PRHODREF,PPABST, PNCN, PMCN,&
PSOLORG, PMI, PSMAX)
!! #######################################
!!
!! PURPOSE
!! -------
!! Input of Abdul-Razzak activation scheme
!! Here we compute the size distribution of aerosols together with the
!! dissociative ions and the soluble fraction need in the parameterization.
!! Dynamical variables are also passed by argument.
!! REFERENCE
!! ---------
!! none
!!
!! AUTHOR
!! ------
!! Pierre TULET (GMEI)
!!
!! MODIFICATIONS
!! -------------
!! Original
!!
!! Empty routines for PHYEX-offline (Q. Rodier 12/2024)
!!
!! IMPLICIT ARGUMENTS
!
!!
IMPLICIT NONE

!! Arguments variables

REAL, DIMENSION(:,:), INTENT(INOUT) :: PSVT ! Aerosol concentration
REAL, DIMENSION(:), INTENT(IN) :: PTEMP ! Air temperature (K)
REAL, DIMENSION(:), INTENT(IN) :: PRHODREF ! Air density (kg/m3)
REAL, DIMENSION(:), INTENT(IN) :: PWT ! Activation vertical velocity (m/s)
REAL, DIMENSION(:), INTENT(IN) :: PTRAD ! Activation cooling radiative tendency (K/s)
REAL, DIMENSION(:), INTENT(IN) :: PPABST! Air pressure (Pa)
REAL, DIMENSION(:,:), INTENT(IN) :: PSOLORG ! ratio of SOA in acqueous phase
REAL, DIMENSION(:,:), INTENT(IN) :: PMI ! Molecular weight (g/mol)
REAL, DIMENSION(:), INTENT(OUT) :: PNCN ! Number of activated aerosol (#/m3)
REAL, DIMENSION(:), INTENT(OUT) :: PMCN ! Mass of activated aerosol (ug/m3)
REAL, DIMENSION(:), INTENT(INOUT) :: PSMAX ! Maximum supersaturation
!
END SUBROUTINE CH_AER_ACTIVATION
42 changes: 33 additions & 9 deletions src/common/micro/lima.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -4,13 +4,14 @@
!MNH_LIC for details. version 1.
!-----------------------------------------------------------------
! #####################################################################
SUBROUTINE LIMA ( D, CST, ICED, ICEP, ELECD, ELECP, BUCONF, TBUDGETS, KBUDGETS, KRR, &
SUBROUTINE LIMA ( D, CST, ICED, ICEP, ELECD, ELECP, BUCONF, TBUDGETS, HACTCCN, KBUDGETS, KRR, &
PTSTEP, OELEC, &
PRHODREF, PEXNREF, PDZZ, PTHVREFZIKB, &
PRHODJ, PPABST, &
NCCN, NIFN, NIMM, &
NCCN, NIFN, NIMM, KCARB, KSOA, KSP, &
ODUST, OSALT, OORILAM, &
ODTHRAD, PDTHRAD, PTHT, PRT, PSVT, PW_NU, &
PTHS, PRS, PSVS, &
PAERO,PSOLORG, PMI,PTHS, PRS, PSVS, &
PINPRC, PINDEP, PINPRR, PINPRI, PINPRS, PINPRG, PINPRH, &
PEVAP3D, PCLDFR, PICEFR, PPRCFR, PFPR, &
PLATHAM_IAGGS, PEFIELDW, PSV_ELEC_T, PSV_ELEC_S )
Expand Down Expand Up @@ -65,7 +66,6 @@ SUBROUTINE LIMA ( D, CST, ICED, ICEP, ELECD, ELECP, BUCONF, TBUDGETS, KBUDGETS,
LFEEDBACKT, NMAXITER, XMRSTEP, XTSTEP_TS, &
LSEDC, LSEDI, XRTMIN, XCTMIN, LDEPOC, XVDEPOC, &
NMOM_C, NMOM_R, NMOM_I, NMOM_S, NMOM_G, NMOM_H

USE MODE_BUDGET_PHY, ONLY: BUDGET_STORE_ADD_PHY, BUDGET_STORE_INIT_PHY, BUDGET_STORE_END_PHY
use mode_tools, only: Countjv

Expand All @@ -89,6 +89,7 @@ SUBROUTINE LIMA ( D, CST, ICED, ICEP, ELECD, ELECP, BUCONF, TBUDGETS, KBUDGETS,
TYPE(ELEC_PARAM_t), INTENT(IN) :: ELECP ! electrical parameters
TYPE(ELEC_DESCR_t), INTENT(IN) :: ELECD ! electrical descriptive csts
TYPE(TBUDGETCONF_t), INTENT(IN) :: BUCONF
CHARACTER(LEN=4), INTENT(IN) :: HACTCCN ! kind of CCN activation
TYPE(TBUDGETDATA), DIMENSION(KBUDGETS), INTENT(INOUT) :: TBUDGETS
INTEGER, INTENT(IN) :: KBUDGETS
INTEGER, INTENT(IN) :: KRR
Expand All @@ -107,15 +108,20 @@ SUBROUTINE LIMA ( D, CST, ICED, ICEP, ELECD, ELECP, BUCONF, TBUDGETS, KBUDGETS,
INTEGER, INTENT(IN) :: NCCN ! for array size declarations
INTEGER, INTENT(IN) :: NIFN ! for array size declarations
INTEGER, INTENT(IN) :: NIMM ! for array size declarations
INTEGER, INTENT(IN) :: KCARB, KSOA, KSP ! for array size declarations
LOGICAL, INTENT(IN) :: ODUST, OSALT, OORILAM
!
LOGICAL, INTENT(IN) :: ODTHRAD ! Use radiative temperature tendency
REAL, DIMENSION(MERGE(D%NIT,0,ODTHRAD), &
MERGE(D%NJT,0,ODTHRAD), &
MERGE(D%NKT,0,ODTHRAD)), INTENT(IN) :: PDTHRAD ! dT/dt due to radiation
REAL, DIMENSION(D%NIT, D%NJT, D%NKT), INTENT(IN) :: PTHT ! Theta at time t
REAL, DIMENSION(D%NIT, D%NJT, D%NKT), INTENT(IN) :: PTHT ! Theta at time t
REAL, DIMENSION(D%NIT, D%NJT, D%NKT, KRR), INTENT(IN) :: PRT ! Mixing ratios at time t
REAL, DIMENSION(D%NIT, D%NJT, D%NKT, NSV), INTENT(IN) :: PSVT ! Concentrations at time t
REAL, DIMENSION(D%NIT, D%NJT, D%NKT), INTENT(IN) :: PW_NU ! w for CCN activation
REAL, DIMENSION(D%NIT, D%NJT, D%NKT), INTENT(IN) :: PW_NU ! w for CCN activation
REAL, DIMENSION(D%NIT, D%NJT, D%NKT ,NSV), INTENT(INOUT) :: PAERO ! Aerosol concentration
REAL, DIMENSION(D%NIT, D%NJT, D%NKT, 10), INTENT(IN) :: PSOLORG ![%] solubility fraction of soa
REAL, DIMENSION(D%NIT, D%NJT, D%NKT, KSP+KCARB+KSOA), INTENT(IN) :: PMI
!
REAL, DIMENSION(D%NIT, D%NJT, D%NKT), INTENT(INOUT) :: PTHS ! Theta source
REAL, DIMENSION(D%NIT, D%NJT, D%NKT, KRR), INTENT(INOUT) :: PRS ! Mixing ratios sources
Expand Down Expand Up @@ -1042,10 +1048,11 @@ SUBROUTINE LIMA ( D, CST, ICED, ICEP, ELECD, ELECP, BUCONF, TBUDGETS, KBUDGETS,
! --------------------
!
CALL LIMA_NUCLEATION_PROCS (D, CST, BUCONF, TBUDGETS, KBUDGETS, &
KCARB, KSOA, KSP, ODUST, OSALT, OORILAM, &
PTSTEP, PRHODJ, &
PRHODREF, ZEXN, PPABST, ZT, PDTHRAD, PW_NU, &
ZTHT, ZRVT, ZRCT, ZRRT, ZRIT, ZRST, ZRGT, ZRHT, &
ZCCT, ZCRT, ZCIT, &
ZCCT, ZCRT, ZCIT, PAERO,PSOLORG, PMI, HACTCCN, &
ZCCNFT, ZCCNAT, ZIFNFT, ZIFNNT, ZIMMNT, ZHOMFT, &
PCLDFR, PICEFR, PPRCFR, &
ZTOT_RV_HENU, ZTOT_RC_HINC, ZTOT_RI_HIND, ZTOT_RV_HONH)
Expand Down Expand Up @@ -1492,21 +1499,39 @@ SUBROUTINE LIMA ( D, CST, ICED, ICEP, ELECD, ELECP, BUCONF, TBUDGETS, KBUDGETS,
WHERE(ZA_RC(:)<-1.E-20 .AND. ZRCT1D(:)>XRTMIN(2))
ZMAXTIME(:)=MIN(ZMAXTIME(:), -(ZB_RC(:)+ZRCT1D(:))/ZA_RC(:))
END WHERE
WHERE(ZA_CC(:)<-1.E-20 .AND. ZCCT1D(:)>XCTMIN(2))
ZMAXTIME(:)=MIN(ZMAXTIME(:), -(ZB_CC(:)+ZCCT1D(:))/ZA_CC(:))
END WHERE
WHERE(ZA_RR(:)<-1.E-20 .AND. ZRRT1D(:)>XRTMIN(3))
ZMAXTIME(:)=MIN(ZMAXTIME(:), -(ZB_RR(:)+ZRRT1D(:))/ZA_RR(:))
END WHERE
WHERE(ZA_CR(:)<-1.E-20 .AND. ZCRT1D(:)>XCTMIN(3))
ZMAXTIME(:)=MIN(ZMAXTIME(:), -(ZB_CR(:)+ZCRT1D(:))/ZA_CR(:))
END WHERE
WHERE(ZA_RI(:)<-1.E-20 .AND. ZRIT1D(:)>XRTMIN(4))
ZMAXTIME(:)=MIN(ZMAXTIME(:), -(ZB_RI(:)+ZRIT1D(:))/ZA_RI(:))
END WHERE
WHERE(ZA_CI(:)<-1.E-20 .AND. ZCIT1D(:)>XCTMIN(4))
ZMAXTIME(:)=MIN(ZMAXTIME(:), -(ZB_CI(:)+ZCIT1D(:))/ZA_CI(:))
END WHERE
WHERE(ZA_RS(:)<-1.E-20 .AND. ZRST1D(:)>XRTMIN(5))
ZMAXTIME(:)=MIN(ZMAXTIME(:), -(ZB_RS(:)+ZRST1D(:))/ZA_RS(:))
END WHERE
WHERE(ZA_CS(:)<-1.E-20 .AND. ZCST1D(:)>XCTMIN(5))
ZMAXTIME(:)=MIN(ZMAXTIME(:), -(ZB_CS(:)+ZCST1D(:))/ZA_CS(:))
END WHERE
WHERE(ZA_RG(:)<-1.E-20 .AND. ZRGT1D(:)>XRTMIN(6))
ZMAXTIME(:)=MIN(ZMAXTIME(:), -(ZB_RG(:)+ZRGT1D(:))/ZA_RG(:))
END WHERE
WHERE(ZA_CG(:)<-1.E-20 .AND. ZCGT1D(:)>XCTMIN(6))
ZMAXTIME(:)=MIN(ZMAXTIME(:), -(ZB_CG(:)+ZCGT1D(:))/ZA_CG(:))
END WHERE
WHERE(ZA_RH(:)<-1.E-20 .AND. ZRHT1D(:)>XRTMIN(7))
ZMAXTIME(:)=MIN(ZMAXTIME(:), -(ZB_RH(:)+ZRHT1D(:))/ZA_RH(:))
END WHERE
WHERE(ZA_CH(:)<-1.E-20 .AND. ZCHT1D(:)>XCTMIN(7))
ZMAXTIME(:)=MIN(ZMAXTIME(:), -(ZB_CH(:)+ZCHT1D(:))/ZA_CH(:))
END WHERE

! We stop when the end of the timestep is reached
WHERE(PTSTEP-ZTIME1D(:)-ZMAXTIME(:)<=0.)
Expand Down Expand Up @@ -2297,8 +2322,7 @@ SUBROUTINE LIMA ( D, CST, ICED, ICEP, ELECD, ELECP, BUCONF, TBUDGETS, KBUDGETS,
call BUDGET_STORE_ADD_PHY(D, TBUDGETS(NBUDGET_RR), 'ACCR', -ztot_rc_accr(:, :, :) * zrhodjontstep(:, :, :) )
call BUDGET_STORE_ADD_PHY(D, TBUDGETS(NBUDGET_RR), 'REVA', ztot_rr_evap(:, :, :) * zrhodjontstep(:, :, :) )
call BUDGET_STORE_ADD_PHY(D, TBUDGETS(NBUDGET_RR), 'HONR', ztot_rr_honr(:, :, :) * zrhodjontstep(:, :, :) )
! call BUDGET_STORE_ADD_PHY(D, TBUDGETS(NBUDGET_RR), 'ACC', ztot_rr_acc (:, :, :) * zrhodjontstep(:, :, :) )
call BUDGET_STORE_ADD_PHY(D, TBUDGETS(NBUDGET_RR), 'ACC', (ztot_rc_rimss(:, :, :) + ztot_rc_rimsg(:, :, :)) &
call BUDGET_STORE_ADD_PHY(D, TBUDGETS(NBUDGET_RR), 'ACC', (-ztot_rr_accss(:, :, :) - ztot_rr_accsg(:, :, :)) &
* zrhodjontstep(:, :, :) )
call BUDGET_STORE_ADD_PHY(D, TBUDGETS(NBUDGET_RR), 'CFRZ', ztot_rr_cfrz(:, :, :) * zrhodjontstep(:, :, :) )
call BUDGET_STORE_ADD_PHY(D, TBUDGETS(NBUDGET_RR), 'WETG', ztot_rr_wetg(:, :, :) * zrhodjontstep(:, :, :) )
Expand Down
Loading

0 comments on commit 669c0ef

Please sign in to comment.