
'reinit'
'open press_24_00.ctl'
'set t 1 72'
'set lev 750'
'set lat 38.3540'
'g = 9.801'
'omega = 7.292e-5'
'pi = 3.1415265'
'phi = (pi * 38.3540) / 180'
'f = 2*omega*sin(phi)'

'set lon -7.4427'
'ZA = geopt'

'set lon -7.4200'
'ZB = geopt'

'dZ = ZA - ZB'
'dy = 2000'

'ug = -(1/f)*(dZ/dy)'

'set lon -7.4314'
'set lat 38.3450'
'ZC = geopt'

'set lat 38.3629'
'ZD = geopt'

'dZZ = ZC - ZD'
'dx = 2000'

'vg = (1/f)*(dZZ/dx)'





