/[MITgcm]/MITgcm/pkg/seaice/cost_ice_test.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/cost_ice_test.F

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

revision 1.1 by heimbach, Thu Sep 1 05:34:30 2005 UTC revision 1.6 by dimitri, Thu Jun 5 13:27:36 2008 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "SEAICE_OPTIONS.h"  #include "SEAICE_OPTIONS.h"
5    
# Line 8  c     ================================== Line 9  c     ==================================
9  c     SUBROUTINE cost_ice_test  c     SUBROUTINE cost_ice_test
10  c     ==================================================================  c     ==================================================================
11  c  c
12  c     o Compute sea-ice cost function.  The following options can  c     o Compute sea-ice cost function.  The following options can be
13  c       be selected with data.cost variable cost_ice_flag:  c       selected with data.seaice (SEAICE_PARM02) variable cost_ice_flag
14  c  c
15  c     cost_ice_flag = 1  c     cost_ice_flag = 1
16  c     - compute mean sea-ice volume  c     - compute mean sea-ice volume
# Line 81  c     == external functions == Line 82  c     == external functions ==
82  c     == end of interface ==  c     == end of interface ==
83    
84        if ( myTime .GT. (endTime - lastinterval) ) then        if ( myTime .GT. (endTime - lastinterval) ) then
85           tempVar = 1. /           tempVar = 1. /
86       &             ( ( 1. + min(endTime-startTime,lastinterval) )       &             ( ( 1. + min(endTime-startTime,lastinterval) )
87       &             / deltaTClock )       &             / deltaTClock )
88    
89  cph(  cph(
90  c      print *, 'ph-ice B ', myiter, theta(18,24,1,1,1),        print *, 'ph-ice B ', myiter, theta(4,4,1,1,1),
91  c     &     area(18,24,1,1,1), heff(18,24,1,1,1)       &     area(4,4,1,1,1), heff(4,4,1,1,1)
92  cph)  cph)
93           if ( cost_ice_flag .eq. 1 ) then           if ( cost_ice_flag .eq. 1 ) then
94  c     sea-ice volume  c     sea-ice volume
# Line 212  c     assumed that measurements are AREA Line 213  c     assumed that measurements are AREA
213           endif           endif
214        endif        endif
215    
216    cph(
217          print *, 'ph-ice C ', myiter, objf_ice(1,1)
218    cph)
219    
220  #endif /* ALLOW_COST_ICE */  #endif /* ALLOW_COST_ICE */
221    
222        end        end
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22