/[lmdze]/trunk/phylmd/Orography/start_init_orog.f
ViewVC logotype

Diff of /trunk/phylmd/Orography/start_init_orog.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC revision 225 by guez, Mon Oct 16 12:35:41 2017 UTC
# Line 17  CONTAINS Line 17  CONTAINS
17         zthe_2d, zpic_2d, zval_2d)         zthe_2d, zpic_2d, zval_2d)
18    
19      use conf_dat2d_m, only: conf_dat2d      use conf_dat2d_m, only: conf_dat2d
20      use comgeom, only: rlatu, rlonv      use dynetat0_m, only: rlatu, rlonv
21      use grid_noro_m, only: grid_noro      use grid_noro_m, only: grid_noro
22      use indicesol, only: epsfra      use indicesol, only: epsfra
23      use netcdf, only: nf90_nowrite      use netcdf, only: nf90_nowrite
# Line 53  CONTAINS Line 53  CONTAINS
53      INTEGER iml_rel      INTEGER iml_rel
54      INTEGER jml_rel      INTEGER jml_rel
55      INTEGER ncid, varid      INTEGER ncid, varid
56      REAL, pointer:: relief(:, :)      REAL, ALLOCATABLE:: relief(:, :)
57      REAL, ALLOCATABLE:: lon_rad(:), lat_rad(:)      REAL, ALLOCATABLE:: lon_rad(:), lat_rad(:)
58      REAL, pointer:: lon_ini(:), lat_ini(:)      REAL, ALLOCATABLE:: lon_ini(:), lat_ini(:)
59    
60      !-----------------------------------      !-----------------------------------
61    
# Line 92  CONTAINS Line 92  CONTAINS
92    
93      CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, relief , &      CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, relief , &
94           interbar=.FALSE.)           interbar=.FALSE.)
     deallocate(lon_ini, lat_ini) ! pointers  
95    
96      print *, 'Compute all the parameters needed for the gravity wave drag code'      print *, 'Compute all the parameters needed for the gravity wave drag code'
97    
98      CALL grid_noro(lon_rad, lat_rad, relief, rlonv, rlatu, phis, zmea_2d, &      CALL grid_noro(lon_rad, lat_rad, relief, rlonv, rlatu, phis, zmea_2d, &
99           zstd_2d, zsig_2d, zgam_2d, zthe_2d, zpic_2d, zval_2d, mask)           zstd_2d, zsig_2d, zgam_2d, zthe_2d, zpic_2d, zval_2d, mask)
     deallocate(relief) ! pointer  
100      phis(:, :) = phis(:, :) * 9.81      phis(:, :) = phis(:, :) * 9.81
101    
102      mask(2:, 1) = mask(1, 1) ! north pole      mask(2:, 1) = mask(1, 1) ! north pole
Legend:
Removed from v.134  
changed lines
  Added in v.225
  ViewVC Help
Powered by ViewVC 1.1.21