Changeset 242
- Timestamp:
- 07/21/14 12:01:59 (10 years ago)
- Location:
- codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/dyn3dpar/init_gcm_lmdz.F
r222 r242 512 512 ! & iflag_phys) 513 513 514 zcufi(nbp)=1.515 zcvfi(nbp)=1.516 517 CALL iniphysiq(klon_glo,llm,daysec,day_ini,dtphys/nsplit_phys,518 & latfi,lonfi,airefi,zcufi,zcvfi,rad,g,r,cpp,519 & iflag_phys)514 ! zcufi(nbp)=1. 515 ! zcvfi(nbp)=1. 516 517 ! CALL iniphysiq(klon_glo,llm,daysec,day_ini,dtphys/nsplit_phys, 518 ! & latfi,lonfi,airefi,zcufi,zcvfi,rad,g,r,cpp, 519 ! & iflag_phys) 520 520 521 521 #endif -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/comgeomphy.F90
r222 r242 1 1 module comgeomphy 2 2 real,save,allocatable :: airephy(:) 3 real,save,allocatable :: cuphy(:)4 real,save,allocatable :: cvphy(:)5 3 real,save,allocatable :: rlatd(:) 6 4 real,save,allocatable :: rlond(:) … … 14 12 15 13 allocate(airephy(klon_omp)) 16 allocate(cuphy(klon_omp))17 allocate(cvphy(klon_omp))18 14 allocate(rlatd(klon_omp)) 19 15 allocate(rlond(klon_omp)) -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/iniphysiq.F90
r227 r242 1 subroutine iniphysiq(ngrid,nlayer, punjours, pdayref,ptimestep, & 2 plat,plon,parea,pcu,pcv, & 3 prad,pg,pr,pcpp,iflag_phys) 1 subroutine iniphysiq(nlayer, punjours, pdayref,ptimestep, & 2 plat,plon,parea, & 3 prad,pg,pr,pcpp, & 4 preff,ap,bp) 4 5 5 6 use dimphy, only : klev ! number of atmospheric levels … … 12 13 klon_mpi 13 14 use comgeomphy, only : airephy, & ! physics grid area (m2) 14 cuphy, & ! cu coeff. (u_covariant = cu * u)15 cvphy, & ! cv coeff. (v_covariant = cv * v)16 15 rlond, & ! longitudes 17 16 rlatd ! latitudes … … 20 19 21 20 implicit none 22 include "dimensions.h"23 include "comvert.h"24 21 25 22 real,intent(in) :: prad ! radius of the planet (m) … … 28 25 real,intent(in) :: pcpp ! specific heat Cp 29 26 real,intent(in) :: punjours ! length (in s) of a standard day 30 integer,intent(in) :: ngrid ! number of horizontal grid points in the physics (full grid)31 27 integer,intent(in) :: nlayer ! number of atmospheric layers 32 28 real,intent(in) :: plat(klon_mpi) ! latitudes of the physics grid 33 29 real,intent(in) :: plon(klon_mpi) ! longitudes of the physics grid 34 30 real,intent(in) :: parea(klon_mpi) ! area (m2) 35 real,intent(in) :: pcu(klon_mpi) ! cu coeff. (u_covariant = cu * u)36 real,intent(in) :: pcv(klon_mpi) ! cv coeff. (v_covariant = cv * v)37 31 integer,intent(in) :: pdayref ! reference day of for the simulation 38 real,intent(in) :: ptimestep !physics time step (s) 39 integer,intent(in) :: iflag_phys ! type of physics to be called 32 real,intent(in) :: ptimestep ! physics time step (s) 33 real,intent(in) :: preff ! reference surface pressure (Pa) 34 real,intent(in) :: ap(nlayer+1) ! hybrid coordinate at interfaces 35 real,intent(in) :: bp(nlayer+1) ! hybrid coordinate at interfaces 40 36 41 integer :: ibegin,iend,offset37 integer :: offset 42 38 character(len=20) :: modname='iniphysiq' 43 39 character(len=80) :: abort_message … … 52 48 ENDIF 53 49 54 IF (ngrid.NE.klon_glo) THEN55 write(*,*) 'STOP in ',trim(modname)56 write(*,*) 'Problem with dimensions :'57 write(*,*) 'ngrid = ',ngrid58 write(*,*) 'klon = ',klon_glo59 abort_message = ''60 CALL abort_gcm (modname,abort_message,1)61 ENDIF62 50 63 !$OMP PARALLEL PRIVATE(ibegin,iend) & 64 !$OMP SHARED(parea,pcu,pcv,plon,plat) 51 !$OMP PARALLEL SHARED(parea,plon,plat) 65 52 66 53 offset=0 67 54 airephy(1:klon_omp)=parea(offset+klon_omp_begin:offset+klon_omp_end) 68 cuphy(1:klon_omp)=pcu(offset+klon_omp_begin:offset+klon_omp_end)69 cvphy(1:klon_omp)=pcv(offset+klon_omp_begin:offset+klon_omp_end)70 55 rlond(1:klon_omp)=plon(offset+klon_omp_begin:offset+klon_omp_end) 71 56 rlatd(1:klon_omp)=plat(offset+klon_omp_begin:offset+klon_omp_end)
Note: See TracChangeset
for help on using the changeset viewer.