| [3661] | 1 | ! |
|---|
| 2 | ! $Id: $ |
|---|
| 3 | ! |
|---|
| 4 | MODULE callphysiq_mod |
|---|
| 5 | |
|---|
| 6 | IMPLICIT NONE |
|---|
| 7 | |
|---|
| 8 | CONTAINS |
|---|
| 9 | |
|---|
| 10 | SUBROUTINE call_physiq(planet_type,klon,llm,nqtot, & |
|---|
| 11 | debut_split,lafin_split) |
|---|
| 12 | |
|---|
| 13 | USE variables_mod |
|---|
| 14 | USE physiq_mod, ONLY: physiq |
|---|
| 15 | use phys_state_var_mod , only : qsurf |
|---|
| 16 | use comm_wrf, only : allocate_comm_wrf |
|---|
| 17 | use tracer_h, only: noms |
|---|
| 18 | |
|---|
| 19 | IMPLICIT NONE |
|---|
| 20 | character(len=10),INTENT(IN) :: planet_type ! planet type ('earth','mars',...) |
|---|
| 21 | INTEGER,INTENT(IN) :: klon ! (local) number of atmospheric columns |
|---|
| 22 | INTEGER,INTENT(IN) :: llm ! number of atmospheric layers |
|---|
| 23 | INTEGER,INTENT(IN) :: nqtot ! number of tracers |
|---|
| 24 | LOGICAL,INTENT(IN) :: debut_split ! .true. if very first call to physics |
|---|
| 25 | LOGICAL,INTENT(IN) :: lafin_split ! .true. if last call to physics |
|---|
| 26 | |
|---|
| 27 | ! Local variables |
|---|
| 28 | ! CHARACTER(len=11) :: modname="call_physiq" |
|---|
| 29 | |
|---|
| 30 | ! Sanity check on physics package type |
|---|
| 31 | IF (debut_split) THEN |
|---|
| 32 | IF (planet_type.ne."titan") THEN |
|---|
| 33 | PRINT*,"wrong planet_type for this physics package" |
|---|
| 34 | STOP |
|---|
| 35 | ENDIF |
|---|
| 36 | ENDIF |
|---|
| 37 | |
|---|
| 38 | !qsurf(:,:)=zqfi_omp(:,1,:) |
|---|
| 39 | !JL22 what is this, why in the world would we want to put the value of the last layer at the surface |
|---|
| 40 | ! these two quantities do not even have the same units (kg/kg and kg/m^2)!!! |
|---|
| 41 | ! seems it is the same for some other physics. It is probably a bug. |
|---|
| 42 | !EMo25 zqfi_omp in kg/kghumidair? |
|---|
| 43 | |
|---|
| 44 | call allocate_comm_wrf(klon,llm,nqtot) |
|---|
| 45 | ! Call physics package with required inputs/outputs |
|---|
| 46 | CALL physiq(klon, & ! ngrid |
|---|
| 47 | llm, & ! nlayer |
|---|
| 48 | nqtot, & ! nq |
|---|
| 49 | noms, & ! nametrac |
|---|
| 50 | debut_split, & ! firstcall |
|---|
| 51 | lafin_split, & ! lastcall |
|---|
| 52 | jD_cur, & ! pday |
|---|
| 53 | jH_cur_split, & ! ptime |
|---|
| 54 | zdt_split, & ! ptimestep |
|---|
| 55 | zplev_omp, & ! pplev |
|---|
| 56 | zplay_omp, & ! pplay |
|---|
| 57 | zphi_omp, & ! pphi |
|---|
| 58 | zufi_omp, & ! pu |
|---|
| 59 | zvfi_omp, & ! pv |
|---|
| 60 | ztfi_omp, & ! pt |
|---|
| 61 | zqfi_omp, & ! pq |
|---|
| 62 | flxwfi_omp, & ! flxw |
|---|
| 63 | zdufi_omp, & ! pdu |
|---|
| 64 | zdvfi_omp, & ! pdv |
|---|
| 65 | zdtfi_omp, & ! pdt |
|---|
| 66 | zdqfi_omp, & ! pdq |
|---|
| 67 | zdpsrf_omp) ! pdpsrf |
|---|
| 68 | |
|---|
| 69 | |
|---|
| 70 | END SUBROUTINE call_physiq |
|---|
| 71 | |
|---|
| 72 | END MODULE callphysiq_mod |
|---|
