CHANGE LOG FOR UNIFIED 3DVAR CODE ------------------------------------------- Author(s) : John Derber Reviewer(s): Russ Treadon Date : 11 May 2004 CVS tag : ncep-gsi-2004_04p1 REASON FOR CHANGES ------------------ Array xhatsm in htoper was incorrectly dimensioned. The correct dimension is (lat1*lon1*(nsig*6+4)). The code incorrectly declared the array to be of size (lat1*lon1*(nsig*6+2)). SPECIAL NOTES ------------- None EXPECTED DIFFERENCES -------------------- No difference in global and regional test cases. However, the code may fail to compile or run on other machines due to this error. GLOBAL TEST RESULTS ------------------- Test case: /nfsuser/g01/wx20rt/global_gsi/cases/global/sigma/t254l64/2004041512 Old code : /nfsuser/g01/wx20rt/global_gsi/sorc/gmao_cvs/global_anl New code : /nfsuser/g01/wx20rt/global_gsi/sorc/gsi_cvs/global_anl Machine : NCEP IBM SP (frost) Setup : 48 mpi tasks, blocking=unlimited Max memory old: 632124 Kbytes new: 632116 Kbytes Wall clock old: 1390.158451 seconds new: 1277.275531 seconds Output from final iteration of minimization old: 0: penalty,grad ,a,b= 3 1 0.180208145768173825E+06 0.294227658014251074E+04 0.124170306385001602E-02 0.226055222032754166E+01 new: 0: penalty,grad ,a,b= 3 1 0.180208145768173825E+06 0.294227658014251074E+04 0.124170306385001602E-02 0.226055222032754166E+01 REGIONAL TEST RESULTS --------------------- Test case: /nfsuser/g01/wx20rt/global_gsi/cases/regional/wrfnmm/central_nest/2004041512 Old code : /nfsuser/g01/wx20rt/global_gsi/sorc/gmao_cvs/global_anl New code : /nfsuser/g01/wx20rt/global_gsi/sorc/gsi_cvs/global_anl Machine : NCEP IBM SP (frost) Setup : 12 mpi tasks, blocking=unlimited Max memory old: 473472 Kbytes new: 473548 Kbytes Wall clock old: 1712.247852 seconds (slow due to system load) new: 1183.959781 seconds Output from final iteration of minimization old: 0: penalty,grad ,a,b= 3 1 0.472878317408652438E+05 0.365295681744780159E+05 0.482728236857825344E-02 0.119700105442843097E+00 new: 0: penalty,grad ,a,b= 3 1 0.472878317408652438E+05 0.365295681744780159E+05 0.482728236857825344E-02 0.119700105442843097E+00 FILES REMOVED ------------- None FILES ADDED ----------- None FILES MODIFIED -------------- htoper.f90 - increase dimension of xhatsm array from (lat1*lon1*(nsig*6+2)) to (lat1*lon1*(nsig*6+4)). Also increase the upper limit for loop which initializes this array to zero. These changes affected lines 60 and 85 of the routine. MAKEFILE CHANGES ---------------- None SCRIPT CHANGES -------------- None FIX FILE CHANGES ---------------- None