| 1 | ! |
|---|
| 2 | ! $Id: $ |
|---|
| 3 | ! |
|---|
| 4 | MODULE callphysiq_mod |
|---|
| 5 | |
|---|
| 6 | IMPLICIT NONE |
|---|
| 7 | |
|---|
| 8 | CONTAINS |
|---|
| 9 | |
|---|
| 10 | 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 | IF (firstcall) THEN |
|---|
| 62 | IF (planet_type.ne."earth") THEN |
|---|
| 63 | CALL abort_gcm(modname,"wrong planet_type for this physics package",1) |
|---|
| 64 | ENDIF |
|---|
| 65 | firstcall=.false. |
|---|
| 66 | ENDIF |
|---|
| 67 | |
|---|
| 68 | |
|---|
| 69 | ! Call physics package with required inputs/outputs |
|---|
| 70 | !#ifdef ISOVERIF |
|---|
| 71 | ! write(*,*) 'callphysiq 71: zqfi_omp(1,1,:)=',zqfi_omp(1,1,:) |
|---|
| 72 | !#endif |
|---|
| 73 | CALL physiq(klon, & |
|---|
| 74 | llm, & |
|---|
| 75 | debut_split, & |
|---|
| 76 | lafin_split, & |
|---|
| 77 | zdt_split, & |
|---|
| 78 | zplev_omp, & |
|---|
| 79 | zplay_omp, & |
|---|
| 80 | zphi_omp, & |
|---|
| 81 | zphis_omp, & |
|---|
| 82 | presnivs_omp, & |
|---|
| 83 | zufi_omp, & |
|---|
| 84 | zvfi_omp, & |
|---|
| 85 | zrfi_omp, & |
|---|
| 86 | ztfi_omp, & |
|---|
| 87 | zqfi_omp, & |
|---|
| 88 | flxwfi_omp, & |
|---|
| 89 | zdufi_omp, & |
|---|
| 90 | zdvfi_omp, & |
|---|
| 91 | zdtfi_omp, & |
|---|
| 92 | zdqfi_omp, & |
|---|
| 93 | zdpsrf_omp) |
|---|
| 94 | !#ifdef ISOVERIF |
|---|
| 95 | ! write(*,*) 'callphysiq 95: zdqfi_omp(1,1,:)=',zdqfi_omp(1,1,:) |
|---|
| 96 | !#endif |
|---|
| 97 | |
|---|
| 98 | |
|---|
| 99 | END SUBROUTINE call_physiq |
|---|
| 100 | |
|---|
| 101 | END MODULE callphysiq_mod |
|---|
