Directory: | ./ |
---|---|
File: | dyn_phys_sub/callphysiq_mod.f90 |
Date: | 2022-01-11 19:19:34 |
Exec | Total | Coverage | |
---|---|---|---|
Lines: | 6 | 7 | 85.7% |
Branches: | 3 | 4 | 75.0% |
Line | Branch | Exec | Source |
---|---|---|---|
1 | ! | ||
2 | ! $Id: $ | ||
3 | ! | ||
4 | MODULE callphysiq_mod | ||
5 | |||
6 | IMPLICIT NONE | ||
7 | |||
8 | CONTAINS | ||
9 | |||
10 | 480 | SUBROUTINE call_physiq(klon,llm,nqtot,tname, & | |
11 | debut_split,lafin_split, & | ||
12 | jD_cur,jH_cur_split,zdt_split, & | ||
13 | zplev_omp,zplay_omp, & | ||
14 | zpk_omp,zphi_omp,zphis_omp, & | ||
15 | presnivs_omp, & | ||
16 | zufi_omp,zvfi_omp,zrfi_omp,ztfi_omp,zqfi_omp, & | ||
17 | flxwfi_omp,pducov, & | ||
18 | zdufi_omp,zdvfi_omp,zdtfi_omp,zdqfi_omp,zdpsrf_omp) | ||
19 | |||
20 | |||
21 | USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat | ||
22 | USE control_mod, ONLY: planet_type | ||
23 | USE physiq_mod, ONLY: physiq | ||
24 | IMPLICIT NONE | ||
25 | |||
26 | INTEGER,INTENT(IN) :: klon ! (local) number of atmospheric columns | ||
27 | INTEGER,INTENT(IN) :: llm ! number of atmospheric layers | ||
28 | INTEGER,INTENT(IN) :: nqtot ! number of tracers | ||
29 | CHARACTER(len=*),INTENT(IN) :: tname(nqtot) ! tracer names | ||
30 | LOGICAL,INTENT(IN) :: debut_split ! .true. if very first call to physics | ||
31 | LOGICAL,INTENT(IN) :: lafin_split ! .true. if last call to physics | ||
32 | REAL,INTENT(IN) :: JD_cur ! Julian day | ||
33 | REAL,INTENT(IN) :: JH_cur_split ! Julian hour (fraction of day) | ||
34 | REAL,INTENT(IN) :: zdt_split ! time step over which the physics are evaluated | ||
35 | REAL,INTENT(IN) :: zplev_omp(klon,llm+1) ! interlayer pressure (Pa) | ||
36 | REAL,INTENT(IN) :: zplay_omp(klon,llm) ! mid-layer pressure (Pa) | ||
37 | REAL,INTENT(IN) :: zpk_omp(klon,llm) ! Exner function | ||
38 | REAL,INTENT(IN) :: zphi_omp(klon,llm) ! geopotential at midlayer | ||
39 | REAL,INTENT(IN) :: zphis_omp(klon) ! surface geopotential | ||
40 | REAL,INTENT(IN) :: presnivs_omp(llm) ! approximate pressure of atm. layers | ||
41 | REAL,INTENT(IN) :: zufi_omp(klon,llm) ! zonal wind (m/s) | ||
42 | REAL,INTENT(IN) :: zvfi_omp(klon,llm) ! meridional wind (m/s) | ||
43 | REAL,INTENT(IN) :: zrfi_omp(klon,llm) ! relative wind vorticity, in s-1 | ||
44 | REAL,INTENT(IN) :: ztfi_omp(klon,llm) ! temperature (K) | ||
45 | REAL,INTENT(IN) :: zqfi_omp(klon,llm,nqtot) ! tracers (*/kg of air) | ||
46 | REAL,INTENT(IN) :: flxwfi_omp(klon,llm) ! Vertical mass flux on lower mesh interfaces (kg/s) | ||
47 | REAL,INTENT(IN) :: pducov(nbp_lon+1,nbp_lat,llm) ! dynamical tendency on ucov | ||
48 | ! tendencies (in */s) from the physics: | ||
49 | REAL,INTENT(OUT) :: zdufi_omp(klon,llm) ! tendency on zonal winds | ||
50 | REAL,INTENT(OUT) :: zdvfi_omp(klon,llm) ! tendency on meridional winds | ||
51 | REAL,INTENT(OUT) :: zdtfi_omp(klon,llm) ! tendency on temperature | ||
52 | REAL,INTENT(OUT) :: zdqfi_omp(klon,llm,nqtot) ! tendency on tracers | ||
53 | REAL,INTENT(OUT) :: zdpsrf_omp(klon) ! tendency on surface pressure | ||
54 | |||
55 | ! Local variables | ||
56 | CHARACTER(len=11) :: modname="call_physiq" | ||
57 | LOGICAL,SAVE :: firstcall=.true. | ||
58 | !$OMP THREADPRIVATE(firstcall) | ||
59 | |||
60 | ! Sanity check on physics package type | ||
61 |
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 479 times.
|
480 | IF (firstcall) THEN |
62 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (planet_type.ne."earth") THEN |
63 | ✗ | CALL abort_gcm(modname,"wrong planet_type for this physics package",1) | |
64 | ENDIF | ||
65 | 1 | firstcall=.false. | |
66 | ENDIF | ||
67 | |||
68 | |||
69 | ! Call physics package with required inputs/outputs | ||
70 | CALL physiq(klon, & | ||
71 | llm, & | ||
72 | debut_split, & | ||
73 | lafin_split, & | ||
74 | zdt_split, & | ||
75 | zplev_omp, & | ||
76 | zplay_omp, & | ||
77 | zphi_omp, & | ||
78 | zphis_omp, & | ||
79 | presnivs_omp, & | ||
80 | zufi_omp, & | ||
81 | zvfi_omp, & | ||
82 | zrfi_omp, & | ||
83 | ztfi_omp, & | ||
84 | zqfi_omp, & | ||
85 | flxwfi_omp, & | ||
86 | zdufi_omp, & | ||
87 | zdvfi_omp, & | ||
88 | zdtfi_omp, & | ||
89 | zdqfi_omp, & | ||
90 | 480 | zdpsrf_omp) | |
91 | |||
92 | |||
93 | 480 | END SUBROUTINE call_physiq | |
94 | |||
95 | END MODULE callphysiq_mod | ||
96 |