/[lmdze]/trunk/Sources/phylmd/Orography/start_init_orog_m.f
ViewVC logotype

Annotation of /trunk/Sources/phylmd/Orography/start_init_orog_m.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 139 - (hide annotations)
Tue May 26 17:46:03 2015 UTC (9 years ago) by guez
File size: 3626 byte(s) 
dynetat0 read rlonu, rlatu, rlonv, rlatv, cu_2d, cv_2d, aire_2d from
"start.nc" and then these variables were overwritten by
inigeom. Corrected this. Now, inigeom does not compute rlonu, rlatu,
rlonv and rlatv. Moreover, cu_2d, cv_2d, aire_2d are not written to
"restart.nc". Since xprimu, xprimv, xprimm025, xprimp025, rlatu1,
rlatu2, yprimu1, yprimu2 are computed at the same time as rlonu,
rlatu, rlonv, rlatv, and since it would not be convenient to separate
those computations, we decide to write xprimu, xprimv, xprimm025,
xprimp025, rlatu1, rlatu2, yprimu1, yprimu2 into "restart.nc", read
them from "start.nc" and not compute them in inigeom. So, in summary,
"start.nc" contains all the coordinates and their derivatives, and
inigeom only computes the 2D-variables.

Technical details:

Moved variables rlatu, rlonv, rlonu, rlatv, xprimu, xprimv from module
comgeom to module dynetat0_m. Upgraded local variables rlatu1,
yprimu1, rlatu2, yprimu2, xprimm025, xprimp025 of procedure inigeom to
variables of module dynetat0_m.

Removed unused local variable yprimu of procedure inigeom and
corresponding argument yyprimu of fyhyp.

Moved variables clat, clon, grossismx, grossismy, dzoomx, dzoomy,
taux, tauy from module serre to module dynetat0_m (since they are read
from "start.nc"). The default values are now defined in read_serre
instead of in the declarations. Changed name of module serre to
read_serre_m, no more module variable here.

The calls to fxhyp and fyhyp are moved from inigeom to etat0.

Side effects in programs other than gcm: etat0 and read_serre write
variables of module dynetat0; the programs test_fxyp and
test_inter_barxy need more source files.

Removed unused arguments len and nd of cv3_tracer. Removed unused
argument PPSOL of LWU.

