WRF ARW User's Guide - MMM - UCAR

POST-PROCESSING
Add the markers NCLFORTSTART and NCLEND to the subroutine as indicated
below.
Note, that
local variables are outside these block markers.
FORTRAN 77 routine called myTK.f, with NCL markers added
C NCLFORTSTART
subroutine compute_tk (tk, pressure, theta, nx, ny, nz)
implicit none
C Variables
integer nx, ny,
nz
real tk( nx,ny,nz) , pressure(nx,ny,nz), theta(nx,ny,nz)
C NCLEND
C Local Variables
integer i, j, k
real pi
DO k=1,nz
DO j=1,ny
DO i=1,nx
pi=(pressure(i,j,k)
/ 1000.)**(287./1004.)
tk( i,j,k) = pi*theta(i,j,k)
ENDDO
ENDDO
ENDDO
return
end
Now compile this code using the NCL script WRAPIT.
NOTE: If WRAPIT
cannot be found, make sure the environment variable
NCARG_ROOT
has been set correctly.
If the subroutine compiles successfully, a new library will be
created, called myTK.so.
This
library can be linked to an NCL
script to calculate TK. See how this is done in the
example below:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl” external myTK "./myTK.so"
begin
WRAPIT myTK.f
t = wrf_user_getvar (a,”T”,5)
theta = t + 300
p = wrf_user_getvar (a,”pressure”,5)
WRF-ARW V3: User’s Guide 9-17