*IDENT MYSGCODE
*D SGDIAGS.1348,SGDIAGS.1467
SUBROUTINE OPSG3D(LPNOW)
Add you own code in here, writing any data to channel a valid
channel number - see below. Code should be written using fortran 90.
Use the existing OPSG3D subroutine as a template to ensure you
access all the necessary variables and data.
RETURN
END
*READ /my_directory_pathname/name_of_file_with_code_changes
DUR/SF + 1 = CUT/SF times Nwhere N is an integer. So if the analysis is from days 0.0 to 29.75 ( a duration of 29.75 days), the sampling frequency is 0.25 days (i.e there are 120 samples) and the cutoff is 3.0 days, then n is 10.
0 1.0 5 5 0.0 0 0will compute the covariance of the zonal wind, while
B 0.5 5 5 0.5 6 6will compute the eddy kinetic energy with only real (i.e. above ground - since the land/air mask is zero underground) data being used in the summations and averages.
*IDENT MYTFCODE
*D TRNFLX.116,TRNFLX.206
SUBROUTINE F3OUT(KF,KPASS,PX,KLEV)
Add you own code in here, writing any data to channel a valid
channel number - see below. Code should be written using fortran 90.
Use the existing TRNFLX subroutine as a template to ensure you
access all the necessary variables and data.
RETURN
END
*READ /my_directory_pathname/name_of_file_with_code_changes
*IDENT MYODCODE
*D OPDYNMC.2399,OPDYNMC.2518
SUBROUTINE OPOD3D(LPNOW)
Add you own code in here, writing any data to channel a valid
channel number - see below. Code should be written using fortran 90.
Use the existing OPOD3D subroutine as a template to ensure you
access all the necessary variables and data.
RETURN
END
*READ /my_directory_pathname/name_of_file_with_code_changes
1.0 5 5 0.0 0 0will compute the product of the zonal wind with the zonal wind (zonal flux of zonal wind), while
1.0 5 18 1.0 6 18will compute the sum of the zonal and meridional fluxes of temperature.
*IDENT PRESUMO15 *READ XXXXX/presumo15bug.nupdHere, the first line ensures that all new code is given an edit identifier (a requirement of nupdate) while the second line will ensure that any current bugfixes (or `official' extensions to the current version of PRESUMO) will all be included in your modified source. Make sure that
XXXXXis altered accordingly. Finally, you will need to insert additional nupdate edit of the form
*READ /my_directory_pathname/name_of_file_with_code_changeswhere the specified pathname and file MUST reside on the local host.
*READ XXXXX/presumo15bug.nupd *IDENT SUMO15Here, the first line will ensure that any current bugfixes (or `official' extensions to the current version of SUMO) will all be included in your modified source, and the second line ensures that all new code is given an edit identifier (a requirement of nupdate) while the second line . Make sure that
XXXXXis altered accordingly.
*READ /my_directory_pathname/name_of_file_with_code_changeswhere the specified pathname and file MUST reside on the remote host.