/[MITgcm]/MITgcm/pkg/cost/cost_test.F
ViewVC logotype

Diff of /MITgcm/pkg/cost/cost_test.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by heimbach, Sun Mar 25 22:33:54 2001 UTC revision 1.4 by heimbach, Sun Mar 24 02:14:21 2002 UTC
# Line 30  C     myThid - Thread number for this in Line 30  C     myThid - Thread number for this in
30    
31  #ifdef ALLOW_COST_TEST  #ifdef ALLOW_COST_TEST
32  C     == Local variables  C     == Local variables
       _RS maskC (1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
33        _RL thetaRef        _RL thetaRef
34    
35        integer i, j, k        integer i, j, k
# Line 43  C     == Local variables Line 42  C     == Local variables
42        itlo = mybxlo(mythid)        itlo = mybxlo(mythid)
43        ithi = mybxhi(mythid)        ithi = mybxhi(mythid)
44    
45        iLocOut =  6        iLocOut =  80
46        jLocOut =  35        jLocOut =  30
47        kLocOut =  1        kLocOut =  1
48    
49  ce    some reference temperature  ce    some reference temperature
# Line 53  ce    some reference temperature Line 52  ce    some reference temperature
52  C--   Calculate mask for tracer cells  (0 => land, 1 => water)  C--   Calculate mask for tracer cells  (0 => land, 1 => water)
53        k=1        k=1
54    
       do bj = jtlo,jthi  
         do bi = itlo,ithi  
           do j=1,sNy  
             do i=1,sNx  
               maskC(i,j) = 1.  
               if (_hFacC(i,j,k,bi,bj) .eq. 0.) maskC(i,j)=0.  
             end do  
           end do  
         end do  
       end do  
   
55  C--   Calculate cost function on tile of this instance  C--   Calculate cost function on tile of this instance
56        do bj = jtlo,jthi        do bj = jtlo,jthi
57          do bi = itlo,ithi          do bi = itlo,ithi
# Line 73  C--   Calculate cost function on tile of Line 61  C--   Calculate cost function on tile of
61              ig = myXGlobalLo-1+(bi-1)*sNx+i              ig = myXGlobalLo-1+(bi-1)*sNx+i
62    
63              if ((ig .eq. iLocOut) .and. (jg .eq. jLocOut)) then              if ((ig .eq. iLocOut) .and. (jg .eq. jLocOut)) then
64  cph               write(0,'(a,3(x,i4),a,4(x,i4))')                 write(*,'(a,3(x,i4),a,4(x,i4))')
65  cph     &              'COST    ',ig,jg,kLocOut,' TILE ',i,j,bi,bj       &              'COST    ',ig,jg,kLocOut,' TILE ',i,j,bi,bj
66                 objf_test(bi,bj) = theta(i,j,kLocOut,bi,bj)                 objf_test(bi,bj) = theta(i,j,kLocOut,bi,bj)
67    cph(
68                   print *, 'ph-cost ', ig, jg, kLocOut,
69         &              theta(i,j,kLocOut,bi,bj)
70    cph)
71              endif              endif
72    
73              end do              end do
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22