igsm/src_chem/chemcheck.F


#include "ctrparam.h"

!       ============================================================
!
!       CHEMCHECK.F:    Subroutine for checking negative values and
!                           forcing them to be passitive ones
!                                  of MIT Global Chemistry Model
!
!       ------------------------------------------------------------
!
!       Author:         Chien Wang
!                       MIT Joint Program on Science and Policy
!                               of Global Change
!
!       ----------------------------------------------------------
! 
!       Revision History:
!       
!       When    Who             What
!       ----    ----------      -------
!       081795  Chien Wang      rev.    
!       080100  Chien Wang      repack based on CliChem3 & add cpp
!       051804  Chien Wang      rev.
!
!       ==========================================================

      Subroutine chemcheck(x00,x11,dta)

#include "chem_para"
#include "BD2G04.COM"

      dimension x00  (nlon,nlat,nlev)

!       ----------------------------------------------------------
#if ( defined CPL_CHEM )

c ===
c === 081295: check for negative values:
c ===

        oneover3 = 1./3.

        do i=1,nlon
        do j=1,nlat

          if(x00(i,j,1).lt.0.0) then
            xmean =  max(0.0,
     &                  (x00(i,j,1) + x00(i,j,2))*0.5 )
            x00(i,j,1) = xmean
            x00(i,j,2) = xmean
          endif

        do k =2,nlev1
          if(x00(i,j,k).lt.0.0)then
             xmeanv = max(0.0,
     &                (x00(i,j,k) 
     &               + x00(i,j,k-1) + x00(i,j,k+1))*oneover3 )
             x00(i,j,k-1) = xmeanv
             x00(i,j,k)   = xmeanv
             x00(i,j,k+1) = xmeanv
          endif
        end do

          if(x00(i,j,nlev).lt.0.0) then
            xmean = max(0.0,
     &                       (x00(i,j,nlev)
     &                      + x00(i,j,nlev1))*0.5 )
            x00(i,j,nlev)  = xmean
            x00(i,j,nlev1) = xmean
          endif

        end do
        end do

#endif

      return
       end

1	jscott	1.1
2			#include "ctrparam.h"
3
4			! ============================================================
5			!
6			! CHEMCHECK.F: Subroutine for checking negative values and
7			! forcing them to be passitive ones
8			! of MIT Global Chemistry Model
9			!
10			! ------------------------------------------------------------
11			!
12			! Author: Chien Wang
13			! MIT Joint Program on Science and Policy
14			! of Global Change
15			!
16			! ----------------------------------------------------------
17			!
18			! Revision History:
19			!
20			! When Who What
21			! ---- ---------- -------
22			! 081795 Chien Wang rev.
23			! 080100 Chien Wang repack based on CliChem3 & add cpp
24			! 051804 Chien Wang rev.
25			!
26			! ==========================================================
27
28			Subroutine chemcheck(x00,x11,dta)
29
30			#include "chem_para"
31			#include "BD2G04.COM"
32
33			dimension x00 (nlon,nlat,nlev)
34
35			! ----------------------------------------------------------
36			#if ( defined CPL_CHEM )
37
38			c ===
39			c === 081295: check for negative values:
40			c ===
41
42			oneover3 = 1./3.
43
44			do i=1,nlon
45			do j=1,nlat
46
47			if(x00(i,j,1).lt.0.0) then
48			xmean = max(0.0,
49			& (x00(i,j,1) + x00(i,j,2))*0.5 )
50			x00(i,j,1) = xmean
51			x00(i,j,2) = xmean
52			endif
53
54			do k =2,nlev1
55			if(x00(i,j,k).lt.0.0)then
56			xmeanv = max(0.0,
57			& (x00(i,j,k)
58			& + x00(i,j,k-1) + x00(i,j,k+1))*oneover3 )
59			x00(i,j,k-1) = xmeanv
60			x00(i,j,k) = xmeanv
61			x00(i,j,k+1) = xmeanv
62			endif
63			end do
64
65			if(x00(i,j,nlev).lt.0.0) then
66			xmean = max(0.0,
67			& (x00(i,j,nlev)
68			& + x00(i,j,nlev1))*0.5 )
69			x00(i,j,nlev) = xmean
70			x00(i,j,nlev1) = xmean
71			endif
72
73			end do
74			end do
75
76			#endif
77
78			return
79			end
80