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

Diff of /trunk/dyn3d/test_fxhyp.f90

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

revision 265 by guez, Tue Mar 20 09:35:59 2018 UTC revision 276 by guez, Thu Jul 12 14:49:20 2018 UTC
# Line 1  Line 1 
1  program test_fxhyp  program test_fxhyp
2    
3    USE dimensions, ONLY: iim    USE dimensions, ONLY: iim
4      use dynetat0_m, only: read_serre
5    use fxhyp_m, only: fxhyp    use fxhyp_m, only: fxhyp
6    use jumble, only: new_unit    use unit_nml_m, only: unit_nml, set_unit_nml
   use read_serre_m, only: read_serre  
   use unit_nml_m, only: unit_nml  
7    
8    implicit none    implicit none
9    
10    REAL, dimension(iim + 1):: xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025    REAL, dimension(iim + 1):: xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025
   integer unit  
11    
12    !--------------------------------------------------------    !--------------------------------------------------------
13    
14    call new_unit(unit_nml)    call set_unit_nml
15    open(unit_nml, file="used_namelists.txt", status="replace", action="write")    open(unit_nml, file="used_namelists.txt", status="replace", action="write")
16    call read_serre    call read_serre
17    call fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025)    call fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025)
18    close(unit_nml)    close(unit_nml)
19    
20    unit = unit_nml    ! We can use the same unit number although we are not writing a namelist:
21    open(unit, file="test_fxhyp_out.txt", status="replace", action="write")    open(unit_nml, file="test_fxhyp_out.txt", status="replace", action="write")
22    write(unit, fmt = *) "xprimm025 = ", xprimm025    write(unit_nml, fmt = *) "xprimm025 = ", xprimm025
23    write(unit, fmt = *) "rlonv = ", rlonv    write(unit_nml, fmt = *) "rlonv = ", rlonv
24    write(unit, fmt = *) "xprimv = ", xprimv    write(unit_nml, fmt = *) "xprimv = ", xprimv
25    write(unit, fmt = *) "rlonu = ", rlonu    write(unit_nml, fmt = *) "rlonu = ", rlonu
26    write(unit, fmt = *) "xprimu = ", xprimu    write(unit_nml, fmt = *) "xprimu = ", xprimu
27    write(unit, fmt = *) "xprimp025 = ", xprimp025    write(unit_nml, fmt = *) "xprimp025 = ", xprimp025
28    close(unit)    close(unit_nml)
29    
30  end program test_fxhyp  end program test_fxhyp
Legend:
Removed from v.265  
changed lines
  Added in v.276
  ViewVC Help
Powered by ViewVC 1.1.21