Intermediate complexity atmospheric physics package developed by Franco Molteni.

---------------------------------------------------------------------------
10-26-2002: From Franco Molteni SPEEDY model (auden:~jmc/mitgcm/speedy_code) :
"Multi-decadal simulations using an atmospheric GCM with simplified
physical parametrizations.  I: Model formulation and climatology ",
by F. Molteni, Climate Dynamics, in press (2002).  
http://www.ictp.trieste.it/~moltenif/speedy-doc.html
---------------------------------------------------------------------------

original files extracted from the SPEEDY code (master/ver23) :

lrwxrwxrwx    1 jmc            16 Oct 26 14:56 atparam.h -> par_horres_t30.h
lrwxrwxrwx    1 jmc            12 Oct 26 14:57 atparam1.h -> par_verres.h
-rw-r-----    1 jmc          1498 Mar  8  2001 cls_inphys.h
-rw-r-----    1 jmc           640 Mar  8  2001 com_cnvcon.h
-rw-r-----    1 jmc           938 Mar  8  2001 com_forcing.h
-rw-r-----    1 jmc           573 Mar  8  2001 com_forcon.h
-rw-r-----    1 jmc           525 Mar  8  2001 com_lflags.h
-rw-r-----    1 jmc           321 Mar  8  2001 com_lsccon.h
-rw-r-----    1 jmc          1071 Mar  8  2001 com_physcon.h
-rw-r-----    1 jmc          3197 Mar  8  2001 com_physvar.h
-rw-r-----    1 jmc          2771 Mar  8  2001 com_radcon.h
-rw-r-----    1 jmc          1361 Mar  8  2001 com_sflcon.h
-rw-r-----    1 jmc           513 Mar  8  2001 com_vdicon.h
-rw-r-----    1 jmc          2197 Mar  8  2001 ini_inphys.f
-rw-r-----    1 jmc           249 Mar  8  2001 par_horres_t30.h
-rw-r-----    1 jmc            97 Mar  8  2001 par_verres.h
-rw-r-----    1 jmc          5107 Mar  8  2001 phy_convmf.f
-rw-r-----    1 jmc          6476 Mar  8  2001 phy_fordate.f
-rw-r-----    1 jmc          2142 Mar  8  2001 phy_lscond.f
-rw-r-----    1 jmc          7305 Mar  8  2001 phy_phypar.f
-rw-r-----    1 jmc         15391 Mar  8  2001 phy_radiat.f
-rw-r-----    1 jmc          2403 Mar  8  2001 phy_shtorh.f
-rw-r-----    1 jmc          6941 May  2  2001 phy_suflux.f
-rw-r-----    1 jmc          3825 Mar  8  2001 phy_vdifsc.f 

-----------------
mv ini_inphys.f --> phy_inphys.f
mv cls_inphys.h --> phy_const.h
 include "cls_inphys.h" with CPP #include "phy_const.h" in phy_inphys.F

phy_phypar.f --> phy_driver.F

* from Sigma to P : 
  change loop [DO K=1,NLEV --> DO K=1,kGrd(J) ] ; 
  SIG & DSIG : *PSG or /PSG where necessary
* rewrite stability condition using Pot.Temp. (phy_convmf.F phy_vdifsc.F)
* U,V_Stress computed in Dynamics directly on C-grid ; and DRAG is 
   computed in SUFLUX as a function of square surface wind speed (VsurfSq).

* split com_radcom.h in 2 files: com_radcom.h & com_radvar.h
* take arrays with horizontal index out off com_physcon.h & com_sflcon.h
   and put them into AIM_GRID.h

* low level physics S/R : phy_convmf.F  phy_radiat.F  phy_suflux.F
                          phy_lscond.F  phy_shtorh.F  phy_vdifsc.F
  no longer use common blocks for grid dependent (2D,3D) variables.
  instead, 2D & 3D arrays are passed as argument only.

---------------------------------------------------------------------------
