IMPORTANT CHANGE: Removed all reference/use of ngridmx (dimphys.h) in
routines (necessary prerequisite to using parallel dynamics); in most
cases this just means adding 'ngrid' as routine argument, and making
local saved variables allocatable (and allocated at first call).
In the process, had to convert many *.h files to equivalent
modules: yomaer.h => yomaer_h.F90 , surfdat.h => surfdat_h.F90 ,
comsaison.h => comsaison_h.F90 , yomlw.h => yomlw_h.F90 ,
comdiurn.h => comdiurn_h.F90 , dimradmars.h => dimradmars_mod.F90 ,
comgeomfi.h => comgeomfi_h.F90, comsoil.h => comsoil_h.F90 ,
slope.h => slope_mod.F90
Also updated EOF routines, everything is now in eofdump_mod.F90