Changeset 642 for codes/icosagcm/devel/Python/test/py/RSW2_MPAS_W02.py
- Timestamp:
- 12/19/17 15:26:51 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/Python/test/py/RSW2_MPAS_W02.py
r641 r642 38 38 #scheme = time_step.RKn_simple(1,caldyn.bwd_fast_slow, dt) # forward Euler scheme 39 39 print dt, scheme.csjl, scheme.cfjl 40 step = unst.Caldyn_step(mesh,scheme) 41 step.setup_TRSW() 40 step = unst.caldyn_step_TRSW(mesh,scheme,nt) 42 41 43 42 u0 = Omega*radius/12. # cf Williamson (1991), p.13 … … 63 62 for i in range(20): 64 63 if True: 65 for j in range(nt): 66 step.next() 64 step.next() 67 65 gh,u = step.mass, step.u 68 66 print i, gh.shape, gh.min(), gh.max(), u.min(), u.max()
Note: See TracChangeset
for help on using the changeset viewer.