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

Diff of /trunk/dyn3d/test_fxhyp.f

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

trunk/Sources/test_fxhyp.f revision 139 by guez, Tue May 26 17:46:03 2015 UTC trunk/dyn3d/test_fxhyp.f revision 277 by guez, Thu Jul 12 15:56:17 2018 UTC
# Line 1  Line 1 
1  program test_fxhyp  program test_fxhyp
2    
3    USE dimens_m, ONLY: iim    use dynetat0_m, only: read_serre, fxhyp, xprimm025, rlonv, xprimv, rlonu, &
4    use fxhyp_m, only: fxhyp         xprimu, xprimp025
5    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  
6    
7    implicit none    implicit none
8    
   REAL, dimension(iim + 1):: xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025  
   integer unit  
   
9    !--------------------------------------------------------    !--------------------------------------------------------
10    
11    call new_unit(unit_nml)    call set_unit_nml
12    open(unit_nml, file="used_namelists.txt", status="replace", action="write")    open(unit_nml, file="used_namelists.txt", status="replace", action="write")
13    call read_serre    call read_serre
14    call fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025)    call fxhyp
15    close(unit_nml)    close(unit_nml)
16    
17    unit = unit_nml    ! We can use the same unit number although we are not writing a namelist:
18    open(unit, file="test_fxhyp_out.txt", status="replace", action="write")    open(unit_nml, file="test_fxhyp_out.txt", status="replace", action="write")
19    write(unit, fmt = *) "xprimm025 = ", xprimm025    write(unit_nml, fmt = *) "xprimm025 = ", xprimm025
20    write(unit, fmt = *) "rlonv = ", rlonv    write(unit_nml, fmt = *) "rlonv = ", rlonv
21    write(unit, fmt = *) "xprimv = ", xprimv    write(unit_nml, fmt = *) "xprimv = ", xprimv
22    write(unit, fmt = *) "rlonu = ", rlonu    write(unit_nml, fmt = *) "rlonu = ", rlonu
23    write(unit, fmt = *) "xprimu = ", xprimu    write(unit_nml, fmt = *) "xprimu = ", xprimu
24    write(unit, fmt = *) "xprimp025 = ", xprimp025    write(unit_nml, fmt = *) "xprimp025 = ", xprimp025
25    close(unit)    close(unit_nml)
26    
27  end program test_fxhyp  end program test_fxhyp
Legend:
Removed from v.139  
changed lines
  Added in v.277
  ViewVC Help
Powered by ViewVC 1.1.21