pkg/fizhi/step_fizhi_fg.F

a456aa407c Andr*0001 #include "FIZHI_OPTIONS.h"
f98d2ec0f4 Andr*0002        subroutine step_fizhi_fg (myTime, myIter, myThid, dt)
e337e4ca8c Andr*0003 c-----------------------------------------------------------------------
                0004 c  Subroutine step_fizhi_fg - 'Wrapper' routine to advance
                0005 c        the physics state and make a 'first guess' at the new
                0006 c        value. At this point, increment with the physics
                0007 c        tendency only.
                0008 c        Also: Set up "bi, bj loop" and some timers and clocks here.
                0009 c Call: step_physics
                0010 c-----------------------------------------------------------------------
                0011        implicit none
                0012 #include "SIZE.h"
                0013 #include "fizhi_SIZE.h"
f4a0368053 Andr*0014 #include "fizhi_land_SIZE.h"
e337e4ca8c Andr*0015 #include "fizhi_coms.h"
e5b3b476cb Andr*0016 #include "gridalt_mapping.h"
e337e4ca8c Andr*0017 #include "EEPARAMS.h"
e5b3b476cb Andr*0018 #include "DYNVARS.h"
                0019 #include "GRID.h"
e337e4ca8c Andr*0020
3768927683 Andr*0021        integer myIter, myThid
                0022        _RL myTime
e337e4ca8c Andr*0023
613fa3996d Andr*0024        integer bi, bj
e337e4ca8c Andr*0025        integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2
f98d2ec0f4 Andr*0026        _RL dt
e337e4ca8c Andr*0027
5c1a81edc4 Andr*0028        _RL tempij(sNx,sNy)
                0029        integer i,j,L
                0030
                0031        idim1 = 1-OLx
                0032        idim2 = sNx+OLx
                0033        jdim1 = 1-OLy
                0034        jdim2 = sNy+OLy
                0035        im1 = 1
                0036        im2 = sNx
                0037        jm1 = 1
                0038        jm2 = sNy
e337e4ca8c Andr*0039
                0040        do bj = myByLo(myThid), myByHi(myThid)
                0041        do bi = myBxLo(myThid), myBxHi(myThid)
                0042
                0043 c  Step forward the physics state using physics tendencies
5c1a81edc4 Andr*0044         call step_physics(uphy,vphy,thphy,sphy,dt,idim1,idim2,
                0045      .      jdim1,jdim2,
e337e4ca8c Andr*0046      .      Nrphys,Nsx,Nsy,1,sNx,1,sNy,bi,bj,duphy,dvphy,dthphy,dsphy)
                0047
8598abfc35 Andr*0048       if(2.eq.1 )then
5c1a81edc4 Andr*0049       print *,' In step fizhi fg, new fizhi fields ',bi,' dt= ',dt
                0050       do L = 1,Nrphys
                0051        do j = jm1,jm2
                0052        do i = im1,im2
                0053         tempij(i,j) = uphy(i,j,L,bi,bj)
                0054        enddo
                0055        enddo
                0056 c      print *,' uphy at level ',l,' ',tempij
                0057       enddo
                0058       do L = 1,Nrphys
                0059        do j = jm1,jm2
                0060        do i = im1,im2
                0061         tempij(i,j) = vphy(i,j,L,bi,bj)
                0062        enddo
                0063        enddo
                0064 c      print *,' vphy at level ',l,' ',tempij
                0065       enddo
                0066       do L = 1,Nrphys
                0067        do j = jm1,jm2
                0068        do i = im1,im2
                0069         tempij(i,j) = thphy(i,j,L,bi,bj)
                0070        enddo
                0071        enddo
                0072        print *,' thphy at level ',l,' ',tempij
                0073       enddo
                0074       do L = 1,Nrphys
                0075        do j = jm1,jm2
                0076        do i = im1,im2
                0077         tempij(i,j) = sphy(i,j,L,bi,bj)
                0078        enddo
                0079        enddo
8598abfc35 Andr*0080        print *,' sphy at level ',l,' ',tempij
5c1a81edc4 Andr*0081       enddo
                0082       endif
                0083
8598abfc35 Andr*0084         call qcheck (idim1,idim2,jdim1,jdim2,Nrphys,Nsx,Nsy,im1,im2,
                0085      .       jm1,jm2,bi,bj,dpphys,sphy)
e5b3b476cb Andr*0086
8598abfc35 Andr*0087       if(2.eq.1 )then
                0088       print *,' In step fizhi fg after qcheck ',bi
                0089       do L = 1,Nrphys
                0090        do j = jm1,jm2
                0091        do i = im1,im2
                0092         tempij(i,j) = sphy(i,j,L,bi,bj)
                0093        enddo
                0094        enddo
                0095        print *,' sphy after qcheck at level ',l,' ',tempij
                0096       enddo
                0097       endif
5c1a81edc4 Andr*0098
e337e4ca8c Andr*0099        enddo
                0100        enddo
                0101
                0102        return
                0103        end