pom2k.change

大名鼎鼎的海洋数值计算模型-pom模型源代码及说明书和示例

***********************************************************************

 A log documenting versions and code changes in pom2k model  

 (Please see comments in pom2k.f for changes from pom98.f)

***********************************************************************

VERSION		COMMENTS

2004-04-12	Order of i,j loops replaced in calc. of dy in 
                file2ic (thanks to user A. Manda). 

2004-01-08	A new test case (iproblem=3) is added, were init. cond.
		are read from file generated by GRID-DATA/GRID.f
		for a simple coastal region with real input data.
		File is read in new subroutine file2ic.
		See runpom2k script to see how to run diff. cases.

2003-09-24	profqWAVES is now the default profq in pom2k.f  
		This version fixes some problems of previous version 
		(thanks to C. Tang) and includes new BCs for q2.
		The inclusion of waves may not have much effect
		in the pom2k.f config. unless vertical res. increases. 

2003-09-09	In main code, vertvl & pom2k.c replace vflux by 
		vfluxb & vfluxf to allow time dep. surf. vertical vel.
		forcing (following Newberger & Hunter).
		Change should not affect most cases [when w(i,j,1)=0].

2003-05-08	The part of profqWAVES.sub pertaining to the length 
		scales, l, & l0, has been relocated in the code as 
		suggested by Priscilla Newberger. It is more consistant 
		with the descibing paper, but the change in the test cases 
		are insignificant. 

2003-04-09	A subroutine profqWAVES.sub is provided here.
		It is an enhanced M-Y turbulence profq to include wave
		breaking following the Craig-Banner scheme (see Mellor
		& Blumberg, 2003). It is not yet implemented in
		pom2k.f until tests are completed.

2002-09-06	Following J. Hunter's advice, loops in advave, advct,
		profq, vertvl changed in their range to make the model
		symmetric in values on boundaries. Also, unused parameters
		lij, lijk... etc., removed from pom2k.c. For detail, see
                2002-09-06_Hunter email in USERScomments dir. 

2002-07-11	Following a suggestion by J. Hunter, "dh" replaced by "h"
		in "profq" (for calc. of p= and boygr=) to be consistent 
		with a change in "dens" made in pom98 in Sep. 2001. 

2002-06-21	Cosmetic changes in prxy, prxyz, prxz, pryz printing S.R. 
          	Scale section before "write(6,1) label" simplified.
                (Should not affect results or printout at all.) 

2002-05-01	Minor changes: 1. e_atmos leftovers removed from "bcond". 
		2. art(i,j) moved from "box" & "seamount" to "areas_masks".
		3. in "prxy", "prxz" etc. set amx=1.e-12 instead of 0.0
		(if scale=0 AND all elemrnts of printed field=0, log(0)) 

2002-04-18	The surf. invert barometer effect e_atmos(i,j) has been
		deletet from the pressure calc. in sub. DENS (since elev.
		is already ignored). The way it was previously written was
		not correct. e-atmos also deleted from pressure calc. 
		in PROFQ (for consistency) and from diagnostic CFL condition. 
		(Note that OzPOM version retains e_amos as well as e_atmos0 
		following Robertson et al., JTECH 2001 correction).

2002-03-26	From this date on pom2k will be maintained by Princeton 
                while a version more confined to the original code,
            	named OzPOM is maintained by John Hunter at:
 		http://www.antcrc.utas.edu.au/~johunter/ozpom.html

		In the Princeton's version of pom2k about 300 lines of 
                comments have been removed. A new switch to disable the netCDF 
		output is now added by setting netcdf_file='nonetcdf'. 
 
2002-03-18	period=(2.e0*pi)/(abs(cor(im/2,jm/2))/86400.e0) replaced by 
           	period=(2.e0*pi)/abs(cor(im/2,jm/2))/86400.e0 as in pom98
		(Bug reported by Mingqiang FANG)

2002-02-26      A bug in printout of variable "lramp" (introduced by Princeton)
		has been fixed- format i10 replaced by l10.

2002-02-02      Code made available to users through pom web page. 
		Changes to original code include modifications of surf. 
		BC for W, format changes to printout and reducing 
		lines of comments for parameters setting.
 
2001-01-21	Original code delivered from J. Hunter to Princeton.