Bug fix in test_inter_barxy: forgotten call to read_serre.
1 guez 3 MODULE start_init_orog_m
2    
3     ! From startvar.F, version 1.4
4     ! 2006/01/27 15:14:22 Fairhead
5    
6 guez 78 use dimens_m, only: iim, jjm
7    
8 guez 3 IMPLICIT NONE
9    
10 guez 78 REAL, SAVE:: mask(iim + 1, jjm + 1) ! interpolated fraction of land
11 guez 3
12 guez 78 private iim, jjm
13    
14 guez 3 CONTAINS
15    
16 guez 78 SUBROUTINE start_init_orog(phis, zmea_2d, zstd_2d, zsig_2d, zgam_2d, &
17     zthe_2d, zpic_2d, zval_2d)
18 guez 3
19     use conf_dat2d_m, only: conf_dat2d
20 guez 139 use dynetat0_m, only: rlatu, rlonv
21 guez 39 use grid_noro_m, only: grid_noro
22 guez 3 use indicesol, only: epsfra
23 guez 43 use netcdf, only: nf90_nowrite
24     use netcdf95, only: nf95_open, nf95_gw_var, nf95_inq_varid, nf95_close
25 guez 78 use nr_util, only: pi, assert
26 guez 3
27 guez 78 REAL, intent(out):: phis(:, :) ! (iim + 1, jjm + 1)
28     ! surface geopotential, in m2 s-2
29 guez 3
30 guez 78 REAL, intent(out):: zmea_2d(:, :) ! (iim + 1, jjm + 1) orographie moyenne
31    
32 guez 42 REAL, intent(out):: zstd_2d(:, :) ! (iim + 1, jjm + 1)
33 guez 3 ! (deviation standard de l'orographie sous-maille)
34    
35 guez 42 REAL, intent(out):: zsig_2d(:, :) ! (iim + 1, jjm + 1)
36 guez 3 ! (pente de l'orographie sous-maille)
37    
38 guez 42 REAL, intent(out):: zgam_2d(:, :) ! (iim + 1, jjm + 1)
39 guez 3 ! (anisotropie de l'orographie sous maille)
40    
41 guez 42 REAL, intent(out):: zthe_2d(:, :) ! (iim + 1, jjm + 1)
42 guez 3 ! (orientation de l'axe oriente dans la direction de plus grande
43     ! pente de l'orographie sous maille)
44    
45 guez 42 REAL, intent(out):: zpic_2d(:, :) ! (iim + 1, jjm + 1)
46     ! hauteur pics de la SSO
47 guez 3
48 guez 42 REAL, intent(out):: zval_2d(:, :) ! (iim + 1, jjm + 1)
49     ! hauteur vallees de la SSO
50    
51 guez 3 ! Local:
52    
53 guez 42 INTEGER iml_rel
54     INTEGER jml_rel
55 guez 43 INTEGER ncid, varid
56 guez 78 REAL, pointer:: relief(:, :)
57 guez 3 REAL, ALLOCATABLE:: lon_rad(:), lat_rad(:)
58 guez 43 REAL, pointer:: lon_ini(:), lat_ini(:)
59 guez 3
60     !-----------------------------------
61    
62     print *, "Call sequence information: start_init_orog"
63    
64 guez 78 call assert((/size(phis, 1), size(zmea_2d, 1), size(zstd_2d, 1), &
65     size(zsig_2d, 1), size(zgam_2d, 1), size(zthe_2d, 1), &
66     size(zpic_2d, 1), size(zval_2d, 1)/) == iim + 1, "start_init_orog iim")
67     call assert((/size(phis, 2), size(zmea_2d, 2), size(zstd_2d, 2), &
68     size(zsig_2d, 2), size(zgam_2d, 2), size(zthe_2d, 2), &
69     size(zpic_2d, 2), size(zval_2d, 2)/) == jjm + 1, "start_init_orog jjm")
70 guez 3
71 guez 43 print *, 'Reading the high resolution orography...'
72 guez 3
73 guez 43 call nf95_open('Relief.nc', nf90_nowrite, ncid)
74 guez 3
75 guez 43 call nf95_inq_varid(ncid, "longitude", varid)
76     call nf95_gw_var(ncid, varid, lon_ini)
77     lon_ini = lon_ini * pi / 180. ! convert to rad
78     iml_rel = size(lon_ini)
79 guez 3
80 guez 43 call nf95_inq_varid(ncid, "latitude", varid)
81     call nf95_gw_var(ncid, varid, lat_ini)
82     lat_ini = lat_ini * pi / 180. ! convert to rad
83     jml_rel = size(lat_ini)
84 guez 3
85 guez 43 call nf95_inq_varid(ncid, "RELIEF", varid)
86 guez 78 call nf95_gw_var(ncid, varid, relief)
87 guez 3
88 guez 43 call nf95_close(ncid)
89 guez 3
90 guez 43 ALLOCATE(lon_rad(iml_rel))
91 guez 3 ALLOCATE(lat_rad(jml_rel))
92    
93 guez 78 CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, relief , &
94 guez 3 interbar=.FALSE.)
95 guez 43 deallocate(lon_ini, lat_ini) ! pointers
96 guez 3
97     print *, 'Compute all the parameters needed for the gravity wave drag code'
98    
99 guez 78 CALL grid_noro(lon_rad, lat_rad, relief, rlonv, rlatu, phis, zmea_2d, &
100 guez 15 zstd_2d, zsig_2d, zgam_2d, zthe_2d, zpic_2d, zval_2d, mask)
101 guez 78 deallocate(relief) ! pointer
102 guez 3 phis(:, :) = phis(:, :) * 9.81
103    
104 guez 15 mask(2:, 1) = mask(1, 1) ! north pole
105     mask(2:, jjm + 1) = mask(1, jjm + 1) ! south pole
106     mask(iim + 1, 2:jjm) = mask(1, 2:jjm) ! Greenwich
107     WHERE (mask < EPSFRA)
108     mask = 0.
109     elsewhere (1. - mask < EPSFRA)
110     mask = 1.
111 guez 3 endwhere
112    
113     END SUBROUTINE start_init_orog
114    
115     END MODULE start_init_orog_m
  ViewVC Help
Powered by ViewVC 1.1.21