/[lmdze]/trunk/libf/test_disvert.f90
ViewVC logotype

Annotation of /trunk/libf/test_disvert.f90

Parent Directory Parent Directory | Revision Log Revision Log


Revision 10 - (hide annotations)
Fri Apr 18 14:45:53 2008 UTC (16 years, 2 months ago) by guez
File size: 798 byte(s)
Added NetCDF directory "/home/guez/include" in "g95.mk" and
"nag_tools.mk".

Added some "intent" attributes in "PVtheta", "advtrac", "caladvtrac",
"calfis", "diagedyn", "dissip", "vlspltqs", "aeropt", "ajsec",
"calltherm", "clmain", "cltrac", "cltracrn", "concvl", "conema3",
"conflx", "fisrtilp", "newmicro", "nuage", "diagcld1", "diagcld2",
"drag_noro", "lift_noro", "SUGWD", "physiq", "phytrac", "radlwsw", "thermcell".

Removed the case "ierr == 0" in "abort_gcm"; moved call to "histclo"
and messages for end of run from "abort_gcm" to "gcm"; replaced call
to "abort_gcm" in "leapfrog" by exit from outer loop.

In "calfis": removed argument "pp" and variable "unskap"; changed
"pksurcp" from scalar to rank 2; use "pressure_var"; rewrote
computation of "zplev", "zplay", "ztfi", "pcvgt" using "dyn_phy";
added computation of "pls".

Removed unused variable in "dynredem0".

In "exner_hyb": changed "dellta" from scalar to rank 1; replaced call
to "ssum" by call to "sum"; removed variables "xpn" and "xps";
replaced some loops by array expressions.

In "leapfrog": use "pressure_var"; deleted variables "p", "longcles".

Converted common blocks "YOECUMF", "YOEGWD" to modules.

Removed argument "pplay" in "cvltr", "diagetpq", "nflxtr".

Created module "raddimlw" from include file "raddimlw.h".

Corrected call to "new_unit" in "test_disvert".

1 guez 3 program test_disvert
2    
3     use dimens_m, only: llm
4     use comvert, only: pa, disvert, ap, bp, preff
5     use comconst, only: initialize
6 guez 7 use new_unit_m, only: new_unit
7 guez 3
8     implicit none
9    
10     integer unit, l
11     REAL sigma(llm+1) ! "sigma(l)" : à l'interface des couches "l" et "l-1"
12     real p(llm + 1) ! pressure at level, in Pa
13    
14     !---------------------
15    
16     print *, "llm = ", llm
17     call initialize
18     pa = 5e4
19     call disvert
20     sigma = ap / pa + bp
21     p = ap + bp * preff
22    
23 guez 10 call new_unit(unit)
24 guez 3 open(unit, file="test_disvert.csv", status="replace", action="write")
25     write(unit, fmt=*) '"ap (Pa)" "bp" "sigma" "p (Pa)"' ! title line
26     do l = 1, llm
27     write(unit, fmt=*) ap(l), bp(l), sigma(l), p(l)
28     end do
29     close(unit)
30     print *, 'The file "test_disvert.csv" has been created.'
31    
32     end program test_disvert

  ViewVC Help
Powered by ViewVC 1.1.21