| 1 | SUBROUTINE diag_slp(nlon,t,pab,pal,pphis,tasfc,tastd,pmer) |
|---|
| 2 | USE dimphy |
|---|
| 3 | USE phys_output_write_mod |
|---|
| 4 | USE phys_output_ctrlout_mod |
|---|
| 5 | USE phys_local_var_mod |
|---|
| 6 | IMPLICIT NONE |
|---|
| 7 | !>====================================================================== |
|---|
| 8 | !! |
|---|
| 9 | !! Auteur(s) I.Musat (LMD/CNRS) date: 20151106 |
|---|
| 10 | !! |
|---|
| 11 | !! Objet: Calcul pression au niveau de la mer cf. Arpege-IFS |
|---|
| 12 | !! ctstar: calcule la temperature de l'air a la surface (tasfc) et |
|---|
| 13 | !! la temperature de l'air standard a la surface (tastd) |
|---|
| 14 | !! pppmer: calcule la slp a partir de tasfc, tastd, de la pression a la surface (pab1) |
|---|
| 15 | !! et du geopotentiel de la surface |
|---|
| 16 | !!====================================================================== |
|---|
| 17 | !! nlon--input-R-temperature au milieu de chaque couche (en K) |
|---|
| 18 | !! t--input-R-temperature au milieu de chaque couche (en K) |
|---|
| 19 | !! pab--input-R-pression pour chaque inter-couche (en Pa) |
|---|
| 20 | !! pal---input-R-pression pour le mileu de chaque couche (en Pa) |
|---|
| 21 | !! pphis---input-R-geopotentiel du sol (en m2/s2) |
|---|
| 22 | !! tasfc---output-R-temperature air au sol (en K) |
|---|
| 23 | !! tastd---output-R-temperature air 'standard' au sol (en K) |
|---|
| 24 | !! pmer---output-R-pression au niveau de la mer (en Pa) |
|---|
| 25 | !!====================================================================== |
|---|
| 26 | INTEGER nlon |
|---|
| 27 | REAL t(nlon,klev) |
|---|
| 28 | REAL pab(nlon,klev+1) |
|---|
| 29 | REAL pal(nlon,klev) |
|---|
| 30 | REAL pphis(nlon) |
|---|
| 31 | REAL tasfc(nlon), tastd(nlon) |
|---|
| 32 | REAL pmer(nlon) |
|---|
| 33 | !!! |
|---|
| 34 | !!! calcul tasfc et tastd |
|---|
| 35 | tal1(:)=t(:,1) |
|---|
| 36 | pab2(:)=pab(:,2) |
|---|
| 37 | pal1(:)=pal(:,1) |
|---|
| 38 | CALL ctstar(nlon,1,nlon,tal1,pab2,pal1,pphis,tasfc,tastd) |
|---|
| 39 | !!! |
|---|
| 40 | !!! calcul slp |
|---|
| 41 | pab1(:)=pab(:,1) |
|---|
| 42 | CALL pppmer(nlon,1,nlon,pab1,pphis,tasfc,tastd,pmer) |
|---|
| 43 | ! |
|---|
| 44 | RETURN |
|---|
| 45 | END |
|---|
