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