Home

Planet Simulator User's Guide

1. DSET puma m gra UNDEF 9e 09 XDEF 64 LINEAR 0 0000 5 6250 OPTIONS YREV YDEF 32 LEVELS 85 7606 80 2688 74 7445 69 2130 63 6786 58 1430 52 6065 41 5325 35 9951 30 4576 24 9199 19 3822 13 8445 8 3067 2 7689 8 3067 13 8445 19 3822 24 9199 30 4576 35 9951 47 0696 52 6065 58 1430 63 6786 69 2130 74 7445 80 2688 ZDEF 1 LINEAR 1 1 TDEF 12 LINEAR 00 00Z01jan0001 1mo VARS 3 c139 0 99 139 0 0 c151 0 99 151 0 0 c175 0 99 175 0 0 ENDVARS Here the line starting with TDEF ends with 1mo since we are handling monthly mean data When the PUMA output is used without averaging this should corre spond to the output interval given by the nafter variable used in the namelist of your PUMA run see section The number of variables depends on how the pumaburner was called In this example only 3 variables were processed i e the surface temperature c139 the sea level pressure c151 and the albedo c175 refer to appendix for a list of codes The GrADS program is started by typing grads in a terminal window Then data is visualised either by typing commands line by line or preferably by using scripts The following script called tglob gs displays the monthly mean surface temperature tglob gs function pass m reinit open puma_m enable print print mf set t m 47 0696 2 7689 41 5325 85 7606 6 1 GRADS 67 o set gxout shaded d c139 273 16 cbar gs
2. set gxout contour d c139 273 16 draw title Surface Temperature deg C month m print disable print Igxps i print mf o tglob m ps The variable m at the beginning of the script defines the month which should be displayed It is passed from the terminal with the script call Note that in this line no quotation marks are present since only GrADS specific commands are framed by quotation marks Script commands like variable definitions if clauses etc are used without quotation marks The script is executed by typing its name without the ending and the number of the month to be shown For example tglob 7 displays the monthly mean surface temperature in July The resulting output file is called tglob7 ps The following script thh displays the time dependent surface temperature of Hamburg Here two variables are passed to GrADS the first and last day to plot note that here the file puma gra is opened which contains data on a daily basis The call thh 91 180 displays the surface temperature of Hamburg for the spring season from April 1st to June 30th thh gs function pass d1 d2 reinit open puma enable print print mf set lat 53 set lon 10 set t d1 d2 0 d c139 273 16 draw title Surface Temperature deg C in Hamburg print disable print Igxps i print mf o thh ps It is possible to have more than one figur
3. 30 3 2 BATCH MODE 37 OUTYEAR INYEAR OUTFTP expr INFTP 1 if OUTMON eq 13 then OUTMON 1 OUTYEAR expr OUTYEAR 1 fi echo OUTYEAR gt EXPDIR saveyear EXP echo OUTMON gt EXPDIR savemon EXP echo OUTDAY gt EXPDIR saveday EXP HHH HHH HHH HHH HHH HHH HHHH HHH HHH if OUTFTP gt FTPINT then cd EXPDIR OUTFTP 1 echo OUTFTP gt EXPDIR saveftp EXP tar tar tar tar tar tar tar tar tar tar tar tar tar rf rf rf rf rf rf rf rf rf rf rf rf rf EXPDIR TARFILE EXPDIR TARFILE EXPDIR TARFILE EXPDIR TARFILE EXPDIR TARFILE EXPDIR TARFILE EXPDIR TARFILE EXPDIR TARFILE EXPDIR TARFILE EXPDIR TARFILE EXPDIR TARFILE EXPDIR TARFILE EXPDIR TARFILE EXP RES EXP LANDRES EXP SEARES NAMLIST exe land_namelist sea_namelist ice_namelist ocean_namelist saveday EXP savemon EXP saveyear EXP savetar EXP saveftp EXP tar rf EXPDIR TARFILE model x tar rf EXPDIR TARFILE runall mv EXPDIR TARFILE DATADIR TARFILE to YY MM if SCHAUER eq 1 then cat gt EXPDIR put_schauer YY MM lt lt EOP set ex cd EXPDIR home larry utils uput2 DATADIR TARFILE to YY MM SCHAUERD
4. destination path of the installation The default paths and the names of the corresponding enviroment variables for the directories are e Java2 binary usr lib java bin JAVA BIN e Perl binary usr bin PERL_BIN e GrADS binary usr local grads bin GRADS_BIN e target directory PlaisirGUI ttp www cpan org http www activestate com http www blackdown org http java sun com http grads iges org grads 41 To accept a default path the input can be confirmed with lt RETURN gt For the directories also relative paths are allowed When you already entered all paths you can correct individual directories or confirm the installation To use GrADS in addition to the binary path two more environment variables have to be set If these variables are not found by the installation script the path names will be queried These two variables are listed in the following the given path names are the standard paths which are used also in the installation script e GrADS font files usr local grads dat GADDIR e GrADS library usr local grads lib GASCRP After completion of the installation script you can change the working directory to the PlaisirGUI directory and the GUI can be started cd PlaisirGUI bin plaisirgui Alternativly you can put the binary directory to the standard path variable PATH E g bash export PATH lt gui_home_dir gt bin PATH Afterwards you can start the GUI program
5. option c option d option g option n option m InputFile OutputFile namelist printout help this output print available codes and names debug mode verbose output Grib output override namelist option NetCDF output override namelist option Mean 1 output override namelist option Planet Simulator or PUMA data file GRIB SERVICE or NetCDF format file redirected lt stdin gt redirected lt stdout gt 5 1 3 Namelist The namelist values control the selection coordinate system and output format of the postprocessed variables Names and values are not case sensitive You can assign values to the following names Name Def Type Description Example HTYPE S char Horizontal type HTYPE G VTYPE S char Vertical type VTYPE P MODLEV 0 int Model levels MODLEV 2 3 4 hPa 0 real Pressure levels hPa 500 1000 CODE 0 int ECMWF field code CODE 130 152 GRIB 0 int GRIB output selector GRIB 1 NETCDF 0 int NetCDF output selector NETCDF 1 MEAN 1 int Compute monthly means MEAN 0 HHMM 1 int Time format in Service format HHMM 0 HEAD7 0 int User parameter HEAD7 0815 MARS O int Use constants for planet Mars MARS 1 MULTI 0 int Process multiple input files MULTI 12 5 1 3 1 HTYPE HTYPE accepts the first character of the following string Following settings are equivalent HTYPE S HTYPE Spherical Harmonics HTYPE Something Blanks and the equal sign are
6. 4 0 4 Troubleshooting and the Programmers Guide part II are more details about logging PlaisirGUl Intel Set model properties Select interactive parameters Set program parameters at Log Window T Advanced Logging Clear Log Files Figure 4 8 Entries in Logging Menu 3 Wind Force Mean Temneratire 4 0 3 4 Running the Simulation Program Figure shows the control panel together with three GrADS frames which display the actual state of temperature and wind force The control panel is divided into three parts Interactive Visualisation and Simulation With the displayed sliders in the first part of the main window model parame ters can be changed interactively during the simulation In order to pass these adjustments on to the simulation the running program must be interrupted and the button Set of the section Interactive has to be pushed When the simulation program has accepted the new parameters the calculations can be continued The button Reset puts all parameters to the initial values After pressing the button Set the old setting is lost 8 og4j http jakarta apache org log4j File Options X PlaisirGUI ol HA Graps 1 7Betag Temperature in Grad Celsius Date 0001 04 27 o Interactive co2 co2 concentration ppm 360 0 324 0 396 0 dz0land roughness length land Ea 2 0 1 8 2 2 albice max albedo for sea ice Ko 0 7 0 63 0
7. 77 RESET Visualisation vi Temperature Edit Visu Parameters v Wind Force Close all visu windows vi Temperature Open all visu windows C Diff Of Mean Temperature um Running step 3725 of 1000 300 100 Wind Force in m ate 0001 04 27 HA GrADS 1 7Beta9 lt 3 gt Temperature in Grad Celsius Dote 0001 04 27 WE 15E 2DE 25 Figure 4 9 Control Panel 49 50 CHAPTER 4 GRAPHICAL USER INTERFACE The second part of the window shows the check boxes for the interactive visu alisation of simulation results Up to four variables can be selected before and during the simulation The graphical display of simulation results starts as soon as the simulation is running There are three buttons on the right side of the visulisation panel The first button opens a dialogbox to edit the visualisation parameters which is described below The middle button opens all GrADS windows and the last button closes all open Gr ADS applications The third sector serves for the control of the simulation where the symbols stand for the following gt starting continuing the simulation u interrupting the simulation E terminating the simulation The START button starts the simulation If variables were already activated for displaying results before the start they are shown now in separate windows To interrupt the si
8. Temperature lon lat lev t Temperature lat lev u Zonal wind component lat lev du Zonal wind component lon lat lev v Meridional wind component lat lev dv Meridional wind component lon lat lev Dialogbox for the Common properties To configure the GrADS output the GUI includes a dialogbox which can be opened with the button Edit Visu Parameter in the Visualisation panel Figure shows this dialog box longitude latitude Meross section along the Greenwich meridian 12atmospheric level 52 CHAPTER 4 GRAPHICAL USER INTERFACE D ox Visu Parameter Title Temperature in Kelvin Europe Du Short Title Temp Europe x Formula dt v 7 Display mean of values vi Display monthly deviation of values Level 1000 I Longitude 15 30 a Latitude 90 90 2 Color Level 250 255 260 265 270 275 280 285 290 295 300 bi Cancel OK Figure 4 11 Dialog Box for the Detailed Properties The dialogbox panel shows four rows with buttons and checkboxes etc and finally a Cancel and an OK button The first item in the row is a checkbox to enable the editing of the corresponding GrADS window properties If this checkbox is not selected it is not possible to edit this item in the dialogbox and you can t start the associated GrADS windows in the visualisation panel of the main window The second item in the row is a pop down list The list title is a short description of
9. data have to be located there e g surface_parameter and all output data will be put there The first time the GUI is started the working directory is set to the starting directory when the program is restarted the working directory is set to the last used path With the menu item File Select simulation executable the simulation program can be adjusted for the next run It will be inquired whether the model con figuration file is to be used whose name is derived from the program i e pro gram cfg Alternatively the pre selected configuration will be kept in use Loading and Saving model configuration Under File Load model configuration file a model configuration is specified which contains the initialisation of all modelling parameters see part II sec 44 CHAPTER 4 GRAPHICAL USER INTERFACE tion It is necessary to configure the model before choosing or modifying the static and interactive modelling parameters see section 4 0 3 3 File Options Select working directory Select simulation executable Load model configuration file le Save model configuration file Save model configuration file as Exit program fame _ Temperature Close all visu wil VA Draccnra EO X___ _ AA Figure 4 2 File Menu The next file menu items save the current model The simple File Save model configuration file item overwrites the old configuration In the item File Save model conf
10. for use on scalar machines instead of mpimod f90 1 2 UNPACKING AND INSTALLATION 9 outmod 90 pumamod f90 puma f90 radmod f90 rainmod f90 seamod_climatology f90 surfmod f90 Data file surface_parameter Namelists puma_namelist land_namelist sea_namelist Scripts and Makefile output module for global variables to be used in other modules atmospheric main programm and dynamic radiation long and short wave moist processes large scale and convective precipitation cloud cover climatological sea surface ocean and sea ice interface to surface modules sea and land description initial fields and climatologies orography land sea mask glacier mask surface temperature sea ice cover AMIP SSTs are used for this climatology description example namelist for puma typical set up example namelist for landmod f90 example namelist for seamod_climatology f90 description 10 Makefile compile_mpp runscript_template Additional pumaburn4 c srv2gra f90 srv2ascii f90 PumaGui install tar gz PumaGui doc pdf README 0 2 1 CHAPTER 1 INSTALLATION makefile to compile PUMA scalar machine script to compile PUMA for MPI see NPRO in PUMAMOD f90 example script to run PUMA for some years unix description pumaburner with netcdf export support compile with standard c compiler e g cc 02 o pumaburn4 D_FILE_OFFSET_BITS 64 pumaburn
11. logfiles 4 0 5 Additional Tools For demonstration purposes data snapshots of a special PUMA run can be saved to show them later as a movie which the GUI can display The look and feel of the visualized data is the same as that of a concrete PUMA run The snapshots are also the basis to create a reference file used to compute and visualise deviations to the actually computed data by PUMA Depending on the value of writeStep default 32 set in the GUI namelist see part II section an output of the actual values of some PUMA variables will be written into the snapshot netCDF file puma_var nc Additionally after each month the mean values over all written values for the preceeding month are written into the netCDF file puma_mean nc These files are read permanently by the visualisation part of the GUI 4 0 5 1 Running a Movie If the parameter makemovie default false set in the GUI namelist see part II section is set to true a copy of each newly written file puma_var nc and puma mean nc is stored in the directory pumamovie which will be cre ated if necessary These copies will have the filenames puma_var nc_000000 and puma_mean nc_000000 where 000000 is replaced by the actual timestep The names of the files are saved in the files puma_var dat and puma_mean dat for use in the puma_movie tool For using this tool it is necessary to set the simulation executable File Se lect simulation executable to the examples exe
12. m I mixing length for momentum m Ls latent heat of sublimation 2 8345 10 Jkg Lsn latent heat of fusion for snow 3 32 10 TEBA Ly latent heat of vapourization 2 5008 10 Jkg e precipitation ms P u associated Legendre function of the first kind p pressure Pa Ps surface pressure Pa Ps scaled surface pressure q specific humidity kg ke Q total heat flux through sea ice Wm Q flux correction heat flux through sea ice Wm Qa total atmospheric heat flux Wm Qe conductive heat flux through sea ice Wm Qf heat flux available for freezing sea ice Wm Qo oceanic heat flux Wm ds surface specific humidity kg ke dsa saturation specific humidity kg ke Ra cloud albedo Ra gas constant for dry air 287 05 Jee PK Ri surface long wave radiation Wm Rs surface short wave radiation W m R gas constant for water vapor 461 51 Jkg7t KU Ro zeroth moment of the temperature distribution Km R first moment of the temperature distribution K m Ri Richardson number Sw salinity of sea water 34 7 psu 74 APPENDIX A LIST OF CONSTANTS AND SYMBOLS Symbol Definition Value Unit t time S t scaled time step T temperature K T temperature anomaly T To _ Ta deep ocean temperature at 400m K Te sea ice surface temperature K Ty freezing temperature 271 25 K Es surface temperature K Tmel melting point 273 16 K Paiz mixed layer temperature K Tmiz climatological mixed layer temperature K Tref asymptotic referenc
13. of hPa are ignored hPa is a real array that accepts pressure values with the units hectoPascal or millibar All output variables will be inter polated to the selected pressure levels There is no extrapolation on the top of the atmosphere For pressure values that are lower than that of the model s top level the top level value of the variable is taken The variables temperature and geopotential height are extrapolated if the selected pressure is higher than the surface pressure All other variables are set to the value of the lowest mode level for this case The outputfile contains the levels in the same order as set in hPa Example hpa 100 300 500 700 850 900 1000 60 CHAPTER 5 POSTPROCESSING 5 1 3 5 MEAN MEAN can be used to compute montly means and or deviations The Pum aburner reads date and time information from the model file and handles different lengths of months and output intervals correctly Setting Description MEAN 0 Do no averaging all terms are processed MEAN 1 Compute and write monthly mean fields Not for spherical har monics Fourier coefficients or zonal means on sigma levels MEAN 2 Compute and write monthly deviations Not for spherical harmon ics Fourier coefficients or zonal means on sigma levels Deviations are not available for NetCDF output MEAN 3 Combination of MEAN 1 and MEAN 2 Each mean field is fol lowed by a deviation field with an identical header record Not fo
14. to convert it to either GRIB for instance with the PINGOs grb copy2 data srv data_with_grib_metainfo grb output grb or NETCDF using the program puma2cdf which is available with the PUMA postprocessing tools Despite of its name this program cannot process raw PUMA output but takes SERVICE format as input It can as well be called as srv2cdf which changes its behaviour oddities of model output such as the existence of February 30 are then no longer removed Once the format is changed proceed from there e your data is in NETCDF format it can easily transformed to Vis5D s native format by means of the program cdf2v2d which is available with the PUMA postprocessing tools e your data is in GRIB format you can find a transformation tool named Grib2V5d at lt http grib2v5d sourceforge net gt which offers various practical fea tures Once the conversion to Vis5D s native format is achieved please follow the in structions from the Vis5D documentation or if Vis5D is already installed on your system try finding your own way by typing visbd my_data v5d 71 72 APPENDIX A LIST OF CONSTANTS AND SYMBOLS Appendix A List of Constants and Symbols A 0 1 List of Constants and Symbols SymbolDefinition Value Unit a earth radius 6371 10 m A surface albedo A D V Vinp Aq cloud absorptivity _ char Charnock constant 0 018 Ch transfer coefficient for heat ce drag coefficient for momentum c s
15. 1 g surface albedo frac 176 1 ga surface solar radiation W m 177 1 ga surface thermal radiation W m 178 1 ga top solar radiation W m 179 1 ga top thermal radiation W m 180 1 ga u stress Pa 181 1 ga v stress Pa 182 1 ga evaporation m s 183 1 g soil temperature K 203 1 ga top solar radiation upward W m 204 1 ga surface solar radiation upward W m 205 1 ga surface thermal radiation upward W m 207 1 g soil temperature level 2 K 208 1 g soil temperature level 3 K 209 1 g soil temperature level 4 K 210 1 g sea ice cover frac 211 1 g sea ice thickness m 212 1 g vegetation cover frac 218 1 g snow melt water equiv m s 221 1 g snow depth change water equiv m s 230 1 ga vertical integrated spec hum kg m 232 1 g glacier cover frac 79 s PUMA spectral field g PUMA grid point field c computed by PUMA burner a accumulated
16. 4 c lm lnetcdf convert service format typical afterburner output to grads input control and data compile it with FORTRAN9O compiler usage srv2gra INPUT will create an INPUT ctl and an INPUT gra file Open INPUT ctl using GRADS on some machines you have to modify the record length converts service format model output to its ascii representation presently using the fortran format 110 for the header and the format 8E12 6 for the data graphical User Interface for PLASIM see Chapter 4 documentation for the graphical user interface README for the graphical user interface 1 3 Compiling and running PUMA To compile PUMA you may take these steps e use make for compiling a set of modules Just type make in the module directory and the executable puma x will be created e The utilities pumaburn and srv2gra need to be compiled separately The package pumaburn2 tar gz contains a Makefile which might have to be 1 3 COMPILING AND RUNNING PUMA 11 customized before invocation The program srv2gra f90 only needs to be compiled with any FORTRAN90 compiler e g f90 o srv2gra srv2gra f90 After compiling create a directory e g run and copy the executable puma x all initial fields _parameter and the namelists _namelist into it Addtion ally you can run PUMA on maschines with more than one processor using the compile_mpp script see file list The latter
17. APPENDIX B PUMA CODES Codes available from PUMA burner adapted from ECHAM Code Levels Type Variable Unit 110 1 g mixed layer depth m 129 1 S surface geopotential m s 130 NLEV s temperature K 131 NLEV c u velocity m s 132 NLEV e v velocity m s 133 NLEV s specific humidity kg kg 135 NLEV c vertical velocity Pa s 138 NLEV s vorticity 1 s 139 1 g surface temperature K 140 1 g soil wetness m 141 1 g snow depth water equi m 142 1 ga large scale precipitation m s 143 1 ga convective precipitation m s 144 1 ga snow fall m s 146 1 ga surface sensible heat flux W m 147 1 ga surface latent heat flux W m 148 NLEV e horizontal streamfunktion m s 149 NLEV c velocity potential m s 151 1 c mean sea level pressure Pa 152 1 S In surface pressure 153 NLEV g cloud liquid water content kg kg 155 NLEV s divergence 1 s 156 NLEV c geopotential height gpm 157 NLEV c relative humidity frac 159 1 g u m s 160 1 ga surface runoff m s Code Levels Type Variable Unit 162 NLEV g cloud cover frac 164 1 ga total cloud cover frac 169 1 ga surface temperature K 170 1 g deep soil temperature K 172 1 g land sea mask O sea 1 land 173 1 g surface roughness m 175
18. D ar e a a e en 59 SL IAS a o a ss ts ns re re 60 ETA MULTI 2 0 2 04 4 00 e A A ee 60 og Serge Y oe ce AN LR ee 60 5 1 11 Troubleshooting 246244 2 4 ae ee ewe eRe 61 63 A ee WA ane aa bo Gee ee ee ee ee 63 Deh Ae A ee oe eS WAUA ee oe ee ee 67 68 CONTENTS A List of Constants and Symbols A 0 1 List of Constants and Symbols B_ Puma Codes CONTENTS Chapter 1 Installation 1 1 Quick start This section is for the very unpatient used to skip manuals and just eager to see results as fast as possible The quick start just guides you till the first data output is produced The next steps the post processing and the graphical presentation are described in chapter 5 6 For a more detailed description of what this package is all about read the paragraph following this quick start type mkdir plasim e type cp plasim tar plasim e type cd plasim e type tar xvf plasim tar e type mkdir run e type cd module e type make e type cp puma x run e type cd e type cp puma x run e type cp data _namelist data _parameter run e type cd run 8 CHAPTER 1 INSTALLATION e type puma x You may need to adjust the Makefile to your environment in particular the fortran compiler command and the c compiler command The last step may take a while to be finished depending on your maschine s processing speed After all you should find four output files
19. DIR model x This first block of the script defines the basic settings i e the directories used the length of the integration etc If SCHAUER is set to 1 all puma output is trans ferred to the specified directory on the schauer thus avoiding the users directory 3 2 BATCH MODE 33 from filling up Otherwise the output is written to the DATADIR directory A temporary directory is created where the model is eventually run cat gt EXPDIR NAMLIST exe lt lt EOX NAME LIST PARAMETER cat gt puma_namelist lt lt EOF amp INP NDAYS 30 NTSPD 32 NRESTART 1 NDIAG 480 NAFTER 32 NEQSIG 1 PSURF 101325 NPACKSP 0 NPACKGP 0 amp END amp MISCPAR amp END amp FLUXPAR amp END amp RADPAR amp END amp RAINPAR amp END amp SURFPAR amp END EOF EOX cat gt EXPDIR land_namelist lt lt EOL amp landpar amp end EOL cat gt EXPDIR sea_namelist lt lt EOO amp seapar amp end E00 cat gt EXPDIR ocean_namelist lt lt E00 34 CHAPTER 3 RUNNING PUMA amp oceanpar amp end E00 cat gt EXPDIR ice_namelist lt lt E00 amp icepar amp end E00 chmod u x EXPDIR NAMLIST exe Now the necessary namelists are generated The puma namelist is defined as an executable which can be called with a parameter setting the restart mode EXPDIR NAMLIST exe 0 cp EXPDIR land_namelist land_namelist cp EXPDIR sea_namelist
20. IR TAF 38 CHAPTER 3 RUNNING PUMA rm DATADIR TARFILE to YY MM EOP chmod u x EXPDIR put_schauer YY MM EXPDIR put_schauer YY MM gt EXPDIR put_ YY MM 2 gt amp 1 amp fi cd TMPDIR fi echo OUTFTP gt EXPDIR saveftp EXP if OUTYEAR gt LASTYEAR then cd EXPDIR rm r TMPDIR exit fi exit EOR cd EXPDIR chmod u x runall runall etime date echo echo gt started stime echo echo gt done etime echo exit Chapter 4 Graphical User Interface 4 0 1 Introduction The Graphical User Interface GUI for the Planet Simulator provides the user with an interactive mode designed to aid in tuning of parameterization and de bugging of the PUMA simulation program The configuration of different experiments is handled easily by the structured graphical input of all model parameters predefined for unexperienced users those who remain static during a run as well as the user defined interactive parameters who are selected and changed at any time with direct effect on the running model An online visualisation of climate variables as they are computed may be dis played during the run the visualisation parameters may also be adapted during the run The GUI supports the control of different versions of the simulation program with different configurations which can run also remotely on ano
21. Planet Simulator User s Guide F Lunkeit K Fraedrich H Jansen E Kirk U Luksch F Sielmann K Cassirer W Joppich M Lob M Preuhs December 17 2003 Contents 5 MEA EN 5 Sore at Seon ee Oe Gi acc ey wa ana ee he Braten 6 SS pags ec Scat hacks ee ok KUNI gg a E 8 131 pumaburner 9 bom ede e a EA TEN 10 AAA A 10 2 Modules 11 o TO ee al eo el a s 13 22 MISCO AA e AA ri 14 28 SUMA e a AA ira RO ES 15 A A TER 16 A A 17 20 SUNNA LIA s s nob one a ce eg 18 2 7 Seaiceandoceanmodules 19 A AA hae Oe RE EERE LEE RPE 23 2 7 2 mtermodatm f90 24 2 7 3 mtermodice f90 25 2TA REMO OO re a aa eR a 26 NO A 27 276 oceanmod50 f90 28 4 CONTENTS 29 3 1 Interactive Console Model o 29 3 2 Batch Mode y e spo i ras a h 30 37 40 1 Introductioni 37 4 0 2 Installing the PlaisirGUI 37 4 0 3 Operatingi 40 4 0 4 Trouble Shootingi 51 4 0 5 Additional Tools 52 5 Postprocessing 55 5 17 Pissbumenl sese 2 40 44 mia seia eee oe AA 55 5 1 1 Inmtroduction 55 KUWA Eco ok da Sade Le Me enge ot nd ee ee 56 5 13 Names ci e e WAA a Aa 24 4 amp 56 FONTES 58 0 1 5 SERVICE format AA 59 ee A PA 59 Aa HEA
22. ad 8 experiment number extracted from filename Example for reading the SERVICE format GRIB 0 NETCDF 0 INTEGER HEAD 8 REAL FIELD 64 32 dimensions for T21 grids READ 10 ERR 888 END 999 HEAD READ 10 ERR 888 END 999 FIELD 888 STOP I 0 ERR 999 STOP EOF 5 1 6 HHMM Setting Description HHMM 0 head 4 shows the time in hours HH HHMM 1 head 4 shows the time in hours and minutes HHMM 5 1 7 HEAD7 The 7th element of the header is reserved for the user It may be used for experiment numbers flags or anything else Setting HEAD7 to a number exports this number to every header record in the output file SERVICE format only 62 CHAPTER 5 POSTPROCESSING 5 1 8 MARS This parameter is used for processing simulations of the Mars atmosphere Setting MARS 1 switches gravity gas constant and planet radius to the correct values for the planet Mars 5 1 9 MULTI The parameter MULTI can bes used to process a series of input data within one run of the pumaburner Setting MULTI to a number n tells the pumaburner to procees n input files The input files must follow one of the following two rules e YYMM rule The last four characters of the filename contain the data in the form YY MM e NNN rule The last four characters of the filename consist of a dot followed ny a 3 digit sequence number Examples Namelist contains MULTI 3 Command pumaburn lt namelist gt printout r
23. c puma_movie x program The working directory File Select working directory must be set to the directory with the formerly created movie files if not changed this is pumamovie When starting the program puma_movie offers the movie files one after the other for visualisation The model configuration file see part II section may include a special group movie described by the entries amp movie movie_namelist movie parameters waitTime 3 1000 0 false wait x msec per time step The waitTime Parameter tells the movie program the number of microseconds to wait before changing to the next snapshot The other entries in the model 59 configuration file have no effect on the movie run except for the progressBar entry in the gui group which is recommended to be changed to a suitable value for the movie run default 200 msec The value should leave the viewer enough time to look at each displayed graphic 4 0 5 2 Creating a Reference Movie data can also be used to create a reference file by the puma_reference tool The program examples exec puma_reference x converts a number of input files like puma_var nc into a file puma_reference nc containing the monthly means for one year of data The number of input files must be dividable by 12 The provided file examples ref puma_reference nc contains the monthly means of the third year of a PUMA run based on the provided model configuration file The current version of the puma_reference
24. continue execution with the previous parameters With OK the selection is confirmed and a window with a table appears which serves for further specification of the selected parameters Beside the name the standard value initialized by the table value of the static parameters the minimum and maximum values initialized by 90 and 110 of the standard value the increment initialized as 1 of the standard value and the maximum change per time step Delta that so far always is equal to the increment are denoted Any data except the name are changeable whereby it is to be noted that the conditions Min lt Default lt Max and Step lt Delta lt The PUMA namelist parameter nrun or ndays must either be the total amount of simulation time steps or days to simulate or has to be set to negative The last case causes an infinite simulation loop until it is ended by the GUTs control panel 47 P ete ectio AE r Select up to 4 parameters radiation parameters vj co2 concentration ppm co2 _ solar constant gsol0 land parameters min albedo for land albland roughness length land dzOland a Ni wetness factor land drhsland u 7 parameters max albedo for sea ice albice albedo for free ocean albsea wetness factor ice drhsice roughness length ice dzOice vi roughness length sea dz0sea OK Reset Cancel Figure 4 6 Selecting Interacti
25. e Ice coupling 32 time steps SSTFILE CHAR 80 file containing climatol ogy surface_parameter Structure Internally seamod f90 uses the FORTRAN 90 module seamod which uses the global common module pumamod from pumamod f90 Subroutine seaini reads the namelist and if the parallel version is used distributes the namelist parameters to the different processes If the run is not started from a restart file NRESTART from namelist inp is 0 the sea surface temperature and the ice cover is read from the surface_parameter file Ice thickness is computed from ice cover Additionally mixed layer depth and the 400 m Levitus temperature is read from the file ocean_parameter Climatology and namelist information is passed to the ice and ocean modules via the external subrou tines iceini in icemod f90 and oceanini in oceanmod f90 or oceanmod50 f90 Every NCPL_ATMOS_ICE time steps seastep calls the ice module via the external subroutine cplexchange_ice defined in intermod_atm f90 At the end of the integration seastop writes the restart information into file sea_restart 26 CHAPTER 2 MODULES 2 7 2 intermodatm f90 General The module intermod_atm f90 contains subroutines that exchange information between the atmospheric module and the sea ice module If an external coupler is used with an independent sea ice ocean model the module is replaced e g by mpccimod_atm f90 which contains the relevant subroutines for the MpCCI c
26. e data directory cp EXPDIR model x model x HHH HHH HHH II 1 while II le 2 do INMON cat EXPDIR savemon EXP INYEAR cat EXPDIR saveyear EXP INDAY cat EXPDIR saveday EXP INFTP cat EXPDIR saveftp EXP TARFILE cat EXPDIR savetar EXP YY INYEAR 36 CHAPTER 3 RUNNING PUMA if INYEAR 1t 10 then YY O YY fi MM INMON if INMON 1t 10 then MM 0 MM fi make namelist EXPDIR NAMLIST exe 1 cp EXPDIR land_namelist land_namelist cp EXPDIR sea_namelist sea_namelist cp EXPDIR ice_namelist ice_namelist cp EXPDIR ocean_namelist ocean_namelist cp EXPDIR EXP RES puma_restart cp EXPDIR EXP LANDRES land_restart cp EXPDIR EXP SEARES sea_restart model x gt EXPDIR EXP PROUT_ YY S MM history and restart saved for further diagnostics mv puma output EXP PUMA_ YY MM if INFTP eq 1 then TARFILE EXP TAR_ YY MM echo TARFILE gt EXPDIR savetar EXP tar cf EXPDIR TARFILE EXP PUMA_ YY MM else tar rf EXPDIR TARFILE EXP PUMA_ YY MM fi mv puma_restart EXPDIR EXP RES mv land_restart EXPDIR EXP LANDRES mv sea_restart EXPDIR EXP SEARES rm EXP PUMA_ YY MM OUTMON expr INMON 1 OUTDAY expr INDAY
27. e in a plot which is illustrated in the following script It plots seasonal means of the sea level pressure The data file 68 CHAPTER 6 GRAPHICS is prepared like this srv selcode 151 puma srv slp srv srv seasmean slp srv slp_sm srv srv2gra slp_sm srv The commands set vpage sets virtual pages inside the graphic window The full window is 11 inch wide and 8 5 inch high so set vpage 0 5 5 4 25 8 5 defines the upper left corner If setlevs 1 is specified the pressure levels as given are used Otherwise GrADS defines contour levels depending on the data set slp_sm gs setlevs 1 reinit open slp_sm enable print print mf set vpage 0 5 5 4 25 8 5 set gxout contour if setlevs 1 set clevs 990 995 1000 1005 1010 1015 1020 endif set ccols 1 set grads off set t 1 d c151 100 draw title SLP hPa yr ny DJF set vpage 5 5 11 4 25 8 5 set gxout contour if setlevs 1 set clevs 990 995 1000 1005 1010 1015 1020 endif set ccols 1 set grads off set t 2 d c151 100 draw title yr ny MAM set vpage 0 5 5 0 4 25 set gxout contour if setlevs 1 set clevs 990 995 1000 1005 1010 1015 1020 endif 6 2 VIS5D 69 set ccols 1 set grads off set t 3 d c151 100 draw title yr ny JJA set vpage 5 5 11 0 4 25 set gxout contour if
28. e temperature K Ti oceanic temperature profile K To melting temperature 273 16 K To reference temperature profile 250 0 K U scaled zonal wind u cosy u zonal wind ms Ux friction velocity ms V scaled meridional wind v cosy v meridional wind ms Uv horizontal wind vector ms W transmitted water vapor path kg m z vertical coordinate m Zo roughness length m At time increment 75 Symbol Definition Value Unit a thermal expansion coefficien T Jato KT scaled vorticity 0 potential temperature K K z Ra Cpa kappa mean heat conductivity in ice and snow Wm K Ki heat conductivity in ice 2 03 Wm K Ks heat conductivity in snow 0 31 W mt K A longitude H sin p ko cosine of the solar zenith angle p density of air kg m Pi density of sea ice 920 kgm Ps density of snow 330 kgm Pw density of sea water 1030 kg m Po density of fresh water 1000 kgm o normalized pressure coordinate p ps ol vertical velocity in o system TF time scale for RF TR time scale for NC TT time scale for temperature flux correction s Th time scale for depth flux correction s geopotential height g z m s scaled geopotential height p latitude X scaled velocity potential Y scaled streamfunction Y angular velocity of the earth 7 292 1074 s7 76 APPENDIX A LIST OF CONSTANTS AND SYMBOLS Appendix B Puma Codes 78
29. elist parameter NICE is set to 1 the subroutine subice is called which calculates ice cover and thickness Otherwise climatologi cal data interpolated to the current time step by iceget are used If an ice cover is present the surface temperature is calculated in skintemp Otherwise the surface temperature is set to the sea surface temperature calculated by the ocean model Every NCPL_ICE_OCEAN defined in sea_namelist time steps the external subroutine cplerchange_ ocean defined in intermod_ice is called to pass the atmospheric forcing to and retrieve oceanic data from the ocean module oceanmod f90 The oceanic data is used for ice calculations in the next time step 2 7 SEA ICE AND OCEAN MODULES 29 2 7 5 oceanmod f90 General The module oceanmod f90 contains a mixed layer ocean model i e subroutines to compute sea surface temperature and mixed layer depth The interface to the main PUMA module is via the module icemod f90 given by the subroutine oceanstep which is called by cplexchange_ocean defined in intermod_ice Input Output oceanmod f90 requires the file ocean_flxcor if NFLX CORRSST or NFLXCORRMLD is set to a negative value If NOUTPUT is set to 1 the output file fort 31 containing global fields of ocean model data in service format is produced for details see the ice modul section of the ref erence guide The module is controlled by the namelist oceanpar in the file ocean_namelist Parameter Type Pu
30. et the standard values will be recovered and can be changed again whereas Cancel leads to rejection of the selected attitudes Subsequently you return to the main window al Model properties Name Value C ModelProperties a inp Z CJ miscpar a 4 fluxpar A Ciradpar A co2 concentration ppm co2 360 0 A solar constant gsol0 1367 0 A daily cycle of insolation ndcycle 0 a 4 rainpar A surfpar A ES landpar A C seapar A max albedo for sea ice albice 0 7 a albedo for free ocean albsea 0 069 7 wetness factor ice drhsice 1 0 A wetness factor sea drhssea 1 0 A roughness length ice dz0ice 0 05 A roughness length sea dz0sea 0 001 A ice model nice 1 7 ocean model nocean 1 9 Dai parts on progressBar progressBar 1000 number of timesteps betw visu steps writeStep 32 save snapshots for visualization makemovie false Giicepar Reset OK Cancel Figure 4 5 Setting Model Parameters 4 0 3 3 2 Selecting and Configuring Interactive Model Parameters The menu entry Options Select interactive parameters made for the interactive PUMA control leads to a list of parameters from which the interactive ones can be selected by clicking see fig 4 6 Up to four interactive parameters are allowed Reset deletes the whole choice and new parameters can be selected Cancel rejects the selection and returns to the main window to
31. here gt Ice Ice Ocean Ice cover Ice thickness Snow thickness Surface temperature Deep sea temperature z Mixed layer depth Net precipitation runoff Salinity Melt and freeze volume Heat fluxes d Heat fluxes dT Radiation Wind stress 1117 TEILT Tas III Table 2 1 Parameters to be exchanged between models Arrows denote the direction in which the parameter is passed e g the atmosphere receives ice cover information from the ice model 22 32 Timesteps Timesteps CHAPTER 2 MODULES surface temperature net precipitation ice cover runoff snow thickness total heatflux ice thickness sensible heatflux radiation wind stress sea surface temperature deep sea temperature mixed layer depth net precipitation salinity runoff total heatflux wind stress freeze and melt volume Figure 2 1 Schematic illustration of the model coupling 2 7 SEA ICE AND OCEAN MODULES FLOW DIAGRAM ATMOSPHERE ICE OCEAN EXCHANGE PUMA MAIN LOOP puma f90 LANDSTEP _ lt 7 landmod f90 SURFSTEP surfmod f90 DR SEASTEP seamod f90 CPLEXCHANGE_ICE intermod_atm f90 ICESTEP icemod f90 IN CPLEXCHANGE_OCEAN CPLEXCHANGE_ATMOS lintermod_ice f90 intermod_ice f90 OCEANSTEP oceanmod 90 Figure 2 2 Subroutine flow when no e
32. iguration file as you can save the configuration in a different file as derived from the previous configuration 4 0 3 3 Setting Optional Parameters The substantial component in the GUI concept are the model parameters which are defined in the model configuration file and which are provided for PUMA by the Fortran namelists see part II section There are two categories of parameters On the one hand there are static ones which remain constant during arun and on the other hand we have parameters that can be interactively changed during the run Before starting a simulation run therefore different settings have to be made e check up or redefinition of the static parameters e selection of the interactive parameters e specification of the interactive parameters These settings are made in the menu Options The execution of the simulation program may be redirected to a remote host by choosing the menu item Options Set Program Parameters A new window appears see fig 4 4 where a remote hostname may be entered after selecting Make Remote Connection After pressing OK the next run of the simulation executable will be started on the specified remote host Cancel rejects this specification This is realized by a 45 File Options Set model properties Select interactive parameters Set program parameters Logging parameters 3 Deviaton Temp Edit Visu Paramete Fi Temne
33. in the other rows neither the short title nor the formula or anything sensible else In the following table you find an overview about the seven elements of the dialogbox and some input examples Label Description Title Full title displayed in the GrADS window Short Title Short title displayed in the visualisation panel of the main window Formula With this formula GrADS displays the NetCDF data e g dt 273 15 or du dv hcurl du dv Level Atmospheric pressure level s e g 0 1000 or 0 1000 Longitude Global longitudes e g Oor 0 360 Latitude Global latitudes e g Oor 90 90 Color Level Level range of colors which GrADS will use e g 40 35 30 25 20 15 10 5 0 5 10 15 20 In the list for level longitude and latitude may only be entries with one or two numbers The third examples for the formula is more complex than the others du dv hcurl du dv The semicolon divides the variables The list for the color levels may include one or more numbers If the input for level longitude and latitude includes more than one number the numbers are divided with space There is a special addon to display the wind vectors There has to be the keyword wind in the title string of the displayed Gr ADS window and the formula mag u v must be used 4 0 3 5 Changing the environment 4 0 4 Trouble Shooting This section provides answers to frequently asked questions about the GUI 54 CHAPTER 4 GRAPHICAL USER INTERFACE 4 0 4 1 Where are the
34. ing data to the relevant variables defined in icemod 28 CHAPTER 2 MODULES 2 7 4 icemod f90 General The module icemod f90 contains subroutines to compute sea ice cover and thickness The interface to the main PUMA module is given by the subrou tine icestep which is called by cplexchange_ice defined in intermod_atm f90 which is called by seastep defined in seamod f90 Input Output icemod f90 requires the file ice_flxcor if NFLXCORR is set to a negative value If NOUTPUT is set to 1 the output files fort 75 containing global fields of ice model data and the file fort 76 containing diagnostic ice data are produced for details see the reference guide Both output files are in service format The module is controlled by the namelist cepar in the file ice_namelist Parameter Type Purpose default NDIAG INTEGER Diagnostic output every NDIAG 160 time steps NOUT INTEGER Model data output every NOUT 32 time steps NOUTPUT INTEGER Icemodel output 0 no 1 yes 1 NFLXCORR INTEGER Time constant for restoring gt 0 360d no flux correction 0 use flux correction from file lt 0 Structure icemod f90 uses the module icemod which is not dependent on the module pumamod Subroutine ceini reads the namelist and when re quired the flux correction from the file ice_flxcor Subroutine icestep calls cplerchange_atmos defined in intermod_ice to get the atmospheric forcing fields If the sea_nam
35. known files You will find more about logging in section 4 0 3 3 3 Logging Menu and in section 4 0 4 Troubleshooting 4 0 3 Operating 4 0 3 1 Starting and Terminating The GUI is started with bin plaisirgui in case you are in the PlaisirGUI main directory or with plaisirgui if the bin directory is located in the standard path Then the main window of the GUI appears figure PlaisirGUI can be termi nated either by the menu option File Exit Program or by the close icon of the main window border 4 0 3 2 Adjusting the Simulation Program and the Working Directory Figure shows the GUI file menu with six menu entries With the first three menu items you can choose the working directory the PUMA executable program and the model configuration file With the next two you can save the model configuration And the last file menu item closes the GUI program terminates the PUMA simulation if still running and the GrADS windows and saves the actual model configuration 43 File Options Interactive Visualisation _ Wind Force Edit visu parameters _ Mean Temperature Open all visu windows Temperature Close all visu windows bin m Running step 0 _ Pressure Simulation Figure 4 1 PlaisirGUI Main Window Selecting the simulation program and the working directory It is possible to select the working directory under File Select working directory All input
36. mally or whether parameter values have been changed In the first case data for visualization are written to a file In the latter case it is determined which parameter values have been changed Not all parameters are yet allowed to be changed via the GUI so if one of these parameters is changed the subroutine err_notimplemented is called displaying an appropriate message Otherwise the PUMA parameter value and associated variables are set to the new value 2 7 SEA ICE AND OCEAN MODULES 21 2 7 Sea ice and ocean modules This section describes the modules that represent sea ice and ocean and the necessary interfaces between these modules and the atmospheric modules Con ceptually the sea ice model lies inbetween the atmosphere model and the ocean model Thus the PUMA main part and the ocean model are both coupled to the sea ice model but not directly to each other The sea ice model decides whether a given gridpoint is covered with ice or not in the latter case it merely functions as passing the ocean fluxes to the atmosphere and vice versa The parameters that are exchanged are listed in Table The sea ice and ocean model use a time step of one day Thus atmospheric coupling to the sea ice model is per formed every 32 time steps while the sea ice and ocean model are coupled every time step The coupling scheme is shown in Fig 2 1 Fig shows how the subroutines are placed when no external coupler is used Parameter Atmosp
37. mulation only for a short interruption PAUSE M is pushed The latest shown visualisation windows remain visible To continue the simula tion the START button has to be pushed again The third button 1 serves for stopping the running simulation The lastly shown graphics remain visible At the bottom of the simulation panel is a progressbar which shows the actual step of the simulation If the namelist parameter nrun resp ndays is set to be negative which means that the running simulation is set to infinite loop the progressbar will restart on each 1000th step This increment can be changed with model parameter progressbar see section part II 4 0 3 4 1 Configuring the Visualisation The gui allows to visualise actually generated simulation output The file format used for the interim output is NetCDF format which can be read and displayed with GrADS It is possible to reconfigure the visualisation properties without any handling in the property files To configure the visualisation it is neccessary to know which data variables the simulation writes into the NetCDF files For the delivered sim ulation program lt gui_home_dir gt example run puma_scalar_static x there are seven variables in the NetCDF output 51 Visu Parameter visu Parameter rl Temp Europe a e Pressure T Es m Tempe Figure 4 10 Visualisation Dialogbox Variable Description Free coordinates po Surface pressure lon lat dt
38. needs the message passing interface MPI to be installed properly This package doesn t include MPI but it can be downloaded from many sites e g http www go dlr de fresh unix sre misc mpich 1 2 4 tar gz Run the PUMA executable by either 1 using a script like the runscript_template example script which has been tested on linux but may also work on other unix machines 2 just typing puma x in the run directory 3 creating your own run setup examples are given in chapter 4 2 A detailed description of model parameters to be defined in the namelists can be found in chapter 3 Finally analyse the output using the pumaburn4 program together with your own software If you like to plot some results you may use the srv2gra program which converts the output to a format suitable for GrADS see chapter 6 1 3 1 pumaburner As mentioned above using the postprocessor pumaburner4 is highly recom mended for extracting variables and levels from the raw packed puma data file puma data On platforms other then linux build the executable by just in voking a standard c compiler e g cc O2 o pumaburn4 pumaburn4 c Im Inetedf In case that netCDF output is not wanted or the netCDF library isn t available please change the line define NETCDF_OUTPUT into undefine NETCDF_ OUTPUT within the source and omit binding the netCDF library Inetedf at compile time When compiling the postprocess
39. nmod f90 Input Output If NFLXCORRSST is set to a negative value oceanmod50 f90 requires the file ocean_lgflxcor denoting long term gradient flux correction Otherwise the file heat_parameter is needed to calculate the flux correction If NOUTPUT is set to 1 the output file fort 31 containing global fields of ocean model data in service format is produced for details see the ice module section of the reference guide The module is controlled by the namelist oceanpar in the file ocean_namelist Parameter Type Purpose default NDIAG INTEGER Diagnostic output every NDIAG 480 time steps NOUT INTEGER Model data output every NOUT 32 time steps NOUTPUT INTEGER Oceanmodel output 1 0 no 1 yes NFLXCORRSST INTEGER Flag for calculating the flux cor 1 rection lt 0 or reading the flux correction from file gt 0 Structure The internal structure is exactly the same as in oceanmod f90 Chapter 3 Running PUMA 3 1 Interactive Console Mode Puma is started from a console by simply typing puma x The following files have to be present in the same directory puma_namelist land_namelist sea_namelist ice_namelist ocean_namelist surface_parameter ocean_parameter All settings like length of the integration special parameterizations etc given in the namelist files The parameter files contain the climatology When the integration is finished successfully the following files have been created p
40. ns nor together with the message passing interface In case of compiling the included tools such as srv2gra you have to add the unix compatibility library libU77 via 1U77 This enables the compiler to find the non standard fortran features such as getarg iarg etc Chapter 2 Modules This is a technical documentation of the PUMA II model In the following the purposes of the individual modules is given and the general structure and possible input and output opportunities namelist files are explained 13 14 CHAPTER 2 MODULES 15 2 1 fluxmod f90 General The module fluxmod f90 contains subroutines to compute the differ ent surface fluxes and to perform the vertical diffusion The interface to the main PUMA module puma f90 see is given by the subroutines fluzins fluxstep and flurstop which are called in puma f90 in the subroutines prolog gridpointd and epilog respectively Input Output fluxmod f90 does not use any extra input file or output file and is controlled by the namelist fluxpar which is part of the puma f90 namelist file Parameter Type Purpose Default NEVAP Integer Switch for surface evaporation 0 1 off 1 on NSHFL Integer Switch for surface sensible heat 1 flux 0 off 1 on NSTRESS Integer Switch for surface wind stress 0 1 off 1 on NVDIFF Integer Switch for vertical diffusion 0 1 off 1 on VDIFF_LAMM Real Tuning parameter for ver
41. optional Possible Values are 5 1 PUMABURNER 99 Setting Description Dimension for T21 resolution HTYPE S Spherical Harmonics 506 22 23 coefficients HTYPE F Fourier Coefficients 32 42 latitudes wavenumber HTYPE Z Zonal Means 32 levels latitudes levels HTYPE G Gauss Grid 64 32 longitudes latitudes 5 1 3 2 VTYPE VTYPE accepts the first character of the following string Following settings are equivalent VTYPE S VTYPE Sigma VTYPE Super Blanks and the equal sign are optional Possible Values are Setting Description Remark VTYPE S Sigma model levels Some derived variables are not available VTYPE P Pressure levels Interpolation to pressure levels 5 1 3 3 MODLEV MODLEV is used in combination with VTYPE If VTYPE is not set to Sigma the contents of MODLEV are ignored MODLEV is an integer array that can get as many values as there are levels in the model output The levels are numbered from top of the atmosphere to the bottom The number of levels and the corresponding sigma values are listed in the pumaburner printout The out putfile orders the level according to the MODLEV values MODLEV 1 2 3 4 5 produces an output file of five model levels sorted from top to bottom while MODLEV 5 4 3 2 1 sorts them from bottom to top 5 1 3 4 hPa hPa is used in combination with VTYPE P If VTYPE is not set to Pres sure the contents
42. or under linux please add the following option to the compile command D_FILE_OFFSET_BITS 64 This enables the handling of files gt 2 GByte which is still a hard limit under older linux kernels lt 2 4 20 For a detailed description of the pumaburner see chapter 5 12 CHAPTER 1 INSTALLATION 1 4 Intel FORTRAN compiler The PUMA model is tested with various brands of fortran compilers There is no special care to be taken when compiling it On default the INTEL compiler enables normal optimization O2 which should be sufficient and even recom mended for the model to run For your convenience just use the included Makefile for compiling see former sections In case of compiling the included tools such as srv2gra you have to add the compiler option Vaxlib This enables the compiler to find the non standard fortran features such as getarg iarg etc 1 5 Absoft FORTRAN compiler This compiler should work well with PUMA Just use the Makefile and watch PUMA being built see former sections The default for optimization is none The recommended level is O1 which corresponds to basic optimization This is sufficient for most applications Please change the Makefile accordingly or compile PUMA manually The next level of optimization may rearrange your code substantially including strength reduction loop invariant removal code hoisting and loop closure This is neither usable with debugging optio
43. oupler Input Output intermod_atm f90 does not use any extra input file or output file Structure The subroutines cplstart cplinit cplstop are dummy routines that are real subroutines only in the case of external coupling The subroutine cplerchange_ice which is called by seastep in module seamod f90 calls the external subroutine icestep defined in icemod f90 It then copies the ice ocean data to the relevant PUMA variables 2 7 SEA ICE AND OCEAN MODULES 27 2 7 3 intermodice f90 General The module intermod_ice f90 contains subroutines that exchange in formation between the sea ice module and the ocean and atmosphere module If an external coupler is used with an independent sea ice ocean model the module is replaced e g by mpccimod_ice f90 which contains the relevant subroutines for the MpCCI coupler Input Output intermod_ice f90 does not use any extra input file or output file Structure The subroutine cplexchange_ocean which is called by icestep in module icemod f90 calls the external subroutine oceanstep defined in ocean mod f90 if the sea_namelist entry NOCEAN is set to 1 Otherwise it calls the subroutine oceanget defined in oceanmod f90 which interpolates the cli matological values to the current time step It then returns the ocean data to the subroutine icestep The subroutine cplexchange_atmos which is also called by cestep in module icemod f90 copies the atmospheric forc
44. pecific heat of moist air at constant pressure Ike K7 Cpd specific heat of dry air at constant pressure 1005 46 Jkg K Cu specific heat of water vapor at constant pressure 1869 46 Jkg K Es specific heat of sea ice 2070 W skg t K Ce specific heat of snow 2090 W skg t K Cp specific heat of sea water 4180 W skg K7 Cw coefficient for the deep ocean heat flux 4 W m K Cw wetness factor D scaled divergence E Evaporation ms Eo extrateristical solar flux density Wm f Coriolis parameter 20 sin py 5 1 Fp tendency of the first moment e Km s Fy tendency of the zeroth moment 2 Kms F surface moisture flux kg m s7 Fr surface sensible heat flux W m Fa surface zonal wind stress Pa F surface meridional wind stress Pa g gravitational acceleration 9 81 m 73 SymbolDefinition Value Unit hmixe mixed layer depth m hmiz Climatological mixed layer depth m H effective mixed layer depth Tra Tey m H reduced center of gravity a m Ja vertical turbulent moisture flux kg m s7 Jr vertical turbulent temperature flux Km s Ju vertical turbulent flux of zonal momentum Pa Las vertical turbulent flux of meridional momentum Pa k von Karman constant 0 4 Kn exchange coefficient for heat Km exchange coefficient for momentum L latent heat Jkg Ly latent heat of fusion L L 3 28 10 Jkg L latent heat of fusion for sea ice 3 28 10 Jkg la mixing length for heat
45. provide the interface to land and the ocean modules respectively Input Output surfmod f90 reads the land sea mask and the orography from the surfaceparameter file see puma surfmod f90 is controlled by the namelist surfpar which is part of the puma f90 namelist file Parameter Type Purpose default NSURF Integer Debug switch not active NOROMAX Integer Resolution of orography NTRU OROSCALE Real Scaling factor for orography 1 Structure Internally surfmod f90 uses the FORTRAN 90 module surfmod which uses the global common module pumamod from pumamod f90 Subrou tine surfini reads the namelist and if the parallel version is used distributes the namelist parameters to the different processes If the run is not started from a restart file NRESTART from namelist np is 0 the land sea mask and the orography is read from the surfaceparameter file According to the namelist input the orography is scaled by OROSCALE transfered into spectral space and truncated to NOROMAX Calls to subroutines landini and seaini are the interfaces to the respective initialization routines contained in the land and ocean modules During the run the interface to land and ocean is given by calls to the external subroutines landstep and seastep which are called by surfstep At the end of the integration interface subroutines landstop and seastop are called by surfstop 18 CHAPTER 2 MODULES 2 4 mod f90 General Input Outp
46. puma output puma restart land restart and sea restart Further information on how to run puma is given in chapter 4 The file puma_output contains the model results and has to be postprocessed using pumaburn4 see chapter 5 The _restart files contain information neccessary to restart the model run from the end of the current integration 1 2 Unpacking and installation This chapter briefly describes which files are included how to install the pack age and which steps have to be taken for a first run The file plasim tar in cludes PUMA module files data and namelist files to run PUMA some additonal programs for postprocessing e g afterburner some scripts a Makefile and a README txt All files are listed below with a short description After unpack ing the tar file in a suitable directory with the command tar xvf plasim tar you should find the following files PUMA Modules file description fftmod f90 fast fourier transformation fluxmod f90 boundary layer fluxes and vertical diffusion surface fluxes of heat and momentum vertical diffusion of temperature moisture and momentum landmod f90 land surface and soil processes soil temperature soil water snow cover and snow temperature river runoff vegetation legmod 90 legendre transformation miscmod f90 miscellaneous correction of negative humidity mpimod f90 interface to MPI mpimod_dummy 90 dummy routines substituting the MPI routines
47. r spherical harmonics Fourier coefficients or zonal means on sigma levels 5 1 4 Format of output data The pumaburner supports three different output formats e GRIB GRiIdded Binary WMO standard for gridded data e NetCDF Network Common Data Format e Service Format for user readable data see below For more detailed descriptions see for example http www nws noaa gov om ord iob NOAAPORT resources Setting Description GRIB 1 NetCDF 0 The output file is written GRIB format This option can be used only for HTYPE Spherical Harmonics or HTYPE Gauss Grid GRIB 0 NetCDF 1 The output file is written in NetCDF format This op tion can be used for HTYPE Gauss Grid only GRIB 0 NetCDF 0 The output file is written in Service format This is the preferred format for user programs For a detailed description see the following section GRIB 1 NetCDF 1 Illegal combination 5 1 PUMABURNER 61 5 1 5 SERVICE format The SERVICE format uses the following structure The whole file consists of pairs of header records and data records The header record is an integer array of 8 elements head 1 ECMWF field code head 2 modellevel or pressure in Pa head 3 date yymmdd yymm00 for monthly means head 4 time hhmm or hh for HHMM 0 head 5 1 dimension of data array head 6 2 dimension of data array head 7 may be set with the parameter HEAD7 he
48. rature Figure 4 3 Options Menu Unix remote shell command rsh so for the successful start the working directory must be remotely mounted the user credentials must be set and the rhosts file must have an entry for the machine where the GUI is running El Parameter Dialog Remote Connection v Make Remote Connection Remote Hostname scllagin Figure 4 4 Program Parameters 4 0 3 3 1 Setting the Initial Model Parameters To initialize the model parameters after pressing Options Set model properties a window is obtained that shows all groups of model parameters the individual group branches can be opened by clicking and will then show a table containing the model parameters and their default values see fig 4 5 In this table the values of the parameters can be edited Values can be inserted both in scalar form and in vector form A comma separates the components in vectorial inputs In case one parameter is assigned to be interactive later on the entered value represents its standard value It is important to notice Any changes have to be confirmed either with lt RETURN gt or by selecting another form field double click in order to be finally transferred 6The simulation executable see section 4 0 3 2 may also be a script to start a parallel run 46 CHAPTER 4 GRAPHICAL USER INTERFACE With OK these settings will be accepted for the next PUMA run with Res
49. rpose default NDIAG INTEGER Diagnostic output every NDIAG 480 time steps NOUT INTEGER Model data output every NOUT 32 time steps NOUTPUT INTEGER Oceanmodel output 1 O n0 1 yes NFLXCORRMLD INTEGER Time constant for restoring 60d mixed layer depth gt 0 no flux correction 0 use fluxcorrec tion from file lt 0 NFLXCORRSST INTEGER Time constant for restoring sea 60d surface temperature gt 0 no flux correction 0 use fluxcor rection from file lt 0 Structure oceanmod f90 uses the module oceanmod which is not dependent on the module pumamod Subroutine oceanini reads the namelist and when required the flux corrections from the file ocean_flxcor Subroutine oceanstep calls mixocean which calculates mixed layer depth and temperature If an ice cover is present mixed layer depth is set to the climatological value and the sea surface temperature is set to the freezing temperature For details of the mixed layer model see the reference guide section ute 30 CHAPTER 2 MODULES 2 7 6 oceanmod50 f90 General The module oceanmod50 f90 contains a mixed layer ocean model with depth fixed to 50m and a SST fluxcorrection that does not stem from restoring to climatological data For details see section ute of the reference guide The module oceanmod50 f90 optionally replaces the module oceanmod f90 so the internal structure and the interface to the main PUMA module is identical to ocea
50. s negative humidity and an optional subroutine for dry convective adjustment is included in miscmod f90 Input Output miscmod f90 does not use any extra input or output file and is controlled by the namelist miscpar which is part of the puma f90 namelist file Parameter Type Purpose default NDCA Integer Switch for convective adjustment 0 0 off 1 on Structure Internally miscmod f90 uses the FORTRAN 90 module miscmod which uses the global common module pumamod from pumamod f90 Subrou tine miscini reads the namelist and if the parallel version is used distributes the namelist parameters to the different processes Subroutine miscstep calls the subroutine fizer to eliminate spurious negative humidity arising from the spectral method and if dry convection is switched on calls the subroutine mkdca to do the dry convective adjustment Subroutine miscstop is a dummy subroutine since there is nothing to do to finalize the computations in misc mod f90 17 2 3 surfmod f90 General The module surfmod f90 deals as an interface between the atmo spheric part of the model and modules or models for the land and the oceans The interface to the main PUMA module puma f90 see is given by the subroutines surfini surfstep and surfstop which are called in puma f90 in the subroutines prolog gridpointd and epilog respectively Calls to subroutines named landini landstep and landstop and seaini seastep and seastop
51. sea_namelist cp EXPDIR ice_namelist ice_namelist cp EXPDIR ocean_namelist ocean_namelist cp EXPDIR model x model x cp SSTFILE surface_parameter cp SURFFILE surface_parameter cp OCEANFILE ocean_parameter model x gt EXPDIR EXP PROUT_0101 mv puma_output EXP PUMA_0101 history and restart saved for further diagnostics tar cf EXPDIR TARFILE EXP PUMA_0101 mv puma_restart EXPDIR EXP RES mv land_restart EXPDIR EXP LANDRES mv sea_restart EXPDIR EXP SEARES echo OUTYEAR gt EXPDIR saveyear EXP echo OUTDAY gt EXPDIR saveday EXP echo OUTMON gt EXPDIR savemon EXP echo OUTFTP gt EXPDIR saveftp EXP echo TARFILE gt EXPDIR savetar EXP 3 2 BATCH MODE 39 cat runall to EXPDIR cat gt EXPDIR runall lt lt EOR usr bin ksh TMPDIR TMPDIR mkdir p TMPDIR cd TMPDIR set ex EXPDIR EXPDIR SCHAUER SCHAUER SCHAUERDIR SCHAUERDIR DATADIR DATADIR EXP EXP LASTYEAR LASTYEAR MONTHS MONTHS FTPINT FTPINT The puma namelist is generated with NRESTART 0 i e the first month is inte grated from climatology The script runall is generated which can be used to restart the run after an interruption The remainder of the script is a loop of one month integrations until the de sired integration time is reached After each year the monthly output is tarred together and moved to the schauer or th
52. setlevs 1 set clevs 990 995 1000 1005 1010 1015 1020 endif set ccols 1 set grads off set t 4 d c151 100 draw title yr ny SON print disable print Igxps c i print mf o slp_sm ps 6 2 Vis5D VisdD is a system for interactive visualization of large 5 D gridded data sets such as those produced by numerical weather models One can make isosurfaces contour line slices colored slices volume ren derings etc of data in a 3 D grid then rotate and animate the images in real time There s also a feature for wind trajectory tracing a way to make text annotations for publications support for interactive data analysis etc from the Vis5D home page lt http www ssec wisc edu billh vis5d html gt This powerful visualisation tool together with its documentation is available through the above home page Vis5D uses its own data format which makes it necessary to transform your data Depending on their format see also sec tion and the flowchart on lt http puma dkrz de puma download map gt you have the following choices If e your data is raw PUMA output you need to process it with the pumaburner postprocessor see section 70 CHAPTER 6 GRAPHICS in order to transform it to either NETCDF option n or namelist parameter NETCDF 1 or GRIB option g or namelist parameter GRIB 1 and proceed from there e your data is in SERVICE format you need
53. the selected GrADS window configuration In this pop down list another box item can be selected to display different NetCDF data with GrADS Selecting another box item will change all associated properties automatically If the second checkbox is selected GrADS shows the deviation of the monthly mean values to saved reference data The next button 2 opens a new dialogbox with more detailed properties of the GrADS configuration This dialogbox is described in detail in the next paragraph If the OK button on the bottom of the panel is pressed the actual changes of the settings are accepted and with the Cancel button these changes can be rejected Dialogbox for the Detailed Properties Figure shows the dialog box to configure visualisation parameters In this dialogbox there are more detailed configurable options for displaying the NetCDF data with GrADS Each pop down list is editable so you can add any new element to the list Removing with the Remove button deletes the actual item of the corresponding 93 pop down list in the row and replaces this item with the last element in the list With the pop down list you can change the actual entry by any other item in the list Rembember when changing anything in this dialogbox the box doesn t check any coherency with other items in the dialogbox This warning means for example if you change the title entry from temperature to pressure the dialogbox won t change anything else
54. ther machine than the graphical tool 4 0 2 Installing the PlaisirGUI This chapter describes the installation of the GUI The installation of the GUI and the GUI itself were so far tested only under SuSE Linux and this installation is written for an Linux installation Installation and use of the GUI require the following software e Perl 5 6 39 40 CHAPTER 4 GRAPHICAL USER INTERFACE e Java RE 1 4 0 we recommend Java RE 1 4 2 and e GrADS 1 8 Perl 5 6 is installed directly by most Linux distributions in the developer pack ages Furthermore installation software for Perl is found in the web at CPAN or ActiveStatd If needed JRE Java Runtime Enviroment software can be down loaded from Blackdown or SUN and GrADS is to be obtained on the GrADS home pagg The installation archive of the GUI is named PlaisirGUI 1 0 install tar gz for the version 1 0 For installation the archive must be unpacked to an arbitrary directory gunzip PlaisirGUI 1 0 install tar gz uncompress the archive tar xf PlaisirGUI 1 0 install tar unpack it or tar xzf PlaisirGUI 1 0 install tar uncompress and unpack the archive if GNU tar is used After that the installation script must be started cd PlaisirGUI 1 0 perl setup or simply setup if Perl is located in usr bin At the beginning of the installation the script asks for the binary paths of Java2 Perl and GrADS and the library paths of GrADS respectivly as well as for the
55. tical dif 160 fusion VDIFF_B Real Tuning parameter for vertical dif 5 fusion VDIFF_C Real Tuning parameter for vertical dif 5 fusion VDIFF_D Real Tuning parameter for vertical dif 5 fusion Structure Internally fluxmod f90 uses the FORTRAN 90 module flurmod which uses the global common module pumamod from pumamod f90 Subrou tine fluxini reads the namelist and if the parallel version is used distributes the namelist parameters to the different processes Subroutine fluxstep calls the subroutine surflx to compute the surface fluxes and calls the subroutine vdiff to do the vertical diffusion Subroutine fluxstop is a dummy subrou tine since there is nothing to do to finalize the computations in fluxmod f90 The computation of the surface fluxes in surflx is spitted into several parts After initializing the stability dependent transfer coefficients the subroutines mkstress mkshfl and mkevap are the computations which are related to the surface wind stress the surface sensible heat flux and the surface evaporation respectively 16 CHAPTER 2 MODULES 2 2 miscmod f90 General The module miscmod f90 contains miscellaneous subroutines which do not fit well to other modules The interface to the main PUMA module puma f90 see is given by the subroutines miscini miscstep and miscstop which are called in puma f90 in the subroutines prolog gridpointd and epilog respectively A subroutine to eliminate spuriou
56. tool has to be adapted if the variables written to puma_var nc change 56 CHAPTER 4 GRAPHICAL USER INTERFACE Chapter 5 Postprocessing 5 1 Pumaburner 5 1 1 Introduction The Pumaburner is a postprocessor for the Planet Simulator and the PUMA model family It s the only interface between raw model data output and diag nostics graphics and user software The output data of the Planet Simulator are stored as packed binary 16 bit values using the model representation Prognostic variables like temperature divergence vorticity pressure and humidity are stored as coefficients of spherical harmonics on o levels Variables like radiation precipitation evaporation clouds and other fields of the parameterization package are stored on Gaussian grids The tasks of the Pumaburner are e Unpack the raw data to full real representation e Transform variables from the model s representation to a user selectable format e g grids zonal mean cross sections fourier coefficients e Calculate diagnostic variables like vertical velocity geopotential height wind components etc e Transfrom variables from levels to user selectable pressure levels e Compute monthly means and standard deviations e Write selected data either in SERVICE GRIB or NetCDF format for fur ther processing 57 98 5 1 2 Usage CHAPTER 5 POSTPROCESSING pumaburn4 options InputFile OutputFile lt namelist gt printout option h
57. ul Internet site for reference and installation instructions is lt http grads iges org grads grads html gt Latest versions of GrADS can handle data in NETCDF format via the command sdfopen GRIB HDF SDS and in its native binary format The native format can conveniently be derived from SERVICE format In the following it is assumed that the PUMA output has been converted to SERVICE format with the pumaburner and the resulting file is called puma srv Monthly mean data is either obtained directly from the pumaburner namelist parameter MEAN 1 see section or via a PINGO command srv monmeans puma srv puma_m srv Information on the PINGO package can be found in DKRZ report 11 at lt http www mad zmaw de Pingo repdl html gt The SERVICE file has to be converted to GrADS s native format by the command srv2gra puma_m srv which results in the files puma_m gra and puma_m ctl The first file contains the data the latter one information on the grid time steps and variable names The program srv2gra is one of the postprocessing tools available at 65 66 CHAPTER 6 GRAPHICS lt http puma dkrz de puma download map gt If you chose to compile it yourself please read the comments in the first few lines of the program text Sometimes the srv2gra tool has difficulties to calculate an appropriate time increment from the date headers of the data records so you should check this In this example the file puma m ct1 should look like this
58. uma_output puma_restart land_restart sea_restart The file puma_output contains the model results and has to be postprocessed using the pumaburner cf Chapter The _restart files contain information necessary to restart the model run from the end of the current integration 31 32 CHAPTER 3 RUNNING PUMA 3 2 Batch Mode For long integrations it is more useful to run puma in batch mode i e start puma by calling a script that manages the model run The following script does just that Since it is quite long it is here split to parts with explanations inbetween usr bin ksh stime date This script runs the atmospheric model PUMA on a linux machine EXP example EXPERIMENT IDENTIFIER EXPDIR castor home user EXP EXPERIMENT DIRECTORY MODEL EXPDIR puma x THE MODEL EXECUTABLE SCHAUER 1 TRANSFER OUTPUT TO SCHAUER 1 YES SCHAUERDIR pf u user_account puma EXP U TREE DIRECTORY FOR OUTPUT DATADIR EXPDIR data OUTPUT DIRECTORY SSTFILE EXPDIR surface_parameter INITIAL DATAFILE PUMA SURFFILE EXPDIR surface_parameter INITIAL DATAFILE PUMA OCEANFILE EXPDIR ocean_parameter INITIAL DATAFILE OCEANMOD LASTYEAR 50 LAST YEAR TO BE SIMULATED FTPINT 12 MONTHS PER TAR FILE TMPDIR MFHOME tmpdir run mkdir p TMPDIR cd TMPDIR set ex mkdir p EXPDIR mkdir p DATADIR OUTYEAR 1 OUTDAY 30 OUTMON 2 OUTFTP 2 TARFILE EXP TAR_0101 cp MODEL EXP
59. un 005 out pumaburn processes the files lt run 005 gt lt run 006 gt lt run 007 gt Namelist contains MULTI 4 Command pumaburn lt namelist gt printout exp0211 out pumaburn processes the files lt exp0211 gt lt exp0212 gt lt exp0301 gt lt exp0302 gt 5 1 10 Namelist example VTYPE Pressure HTYPE Grid CODE 130 131 132 hPa 200 500 700 850 1000 MEAN 0 GRIB 0 NETCDF 0 This namelist will write Temperature 130 u 130 and v 131 on pressure levels 200hPa 500hPa 700hPa 850hPa and 1000hPa The output interval is the same as found on the model data e g every 12 or every 6 hours MEAN 0 The output format is SERVICE format 5 1 PUMABURNER 63 5 1 11 Troubleshooting If the pumaburner reports an error or doesn t produce the expected results try the following e Check your namelist especially for invalid codes types and levels e Run the pumaburner in debug mode by using the option d Example pumaburn lt namelist gt printout d data in data out This will print out some details like parameters and memory allocation during the run The additional information may help to detect the problem e Not all combinations of HTYPE VTYPE and CODE are valid Try to use HTYPE Grid and VTYPE Pressure before switching to exotic parameter combinations 64 CHAPTER 5 POSTPROCESSING Chapter 6 Graphics 6 1 Grads In this section visualisation using the graphics package GrADS is described A usef
60. ut Parameter Type Purpose default Structure 2 5 FILES 2 5 Files 19 20 CHAPTER 2 MODULES 2 6 guimod f90 General The module guimod f90 figures as a link between the graphical user interface GUTI with its subroutines contained in visumod f90 and the PUMA modules The basic subroutines guiini guistep and guistop are called by the subroutine prolog the time loop inside the subroutine master and the subrou tine epilog respectively see module puma f90 Input Output guimod f90 redefines all PUMA namelists to register param eter value changes For example the namelist radpar defines a new namelist radpar_n via CO2_N CO2 GSOLO_N GSOLO etc Thus whenever the value of the solar constant is changed via the graphical interface the new value is first stored into the variable GSOLO_N After checking that the value change is allowed and in a sensible range it is then copied to GSOLO to take effect in the PUMA code Structure Internally guimod f90 uses the FORTRAN 90 module guimod which next to the _n parameters redefines all parameters with a suffix _d These parameters set the maximal amount the base variable may be changed per time step Also for all parameters another one with suffix _set contains information on whether this parameter has been changed via the GUI recently The subroutine guistep is called once every time step It first checks whether the model is supposed to be running nor
61. ve Model Parameters Update the parameters or press OK albland albice Figure 4 7 Specifying Interactive Model Parameters Max Min must be satisfied Further plausibility checks will be made by the simulation program itself Any changes have to be confirmed either with lt RETURN gt or by selecting another form field double click in order to be finally transferred Reset changes any values to the default values Back returns to the selection window Cancel neglects all changes and leads back to the main window With OK the defined values are saved and the main window is reopened 4 0 3 3 3 Logging 48 CHAPTER 4 GRAPHICAL USER INTERFACE The GUI program logs all standard and error output to the starting console and to a log file For logging the logging framework log4j P by the Jakarta Apache Project is used The logging menu Options Logging Parameters includes three entries in a sub menu The first item in this submenu Logging Parameters Log Window opens a logging window from log4j The logging configuration is changed automatically to verbose logging when opening the logging window The second entry is a menu checkbox called Logging Parameters Advanced Logging This checkbox activates only a verbose logging there will be no logging window opened The last menu entry is Logging Parameters Clear Logs and deletes all logging data in the actual log files In section
62. with plaisirgui from wherever you want e g the work ing directory If you want to reconfigure some binary or library paths for tests or permanently you can reset the following enviroment variables JAVA_BIN PERL_BIN GRADS_BIN GADDIR and GASCRP The installation script creates five new subdirectories in the target directory The directory bin contains a link to a Perl script starting the GUI The directory examples contains some example files If you choose a different working directory as our run example remember to copy the contents of the example run directory into your favorite working directory Otherwise it is pos sible that PUMA stops running In the directory lib the necessary Java byte code files are located The directory scripts contains the Perl scripts of the program The standard output and the error messages of the simulation program e g PUMA and of GrADS will be written into log files Normally there are three log files in the directory where you start the GUI 42 CHAPTER 4 GRAPHICAL USER INTERFACE e plaisirgui log e plaisirgui grads log e plaisirgui puma log The first log file includes the logging data from the GUI program the second file includes the logging of GrADS and the third one of PUMA If the logging was not deleted before starting the GUI the old log created on the last running day the files will be renamed by adding the date to the file extension and the new logging will be written into the
63. xternal coupler is used 23 24 CHAPTER 2 MODULES 2 7 SEA ICE AND OCEAN MODULES 25 2 7 1 seamod f90 General The module seamod f90 deals as an interface between the atmospheric part of the model and modules for the ocean and sea ice The basic subroutines seaini seastep and seastop are called by the the subroutines surfini surfstep and surfstop respectively see module surfmod f90 See the reference guide section heiko coupling for a visualization of the module coupling structure Input Output seamod f90 needs the parameter file surface_parameter to read the climatological sea surface temperature and ice cover as well as ocean_parameter which contains climatological mixed layer depth and the Levitus 400 m temperature As output data the file sea_restart is produced at the end of a run In the case of a restart this file is required to be read in by the module The namelist seapar which is contained in the file sea_namelist is defined as Parameter Type Purpose default ALBSEA REAL Albedo for open water 0 069 ALBICE REAL Max albedo for sea ice 0 7 DZOSEA REAL Roughness length sea 1 5 10 m DZOICE REAL Roughness length ice 1 0 107m DRHSSEA REAL Wetness factor sea 1 0 DRHSICE REAL Wetness factor ice 1 0 NOCEAN INTEGER Ocean model 1 or cli 1 matology 0 NICE INTEGER Sea ice model 1 or cli 1 matology 0 NCPL_ICE_OCEAN INTEGER Ice Ocean coupling time 1 steps NCPL_ATMOS_ICE INTEGER Atmospher

Planet Simulator User's Guide

Contents

Download Pdf Manuals

Related Search

Related Contents