/[MITgcm]/MITgcm/pkg/kpp/kpp_routines.F
ViewVC logotype

Diff of /MITgcm/pkg/kpp/kpp_routines.F

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

revision 1.10 by heimbach, Sat Jul 13 03:12:30 2002 UTC revision 1.11 by mlosch, Wed Sep 18 16:38:02 2002 UTC
# Line 1340  c     work1, work2 - work arrays for hol Line 1340  c     work1, work2 - work arrays for hol
1340  c calculate density, alpha, beta in surface layer, and set dbsfc to zero  c calculate density, alpha, beta in surface layer, and set dbsfc to zero
1341    
1342        call FIND_RHO(        call FIND_RHO(
1343       I     bi, bj, ibot, itop, jbot, jtop, 1, 1, eosType,       I     bi, bj, ibot, itop, jbot, jtop, 1, 1,
1344       I     theta, salt,       I     theta, salt,
1345       O     WORK1,       O     WORK1,
1346       I     myThid )       I     myThid )
1347    
1348        call FIND_ALPHA(        call FIND_ALPHA(
1349       I     bi, bj, ibot, itop, jbot, jtop, 1, 1, eosType,       I     bi, bj, ibot, itop, jbot, jtop, 1, 1,
1350       O     WORK2 )       O     WORK2 )
1351    
1352        call FIND_BETA(        call FIND_BETA(
1353       I     bi, bj, ibot, itop, jbot, jtop, 1, 1, eosType,       I     bi, bj, ibot, itop, jbot, jtop, 1, 1,
1354       O     WORK3 )       O     WORK3 )
1355    
1356        DO J = jbot, jtop        DO J = jbot, jtop
# Line 1368  CHPF$  INDEPENDENT, NEW (RHOK,RHOKM1,RHO Line 1368  CHPF$  INDEPENDENT, NEW (RHOK,RHOKM1,RHO
1368        DO K = 2, Nr        DO K = 2, Nr
1369    
1370           call FIND_RHO(           call FIND_RHO(
1371       I        bi, bj, ibot, itop, jbot, jtop, K, K, eosType,       I        bi, bj, ibot, itop, jbot, jtop, K, K,
1372       I        theta, salt,       I        theta, salt,
1373       O        RHOK,       O        RHOK,
1374       I        myThid )       I        myThid )
1375    
1376           call FIND_RHO(           call FIND_RHO(
1377       I        bi, bj, ibot, itop, jbot, jtop, K-1, K, eosType,       I        bi, bj, ibot, itop, jbot, jtop, K-1, K,
1378       I        theta, salt,       I        theta, salt,
1379       O        RHOKM1,       O        RHOKM1,
1380       I        myThid )       I        myThid )
1381    
1382           call FIND_RHO(           call FIND_RHO(
1383       I        bi, bj, ibot, itop, jbot, jtop, 1, K, eosType,       I        bi, bj, ibot, itop, jbot, jtop, 1, K,
1384       I        theta, salt,       I        theta, salt,
1385       O        RHO1K,       O        RHO1K,
1386       I        myThid )       I        myThid )
1387    
1388           call FIND_ALPHA(           call FIND_ALPHA(
1389       I        bi, bj, ibot, itop, jbot, jtop, K, K, eosType,       I        bi, bj, ibot, itop, jbot, jtop, K, K,
1390       O        WORK1 )       O        WORK1 )
1391    
1392           call FIND_BETA(           call FIND_BETA(
1393       I        bi, bj, ibot, itop, jbot, jtop, K, K, eosType,       I        bi, bj, ibot, itop, jbot, jtop, K, K,
1394       O        WORK2 )       O        WORK2 )
1395    
1396           DO J = jbot, jtop           DO J = jbot, jtop
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22