/[lmdze]/trunk/Sources/dyn3d/PPM3d/ppm3d.txt
ViewVC logotype

Contents of /trunk/Sources/dyn3d/PPM3d/ppm3d.txt

Parent Directory Parent Directory | Revision Log Revision Log

Revision 166 - (show annotations)
Wed Jul 29 14:32:55 2015 UTC (8 years, 10 months ago) by guez
File MIME type: text/plain
File size: 2223 byte(s) 
Split ppm3d.f into files containing a single procedure.

Factorized computations of filtering matrices into a procedure
inifilr_hemisph. Had then to change the matrices from allocatable to
pointer and from customized lower bound to lower bound 1. The change
in lower bounds does not matter because the matrices are only used as
a whole as actual arguments.

Also, in infilr, instead of finding jfilt[ns][uv] from approximately
jjm /2, start at index j1 that corresponds to the equator. This is not
the same if there is a zoom in latitude.

Also, the test (rlamda(modfrst[ns][uv](j)) * cos(rlat[uv](j)) < 1) in
the loops on filtered latitudes is not useful now that we start from
j1: it is necessarily true. See notes.

Just encapsulated lwvn into a module and removed unused argument ktraer.
1
2 ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/ppm3d.F,v 1.1.1.1 2004/05/19
3 ! 12:53:07 lmdzadmin Exp $
4
5
6 ! From lin@explorer.gsfc.nasa.gov Wed Apr 15 17:44:44 1998
7 ! Date: Wed, 15 Apr 1998 11:37:03 -0400
8 ! From: lin@explorer.gsfc.nasa.gov
9 ! To: Frederic.Hourdin@lmd.jussieu.fr
10 ! Subject: 3D transport module of the GSFC CTM and GEOS GCM
11
12
13 ! This code is sent to you by S-J Lin, DAO, NASA-GSFC
14
15 ! Note: this version is intended for machines like CRAY
16 ! -90. No multitasking directives implemented.
17
18
19 ! ********************************************************************
20
21 ! TransPort Core for Goddard Chemistry Transport Model (G-CTM), Goddard
22 ! Earth Observing System General Circulation Model (GEOS-GCM), and Data
23 ! Assimilation System (GEOS-DAS).
24
25 ! ********************************************************************
26
27 ! Purpose: given horizontal winds on a hybrid sigma-p surfaces,
28 ! one call to tpcore updates the 3-D mixing ratio
29 ! fields one time step (NDT). [vertical mass flux is computed
30 ! internally consistent with the discretized hydrostatic mass
31 ! continuity equation of the C-Grid GEOS-GCM (for IGD=1)].
32
33 ! Schemes: Multi-dimensional Flux Form Semi-Lagrangian (FFSL) scheme based
34 ! on the van Leer or PPM.
35 ! (see Lin and Rood 1996).
36 ! Version 4.5
37 ! Last modified: Dec. 5, 1996
38 ! Major changes from version 4.0: a more general vertical hybrid sigma-
39 ! pressure coordinate.
40 ! Subroutines modified: xtp, ytp, fzppm, qckxyz
41 ! Subroutines deleted: vanz
42
43 ! Author: Shian-Jiann Lin
44 ! mail address:
45 ! Shian-Jiann Lin*
46 ! Code 910.3, NASA/GSFC, Greenbelt, MD 20771
47 ! Phone: 301-286-9540
48 ! E-mail: lin@dao.gsfc.nasa.gov
49
50 ! *affiliation:
51 ! Joint Center for Earth Systems Technology
52 ! The University of Maryland Baltimore County
53 ! NASA - Goddard Space Flight Center
54 ! References:
55
56 ! 1. Lin, S.-J., and R. B. Rood, 1996: Multidimensional flux form semi-
57 ! Lagrangian transport schemes. Mon. Wea. Rev., 124, 2046-2070.
58
59 ! 2. Lin, S.-J., W. C. Chao, Y. C. Sud, and G. K. Walker, 1994: A class of
60 ! the van Leer-type transport schemes and its applications to the moist-
61 ! ure transport in a General Circulation Model. Mon. Wea. Rev., 122,
62 ! 1575-1593.
63
64 ! ****6***0*********0*********0*********0*********0*********0**********72
65
  ViewVC Help
Powered by ViewVC 1.1.21