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
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