/[lmdze]/trunk/dyn3d/test_fxhyp.f
ViewVC logotype

Annotation of /trunk/dyn3d/test_fxhyp.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 277 - (hide annotations)
Thu Jul 12 15:56:17 2018 UTC (5 years, 11 months ago) by guez
File size: 883 byte(s)
Move fxhyp and fyhyp to module dynetat0_m to avoid side effect on
variables of module dynetat0_m. A downside is that we need to link
heavyside, coefpoly and tanh_cautious into the gcm and test_fxhyp
executables.

We must move invert_zoom_x and principal_cshift to module dynetat0_m
to avoid circular dependency.

Move definition of rlatu(1) and rlatu(jjm + 1) inside fyhyp to avoid
side effect on rlatu.

1 guez 123 program test_fxhyp
2    
3 guez 277 use dynetat0_m, only: read_serre, fxhyp, xprimm025, rlonv, xprimv, rlonu, &
4     xprimu, xprimp025
5 guez 276 use unit_nml_m, only: unit_nml, set_unit_nml
6 guez 123
7     implicit none
8    
9     !--------------------------------------------------------
10    
11 guez 276 call set_unit_nml
12 guez 123 open(unit_nml, file="used_namelists.txt", status="replace", action="write")
13     call read_serre
14 guez 277 call fxhyp
15 guez 123 close(unit_nml)
16    
17 guez 276 ! We can use the same unit number although we are not writing a namelist:
18     open(unit_nml, file="test_fxhyp_out.txt", status="replace", action="write")
19     write(unit_nml, fmt = *) "xprimm025 = ", xprimm025
20     write(unit_nml, fmt = *) "rlonv = ", rlonv
21     write(unit_nml, fmt = *) "xprimv = ", xprimv
22     write(unit_nml, fmt = *) "rlonu = ", rlonu
23     write(unit_nml, fmt = *) "xprimu = ", xprimu
24     write(unit_nml, fmt = *) "xprimp025 = ", xprimp025
25     close(unit_nml)
26 guez 123
27     end program test_fxhyp

  ViewVC Help
Powered by ViewVC 1.1.21