| 1 | subroutine iniphysiq(ngrid,nlayer, punjours, pdayref,ptimestep, & |
|---|
| 2 | plat,plon,parea,pcu,pcv, & |
|---|
| 3 | prad,pg,pr,pcpp,iflag_phys) |
|---|
| 4 | |
|---|
| 5 | use dimphy, only : klev ! number of atmospheric levels |
|---|
| 6 | use mod_grid_phy_lmdz, only : klon_glo ! number of atmospheric columns |
|---|
| 7 | ! (on full grid) |
|---|
| 8 | use mod_phys_lmdz_para, only : klon_omp, & ! number of columns (on local omp grid) |
|---|
| 9 | klon_omp_begin, & ! start index of local omp subgrid |
|---|
| 10 | klon_omp_end, & ! end index of local omp subgrid |
|---|
| 11 | klon_mpi_begin ! start indes of columns (on local mpi grid) |
|---|
| 12 | use comgeomphy, only : airephy, & ! physics grid area (m2) |
|---|
| 13 | cuphy, & ! cu coeff. (u_covariant = cu * u) |
|---|
| 14 | cvphy, & ! cv coeff. (v_covariant = cv * v) |
|---|
| 15 | rlond, & ! longitudes |
|---|
| 16 | rlatd ! latitudes |
|---|
| 17 | use infotrac, only : nqtot ! number of advected tracers |
|---|
| 18 | |
|---|
| 19 | implicit none |
|---|
| 20 | |
|---|
| 21 | real,intent(in) :: prad ! radius of the planet (m) |
|---|
| 22 | real,intent(in) :: pg ! gravitational acceleration (m/s2) |
|---|
| 23 | real,intent(in) :: pr ! ! reduced gas constant R/mu |
|---|
| 24 | real,intent(in) :: pcpp ! specific heat Cp |
|---|
| 25 | real,intent(in) :: punjours ! length (in s) of a standard day |
|---|
| 26 | integer,intent(in) :: ngrid ! number of horizontal grid points in the physics (full grid) |
|---|
| 27 | integer,intent(in) :: nlayer ! number of atmospheric layers |
|---|
| 28 | real,intent(in) :: plat(ngrid) ! latitudes of the physics grid |
|---|
| 29 | real,intent(in) :: plon(ngrid) ! longitudes of the physics grid |
|---|
| 30 | real,intent(in) :: parea(klon_glo) ! area (m2) |
|---|
| 31 | real,intent(in) :: pcu(klon_glo) ! cu coeff. (u_covariant = cu * u) |
|---|
| 32 | real,intent(in) :: pcv(klon_glo) ! cv coeff. (v_covariant = cv * v) |
|---|
| 33 | integer,intent(in) :: pdayref ! reference day of for the simulation |
|---|
| 34 | real,intent(in) :: ptimestep !physics time step (s) |
|---|
| 35 | integer,intent(in) :: iflag_phys ! type of physics to be called |
|---|
| 36 | |
|---|
| 37 | integer :: ibegin,iend,offset |
|---|
| 38 | character(len=20) :: modname='iniphysiq' |
|---|
| 39 | character(len=80) :: abort_message |
|---|
| 40 | |
|---|
| 41 | IF (nlayer.NE.klev) THEN |
|---|
| 42 | write(*,*) 'STOP in ',trim(modname) |
|---|
| 43 | write(*,*) 'Problem with dimensions :' |
|---|
| 44 | write(*,*) 'nlayer = ',nlayer |
|---|
| 45 | write(*,*) 'klev = ',klev |
|---|
| 46 | abort_message = '' |
|---|
| 47 | CALL abort_gcm (modname,abort_message,1) |
|---|
| 48 | ENDIF |
|---|
| 49 | |
|---|
| 50 | IF (ngrid.NE.klon_glo) THEN |
|---|
| 51 | write(*,*) 'STOP in ',trim(modname) |
|---|
| 52 | write(*,*) 'Problem with dimensions :' |
|---|
| 53 | write(*,*) 'ngrid = ',ngrid |
|---|
| 54 | write(*,*) 'klon = ',klon_glo |
|---|
| 55 | abort_message = '' |
|---|
| 56 | CALL abort_gcm (modname,abort_message,1) |
|---|
| 57 | ENDIF |
|---|
| 58 | |
|---|
| 59 | !$OMP PARALLEL PRIVATE(ibegin,iend) |
|---|
| 60 | !$OMP+ SHARED(parea,pcu,pcv,plon,plat) |
|---|
| 61 | |
|---|
| 62 | offset=klon_mpi_begin-1 |
|---|
| 63 | airephy(1:klon_omp)=parea(offset+klon_omp_begin:offset+klon_omp_end) |
|---|
| 64 | cuphy(1:klon_omp)=pcu(offset+klon_omp_begin:offset+klon_omp_end) |
|---|
| 65 | cvphy(1:klon_omp)=pcv(offset+klon_omp_begin:offset+klon_omp_end) |
|---|
| 66 | rlond(1:klon_omp)=plon(offset+klon_omp_begin:offset+klon_omp_end) |
|---|
| 67 | rlatd(1:klon_omp)=plat(offset+klon_omp_begin:offset+klon_omp_end) |
|---|
| 68 | |
|---|
| 69 | ! copy some fundamental parameters to physics |
|---|
| 70 | ! and do some initializations |
|---|
| 71 | call inifis(klon_omp,nlayer,nqtot,pdayref,punjours,ptimestep, & |
|---|
| 72 | rlatd,rlond,airephy,prad,pg,pr,pcpp) |
|---|
| 73 | |
|---|
| 74 | !$OMP END PARALLEL |
|---|
| 75 | |
|---|
| 76 | |
|---|
| 77 | end subroutine iniphysiq |
|---|
