File indexing completed on 2022-02-22 06:09:51 UTC
view on githubraw file Latest commit 16cc32c7 on 2022-02-21 18:12:34 UTC
16cc32c739 Mart*0001 #include "COST_OPTIONS.h"
0002 #ifdef ALLOW_STREAMICE
0003 # include "STREAMICE_OPTIONS.h"
0004 #endif
fb1e8c9704 Dani*0005
0006 subroutine cost_test( myThid )
16cc32c739 Mart*0007
fb1e8c9704 Dani*0008
0009
0010
16cc32c739 Mart*0011
fb1e8c9704 Dani*0012
16cc32c739 Mart*0013
fb1e8c9704 Dani*0014 IMPLICIT NONE
0015
0016
0017 #include "SIZE.h"
0018 #include "EEPARAMS.h"
0019 #include "PARAMS.h"
0020 #include "DYNVARS.h"
0021 #include "GRID.h"
16cc32c739 Mart*0022 #include "cost.h"
fb1e8c9704 Dani*0023 #ifdef ALLOW_STREAMICE
0024 # include "STREAMICE.h"
0025 #endif
0026
0027
0028
0029 integer bi, bj
0030 integer myThid
0031
0032 #ifdef ALLOW_COST_TEST
0033
16cc32c739 Mart*0034 integer i, j
0035
fb1e8c9704 Dani*0036 INTEGER ILNBLNK
0037 EXTERNAL ILNBLNK
16cc32c739 Mart*0038
0039
0040
0041
fb1e8c9704 Dani*0042
16cc32c739 Mart*0043
fb1e8c9704 Dani*0044
16cc32c739 Mart*0045
0046
0047
0048
0049
0050
fb1e8c9704 Dani*0051
16cc32c739 Mart*0052
fb1e8c9704 Dani*0053
16cc32c739 Mart*0054
0055
0056
0057
fb1e8c9704 Dani*0058
0059
16cc32c739 Mart*0060
fb1e8c9704 Dani*0061
0062
16cc32c739 Mart*0063 do bj = myByLo(myThid),myByHi(myThid)
0064 do bi = myBxLo(myThid),myBxHi(myThid)
fb1e8c9704 Dani*0065 do j=1,sNy
0066 do i=1,sNx
0067
0068 if (streamice_hmask(i,j,bi,bj).eq.1.0) then
0069 objf_test (bi,bj) = objf_test (bi,bj) +
58ae2ae163 Patr*0070 & u_streamice(i,j,bi,bj)**2+v_streamice(i,j,bi,bj)**2+
0071 & h_streamice(i,j,bi,bj)**2
fb1e8c9704 Dani*0072 endif
0073
16cc32c739 Mart*0074 enddo
0075 enddo
0076 enddo
0077 enddo
fb1e8c9704 Dani*0078
16cc32c739 Mart*0079 #endif /* ALLOW_COST_TEST */
fb1e8c9704 Dani*0080
0081 RETURN
0082 END