| 2 |
|
|
| 3 |
implicit none |
implicit none |
| 4 |
|
|
| 5 |
REAL:: clon= 0. ! longitude en degres du centre du zoom |
REAL:: clon = 0. ! longitude of the center of the zoom, in rad |
| 6 |
|
real:: clat = 0. ! latitude of the center of the zoom, in rad |
| 7 |
|
|
| 8 |
real:: clat= 0. ! latitude en degres du centre du zoom |
real:: grossismx = 1., grossismy = 1. |
| 9 |
|
! facteurs de grossissement du zoom, selon la longitude et la latitude |
| 10 |
|
! = 2 si 2 fois, = 3 si 3 fois, etc. |
| 11 |
|
|
| 12 |
real, save:: transx,transy |
real:: dzoomx = 0.2, dzoomy = 0.2 |
| 13 |
real, save:: alphax, alphay ! anciennes formulations des grossissements |
! extensions en longitude et latitude de la zone du zoom (fractions |
| 14 |
real, save:: pxo,pyo |
! de la zone totale) |
| 15 |
|
|
| 16 |
real:: grossismx= 1.0 ! facteur de grossissement du zoom, selon la longitude |
real:: taux = 3., tauy = 3. |
| 17 |
|
! raideur de la transition de l'intérieur à l'extérieur du zoom |
| 18 |
|
|
| 19 |
|
contains |
| 20 |
|
|
| 21 |
real:: grossismy= 1.0 ! facteur de grossissement du zoom, selon la latitude |
subroutine read_serre |
| 22 |
|
|
| 23 |
real:: dzoomx= 0.0 |
use unit_nml_m, only: unit_nml |
| 24 |
! (extension en longitude de la zone du zoom (fraction de la zone totale)) |
use nr_util, only: assert, pi |
| 25 |
|
|
| 26 |
real:: dzoomy= 0.0 |
REAL:: clon_deg = 0. ! longitude of the center of the zoom, in degrees |
| 27 |
! extension en latitude de la zone du zoom |
real:: clat_deg = 0. ! latitude of the center of the zoom, in degrees |
|
! (fraction de la zone totale) |
|
| 28 |
|
|
| 29 |
real:: taux= 3.0 ! raideur du zoom en X |
namelist /serre_nml/ clon_deg, clat_deg, grossismx, grossismy, dzoomx, & |
| 30 |
|
dzoomy, taux, tauy |
| 31 |
|
|
| 32 |
real:: tauy= 3.0 ! raideur du zoom en Y |
!------------------------------------------------- |
| 33 |
|
|
| 34 |
|
print *, "Enter namelist 'serre_nml'." |
| 35 |
|
read(unit=*, nml=serre_nml) |
| 36 |
|
write(unit_nml, nml=serre_nml) |
| 37 |
|
|
| 38 |
|
call assert(grossismx >= 1. .and. grossismy >= 1., "read_serre grossism") |
| 39 |
|
call assert(dzoomx > 0., dzoomx < 1., dzoomy < 1., & |
| 40 |
|
"read_serre dzoomx dzoomy") |
| 41 |
|
clon = clon_deg / 180. * pi |
| 42 |
|
clat = clat_deg / 180. * pi |
| 43 |
|
|
| 44 |
|
end subroutine read_serre |
| 45 |
|
|
| 46 |
end module serre |
end module serre |
Parent Directory
| 
Revision Log
| 
Patch