Skip to content
Snippets Groups Projects
Commit b05e7f9d authored by RODIER Quentin's avatar RODIER Quentin
Browse files

Quentin 28/02/2022: Minor, adapt mf_turb_greyzone to reversed args in shuman

parent a76a1da9
No related branches found
No related tags found
No related merge requests found
...@@ -264,7 +264,7 @@ CALL TRIDIAG_MASSFLUX(KKA,KKB,KKE,KKU,KKL,PTHLM,PFLXZTHMF,-PEMF,PTSTEP,PIMPL, & ...@@ -264,7 +264,7 @@ CALL TRIDIAG_MASSFLUX(KKA,KKB,KKE,KKU,KKL,PTHLM,PFLXZTHMF,-PEMF,PTSTEP,PIMPL, &
! Est ce que ZVARS est au point de masse pour qu'il doivent être remis au point ! Est ce que ZVARS est au point de masse pour qu'il doivent être remis au point
! de flux ? ! de flux ?
!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!
PFLXZTHMF(:,:) = PEMF(:,:)*(PTHL_UP(:,:)-MZM_MF(KKA,KKU,KKL,ZVARS(:,:))) PFLXZTHMF(:,:) = PEMF(:,:)*(PTHL_UP(:,:)-MZM_MF(ZVARS(:,:),KKA,KKU,KKL))
!!! compute THL tendency !!! compute THL tendency
! !
...@@ -276,7 +276,7 @@ PTHLDT(:,:)= (ZVARS(:,:)-PTHLM(:,:))/PTSTEP ...@@ -276,7 +276,7 @@ PTHLDT(:,:)= (ZVARS(:,:)-PTHLM(:,:))/PTSTEP
CALL TRIDIAG_MASSFLUX(KKA,KKB,KKE,KKU,KKL,PRTM(:,:),PFLXZRMF,-PEMF,PTSTEP,PIMPL, & CALL TRIDIAG_MASSFLUX(KKA,KKB,KKE,KKU,KKL,PRTM(:,:),PFLXZRMF,-PEMF,PTSTEP,PIMPL, &
PDZZ,PRHODJ,ZVARS ) PDZZ,PRHODJ,ZVARS )
! compute new flux ! compute new flux
PFLXZRMF(:,:) = PEMF(:,:)*(PRT_UP(:,:)-MZM_MF(KKA,KKU,KKL,ZVARS(:,:))) PFLXZRMF(:,:) = PEMF(:,:)*(PRT_UP(:,:)-MZM_MF(ZVARS(:,:),KKA,KKU,KKL))
!!! compute RT tendency !!! compute RT tendency
PRTDT(:,:) = (ZVARS(:,:)-PRTM(:,:))/PTSTEP PRTDT(:,:) = (ZVARS(:,:)-PRTM(:,:))/PTSTEP
...@@ -291,7 +291,7 @@ IF (OMIXUV) THEN ...@@ -291,7 +291,7 @@ IF (OMIXUV) THEN
CALL TRIDIAG_MASSFLUX(KKA,KKB,KKE,KKU,KKL,PUM,PFLXZUMF,-PEMF,PTSTEP,PIMPL, & CALL TRIDIAG_MASSFLUX(KKA,KKB,KKE,KKU,KKL,PUM,PFLXZUMF,-PEMF,PTSTEP,PIMPL, &
PDZZ,PRHODJ,ZVARS ) PDZZ,PRHODJ,ZVARS )
! compute new flux ! compute new flux
PFLXZUMF(:,:) = PEMF(:,:)*(PU_UP(:,:)-MZM_MF(KKA,KKU,KKL,ZVARS(:,:))) PFLXZUMF(:,:) = PEMF(:,:)*(PU_UP(:,:)-MZM_MF(ZVARS(:,:),KKA,KKU,KKL))
! compute U tendency ! compute U tendency
PUDT(:,:)= (ZVARS(:,:)-PUM(:,:))/PTSTEP PUDT(:,:)= (ZVARS(:,:)-PUM(:,:))/PTSTEP
...@@ -305,7 +305,7 @@ IF (OMIXUV) THEN ...@@ -305,7 +305,7 @@ IF (OMIXUV) THEN
CALL TRIDIAG_MASSFLUX(KKA,KKB,KKE,KKU,KKL,PVM,PFLXZVMF,-PEMF,PTSTEP,PIMPL, & CALL TRIDIAG_MASSFLUX(KKA,KKB,KKE,KKU,KKL,PVM,PFLXZVMF,-PEMF,PTSTEP,PIMPL, &
PDZZ,PRHODJ,ZVARS ) PDZZ,PRHODJ,ZVARS )
! compute new flux ! compute new flux
PFLXZVMF(:,:) = PEMF(:,:)*(PV_UP(:,:)-MZM_MF(KKA,KKU,KKL,ZVARS(:,:))) PFLXZVMF(:,:) = PEMF(:,:)*(PV_UP(:,:)-MZM_MF(ZVARS(:,:),KKA,KKU,KKL))
! compute V tendency ! compute V tendency
PVDT(:,:)= (ZVARS(:,:)-PVM(:,:))/PTSTEP PVDT(:,:)= (ZVARS(:,:)-PVM(:,:))/PTSTEP
...@@ -321,7 +321,7 @@ DO JSV=1,ISV ...@@ -321,7 +321,7 @@ DO JSV=1,ISV
!* compute mean flux of scalar variables at time t-dt !* compute mean flux of scalar variables at time t-dt
! ( Resulting fluxes are in flux level (w-point) as PEMF and PTHL_UP ) ! ( Resulting fluxes are in flux level (w-point) as PEMF and PTHL_UP )
PFLXZSVMF(:,:,JSV) = PEMF(:,:)*(PSV_UP(:,:,JSV)-MZM_MF(KKA,KKU,KKL,PSVM(:,:,JSV))) PFLXZSVMF(:,:,JSV) = PEMF(:,:)*(PSV_UP(:,:,JSV)-MZM_MF(PSVM(:,:,JSV),KKA,KKU,KKL))
! !
! 3.5 Compute the tendency for scalar variables ! 3.5 Compute the tendency for scalar variables
...@@ -330,7 +330,7 @@ DO JSV=1,ISV ...@@ -330,7 +330,7 @@ DO JSV=1,ISV
CALL TRIDIAG_MASSFLUX(KKA,KKB,KKE,KKU,KKL,PSVM(:,:,JSV),PFLXZSVMF(:,:,JSV),& CALL TRIDIAG_MASSFLUX(KKA,KKB,KKE,KKU,KKL,PSVM(:,:,JSV),PFLXZSVMF(:,:,JSV),&
-PEMF,PTSTEP,PIMPL,PDZZ,PRHODJ,ZVARS ) -PEMF,PTSTEP,PIMPL,PDZZ,PRHODJ,ZVARS )
! compute new flux ! compute new flux
PFLXZSVMF(:,:,JSV) = PEMF(:,:)*(PSV_UP(:,:,JSV)-MZM_MF(KKA,KKU,KKL,ZVARS)) PFLXZSVMF(:,:,JSV) = PEMF(:,:)*(PSV_UP(:,:,JSV)-MZM_MF(ZVARS,KKA,KKU,KKL))
! compute Sv tendency ! compute Sv tendency
PSVDT(:,:,JSV)= (ZVARS(:,:)-PSVM(:,:,JSV))/PTSTEP PSVDT(:,:,JSV)= (ZVARS(:,:)-PSVM(:,:,JSV))/PTSTEP
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment