CHANGE LOG FOR UNIFIED 3DVAR CODE
-------------------------------------------
Author(s)  :  Dave Parrish, John Derber
Reviewer(s):  Russ Treadon                
Date       :  3/23/04                 

New CVS version/tag:   none


REASON FOR CHANGES
------------------
   - Remove bug in premature deallocation of ggrid_g3 array in regional mode
   - Clean up code.  
   - Simplify hopers for global runs
   - Change approximation for saturation q at lower pressures.  Change from
     calculation using 10hPa to a maximum of 5% atmospheric pressure.





EXPECTED DIFFERENCES
--------------------
   - prior to ggrid_g3 bug fix the regional analysis file was filled with junk
   - minor differences to fields due to different q sat. at higher levels
                                                       





GLOBAL TEST RESULTS
-------------------
Test case:   /nfsuser/g01/wx20rt/global_gsi/cases/global/sigma/t62l28/2004032212
Old code :   /nfsuser/g01/wx20rt/global_gsi/sorc/gsi_cvs/ 20040319_emc.change
Machine  :   IBM
Setup    :   12 mpi tasks, blocking=unlimited

Max memory
   old:  158700 Kbytes
   new:  158696 Kbytes

Wall clock
   old:  440.835214 seconds
   new:  430.938181 seconds

Output from final iteration of minimization
old:   0: penalty,grad ,a,b=   3   1 0.280488175507800304E+06 0.816405805249112746E+05 0.193098207243191369E-02 0.103884836063740194E+01
new:   0: penalty,grad ,a,b=   3   1 0.280491846344094025E+06 0.822264402116889396E+05 0.191058956171423510E-02 0.105338653541229688E+01





REGIONAL TEST RESULTS
---------------------
Test case:   /nfsuser/g01/wx20rt/global_gsi/cases/regional/wrfnmm/east_nest/2004021718
Old code :   /nfsuser/g01/wx20rt/global_gsi/sorc/gsi_cvs/ 20040319_emc.change
Machine  :   IBM
Setup    :   12 mpi tasks, blocking=unlimited

Max memory
   old: 293292 Kbytes
   new: 293244 Kbytes

Wall clock
   old: 1189.149560 seconds
   new: 1162.830139 seconds

Output from final iteration of minimization
old:   0: penalty,grad ,a,b=   3   1 0.173637932371515672E+05 0.557860940907702898E+05 0.956470667036052433E-03 0.172686821105438953E+00
new:   0: penalty,grad ,a,b=   3   1 0.173637932371515672E+05 0.557860940907702898E+05 0.956470667036052433E-03 0.172686821105438953E+00





FILES REMOVED
-------------
   None





FILES ADDED
-----------
   None





FILES MODIFIED
--------------
   genqsat.f90     - set maximum vapor pressure to be 5% total pressure
   glbsoi.f90      - deallocate ggrid_g3 after write_all for regional run
   guess_grids.f90 - deallocate ggrid_g3 in destroy_atm_grids only for global run
   hopers.f90      - add increment to spectral coefficients
   locatelat_reg.f90 - change name of regional background error stats from "wrfstat"
                       to "berror_stats".  name is now consistent with global mode
   pcgsoi.f90      - move update of spectral coefficients (global mode) into hopers
   read_guess.f90  - code added to generate q<0 (limq) statistics, CODE COMMENTED OUT
   read_lidar.f90  - change formating of write(6,8000) at end of routine 
   read_nmm_guess.f90 - code added to generate q<0 (limq) statistics, CODE COMMENTED OUT
   read_prepbufr.f90 - change formating of write(6,8000) at end of routine
   read_radar.f90  - change formating of write(6,8000) at end of routine
   read_wrf_nmm_guess.f90 - code added to generate q<0 (limq) statistics, CODE COMMENTED OUT


MAKEFILE CHANGES
----------------
   None





SCRIPT CHANGES
--------------
   Change local name of regional mode background error statistics from "wrfstat"
   to "berror_stats"




FIX FILE CHANGES
----------------
   Background error statistics should be calculated using modified qsat. 
   Little change expected.