3 |
michaesp |
1 |
Tropopause folding calculator:
|
|
|
2 |
(*original reference to be added*)
|
|
|
3 |
|
|
|
4 |
original work by B.Skerlak (bojan.skerlak@env.ethz.ch)
|
|
|
5 |
extended for netcdf data by A.Pozzer (andrea.pozzer@mpic.de)
|
|
|
6 |
|
|
|
7 |
------------------------------------------------------------------
|
|
|
8 |
To use the algorithm:
|
|
|
9 |
|
|
|
10 |
1) Extract the zip file
|
|
|
11 |
2) Modify the Makefile to the current machine set-up
|
|
|
12 |
3) Compile to code with the command "make".
|
|
|
13 |
Other options possible with "make help".
|
|
|
14 |
4) Create a namelist file (filename.nml) with the specific of the
|
|
|
15 |
netcdf files to be analyzed (see EMAC.nml for example).
|
|
|
16 |
5) run the code with the command "3d_labelling_and_fold_id.exe filename.nml"
|
|
|
17 |
|
|
|
18 |
-------------------------------------------------------------------
|
|
|
19 |
NML format:
|
|
|
20 |
the namelist file should contain the following fields:
|
|
|
21 |
------------------------------------------------------------------
|
|
|
22 |
&CTRL
|
|
|
23 |
! input file (path included)
|
|
|
24 |
file_input = '/ptmp/andrep/ZANIS/ZANIS__________20000801_0000_fold.nc'
|
|
|
25 |
! output file (path included)
|
|
|
26 |
file_output = 'test.nc'
|
|
|
27 |
! name of the longitude dimension in the file
|
|
|
28 |
X_name = 'lon'
|
|
|
29 |
! name of the latitude dimension in the file
|
|
|
30 |
Y_name = 'lat'
|
|
|
31 |
! name of the level dimension in the file
|
|
|
32 |
Z_name = 'lev'
|
|
|
33 |
! name of the time dimension in the file
|
|
|
34 |
T_name = 'time'
|
|
|
35 |
! name of the surface pressure
|
|
|
36 |
APS_name = 'aps'
|
|
|
37 |
! name of the hybrid coefficient A
|
|
|
38 |
HYAM_name = 'hyam'
|
|
|
39 |
! name of the hybrid coefficient B
|
|
|
40 |
HYBM_name = 'hybm'
|
|
|
41 |
! name of the specific humidity field
|
|
|
42 |
Q_name = 'qm1'
|
|
|
43 |
! name of the potential vorticity field
|
|
|
44 |
PV_name = 'PV'
|
|
|
45 |
! name of the potential temperature field
|
|
|
46 |
PT_name = 'tpot'
|
|
|
47 |
/
|
|
|
48 |
------------------------------------------------------------------
|
|
|
49 |
WARNINGS:
|
|
|
50 |
1) The code works automatically only if all the needed data (i.e. Q,PV,PT)
|
|
|
51 |
are included in the same input file.
|
|
|
52 |
2) Only hybrid pressure coordinates included in the pressure
|
|
|
53 |
calculation.
|
|
|
54 |
------------------------------------------------------------------
|