Parent Directory | Revision Log
Split "orografi.f": one file for each procedure. Put the created files in new directory "Orography". Removed argument "vcov" of procedure "sortvarc". Removed arguments "itau" and "time" of procedure "caldyn0". Removed arguments "itau", "time" and "vcov" of procedure "sortvarc0". Removed argument "time" of procedure "dynredem1". Removed NetCDF variable "temps" in files "start.nc" and "restart.nc", because its value is always 0. Removed argument "nq" of procedures "iniadvtrac" and "leapfrog". The number of "tracers read in "traceur.def" must now be equal to "nqmx", or "nqmx" must equal 4 if there is no file "traceur.def". Replaced variable "nq" by constant "nqmx" in "leapfrog". NetCDF variable for ozone field in "coefoz.nc" must now be called "tro3" instead of "r". Fixed bug in "zenang".
1 | program test_disvert |
2 | |
3 | ! This is a program in Fortran 95. |
4 | ! This program writes the vertical coordinates at level interfaces. |
5 | |
6 | use dimens_m, only: llm |
7 | use comvert, only: pa, disvert, ap, bp, preff |
8 | use new_unit_m, only: new_unit |
9 | |
10 | implicit none |
11 | |
12 | integer unit, l |
13 | |
14 | REAL s(llm+1) ! hybrid sigma-pressure coordinate |
15 | ! "s(l)" : à l'interface des couches "l" et "l-1" |
16 | |
17 | real p(llm + 1) ! pressure at level, in Pa |
18 | |
19 | !--------------------- |
20 | |
21 | print *, "llm = ", llm |
22 | pa = 5e4 |
23 | call disvert |
24 | s = ap / pa + bp |
25 | p = ap + bp * preff |
26 | |
27 | call new_unit(unit) |
28 | open(unit, file="test_disvert.csv", status="replace", action="write") |
29 | write(unit, fmt=*) '"ap (hPa)" "bp" "s" "p (hPa)"' ! title line |
30 | do l = 1, llm |
31 | write(unit, fmt=*) ap(l) / 100., bp(l), s(l), p(l) / 100. |
32 | end do |
33 | close(unit) |
34 | print *, 'The file "test_disvert.csv" has been created.' |
35 | |
36 | end program test_disvert |
ViewVC Help | |
Powered by ViewVC 1.1.21 |