[MITgcm-cvs] MITgcm/pkg/cheapaml CVS Commit

[bderembl at forge.csail.mit.edu]

Update of /u/gcmpack/MITgcm/pkg/cheapaml
In directory forge:/tmp/cvs-serv21571

Modified Files:
	cheapaml_readparms.F CHEAPAML.h cheapaml_fields_load.F 
	cheapaml.F 
Log Message:



modified files:
cheapaml_readparms.F
CHEAPAML.h
cheapaml_fields_load.F
cheapaml.F


-there was a bug in cheapaml_fields_load.F:
      _RL  cheaph0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  cheaph1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  cheapcl0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  cheapcl1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  cheaplw0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  cheaplw1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  cheappr0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  cheappr1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)

were not in the common block

- there was another bug in cheapaml_fields_load.F:
  the fields were loaded only if a solarfile was present

IF ( SolarFile .EQ. ' '  )

I removed this condition

- last bug: in the bulk formula for long waves, I added a ABS(qair)
  inside the SQRT() for the computation of xlwnet. qair is never
  supposed to be negative but it may happen.
  

-added a flag for input precip (usePrecip set to false by default)
and the corresponding precip file (cheap_prFile) and fields (CheapPrecip)