/[MITgcm]/MITgcm_contrib/torge/itd/code/seaice_growth.F
ViewVC logotype

Diff of /MITgcm_contrib/torge/itd/code/seaice_growth.F

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

revision 1.3 by torge, Wed Sep 26 17:50:17 2012 UTC revision 1.4 by torge, Fri Sep 28 19:01:17 2012 UTC
# Line 1163  C     Negative sublimation is resublimat Line 1163  C     Negative sublimation is resublimat
1163          ENDDO          ENDDO
1164          DO J=1,sNy          DO J=1,sNy
1165           DO I=1,sNx           DO I=1,sNx
1166              a_QbyATM_open(I,J)    = a_QbyATM_open(I,J)
1167         &         * convertQ2HI * ( ONE - AREApreTH(I,J) )
1168              a_QSWbyATM_open(I,J)  = a_QSWbyATM_open(I,J)
1169         &         * convertQ2HI * ( ONE - AREApreTH(I,J) )
1170  C and initialize r_QbyATM_open  C and initialize r_QbyATM_open
1171            r_QbyATM_open(I,J)=a_QbyATM_open(I,J)            r_QbyATM_open(I,J)=a_QbyATM_open(I,J)
1172           ENDDO           ENDDO
# Line 1272  c available turbulent flux Line 1276  c available turbulent flux
1276            a_QbyOCN(i,j) =            a_QbyOCN(i,j) =
1277       &         tmpscal1 * tmpscal2 * MixedLayerTurbulenceFactor       &         tmpscal1 * tmpscal2 * MixedLayerTurbulenceFactor
1278            r_QbyOCN(i,j) = a_QbyOCN(i,j)            r_QbyOCN(i,j) = a_QbyOCN(i,j)
 ctm  
           if (i.eq.20 .and. j.eq.20) then  
             print *, HeatCapacity_Cp  
             print *, rhoConst  
             print *, recip_QI  
             print *, theta(20,20,kSurface,bi,bj)  
             print *, tempFrz  
             print *, SEAICE_deltaTtherm  
             print *, maskC(20,20,kSurface,bi,bj)  
             print *, tmpscal2  
             print *, a_QbyOCN(20,20)  
           endif  
 ctm  
1279           ENDDO           ENDDO
1280          ENDDO          ENDDO
1281    
# Line 1378  C       then update totals Line 1369  C       then update totals
1369          ENDDO          ENDDO
1370  c ToM<<< debug seaice_growth  c ToM<<< debug seaice_growth
1371          WRITE(msgBuf,'(A,7F6.2)')          WRITE(msgBuf,'(A,7F6.2)')
1372    #ifdef SEAICE_ITD
1373       &    ' SEAICE_GROWTH: Heff increments 1, HEFFITD = ',       &    ' SEAICE_GROWTH: Heff increments 1, HEFFITD = ',
1374       &     HEFFITD(20,20,:,bi,bj)       &     HEFFITD(20,20,:,bi,bj)
1375    #else
1376         &    ' SEAICE_GROWTH: Heff increments 1, HEFF = ',
1377         &     HEFF(20,20,bi,bj)
1378    #endif
1379          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
1380       &    SQUEEZE_RIGHT , myThid)       &    SQUEEZE_RIGHT , myThid)
1381  c ToM>>>  c ToM>>>
# Line 1410  C          and hence weighted by fractio Line 1406  C          and hence weighted by fractio
1406            ENDDO            ENDDO
1407           ENDDO           ENDDO
1408          ENDDO          ENDDO
 c ToM<<< debug seaice_growth  
         WRITE(msgBuf,'(A,7F9.6)')  
      &    ' SEAICE_GROWTH: d_HEFFbyOCNonICE w/ITD: ',  
      &     d_HEFFbyOCNonICE(20,20)  
         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
      &    SQUEEZE_RIGHT , myThid)  
 c ToM>>>  
1409  #else /* SEAICE_ITD */  #else /* SEAICE_ITD */
1410          DO J=1,sNy          DO J=1,sNy
1411           DO I=1,sNx           DO I=1,sNx
# Line 1428  c ToM>>> Line 1417  c ToM>>>
1417  #endif  #endif
1418           ENDDO           ENDDO
1419          ENDDO          ENDDO
 c ToM<<< debug seaice_growth  
         WRITE(msgBuf,'(A,7F9.6)')  
      &    ' SEAICE_GROWTH: d_HEFFbyOCNonICE w/o ITD: ',  
      &     d_HEFFbyOCNonICE(20,20)  
         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
      &    SQUEEZE_RIGHT , myThid)  
 c ToM>>>  
1420  #endif /* SEAICE_ITD */  #endif /* SEAICE_ITD */
1421  c ToM<<< debug seaice_growth  c ToM<<< debug seaice_growth
1422          WRITE(msgBuf,'(A,7F6.2)')          WRITE(msgBuf,'(A,7F6.2)')
1423    #ifdef SEAICE_ITD
1424       &    ' SEAICE_GROWTH: Heff increments 2, HEFFITD = ',       &    ' SEAICE_GROWTH: Heff increments 2, HEFFITD = ',
1425       &     HEFFITD(20,20,:,bi,bj)       &     HEFFITD(20,20,:,bi,bj)
1426    #else
1427         &    ' SEAICE_GROWTH: Heff increments 2, HEFF = ',
1428         &     HEFF(20,20,bi,bj)
1429    #endif
1430          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
1431       &    SQUEEZE_RIGHT , myThid)       &    SQUEEZE_RIGHT , myThid)
1432  c ToM>>>  c ToM>>>
# Line 1492  Cgf no additional dependency through sno Line 1479  Cgf no additional dependency through sno
1479  #endif /* SEAICE_ITD */  #endif /* SEAICE_ITD */
1480  c ToM<<< debug seaice_growth  c ToM<<< debug seaice_growth
1481          WRITE(msgBuf,'(A,7F6.2)')          WRITE(msgBuf,'(A,7F6.2)')
1482    #ifdef SEAICE_ITD
1483       &    ' SEAICE_GROWTH: Heff increments 3, HEFFITD = ',       &    ' SEAICE_GROWTH: Heff increments 3, HEFFITD = ',
1484       &     HEFFITD(20,20,:,bi,bj)       &     HEFFITD(20,20,:,bi,bj)
1485    #else
1486         &    ' SEAICE_GROWTH: Heff increments 3, HEFF = ',
1487         &     HEFF(20,20,bi,bj)
1488    #endif
1489          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
1490       &    SQUEEZE_RIGHT , myThid)       &    SQUEEZE_RIGHT , myThid)
1491  c ToM>>>  c ToM>>>
# Line 1562  c         Limit ice growth by potential Line 1554  c         Limit ice growth by potential
1554  #endif /* SEAICE_ITD */  #endif /* SEAICE_ITD */
1555  c ToM<<< debug seaice_growth  c ToM<<< debug seaice_growth
1556          WRITE(msgBuf,'(A,7F6.2)')          WRITE(msgBuf,'(A,7F6.2)')
1557    #ifdef SEAICE_ITD
1558       &    ' SEAICE_GROWTH: Heff increments 4, HEFFITD = ',       &    ' SEAICE_GROWTH: Heff increments 4, HEFFITD = ',
1559       &     HEFFITD(20,20,:,bi,bj)       &     HEFFITD(20,20,:,bi,bj)
1560    #else
1561         &    ' SEAICE_GROWTH: Heff increments 4, HEFF = ',
1562         &     HEFF(20,20,bi,bj)
1563    #endif
1564          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
1565       &    SQUEEZE_RIGHT , myThid)       &    SQUEEZE_RIGHT , myThid)
1566  c ToM>>>  c ToM>>>
# Line 1617  Cgf rain to snow is not a surface proces Line 1614  Cgf rain to snow is not a surface proces
1614  #endif /* ALLOW_ATM_TEMP */  #endif /* ALLOW_ATM_TEMP */
1615  c ToM<<< debug seaice_growth  c ToM<<< debug seaice_growth
1616          WRITE(msgBuf,'(A,7F6.2)')          WRITE(msgBuf,'(A,7F6.2)')
1617    #ifdef SEAICE_ITD
1618       &    ' SEAICE_GROWTH: Heff increments 5, HEFFITD = ',       &    ' SEAICE_GROWTH: Heff increments 5, HEFFITD = ',
1619       &     HEFFITD(20,20,:,bi,bj)       &     HEFFITD(20,20,:,bi,bj)
1620    #else
1621         &    ' SEAICE_GROWTH: Heff increments 5, HEFF = ',
1622         &     HEFF(20,20,bi,bj)
1623    #endif
1624          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
1625       &    SQUEEZE_RIGHT , myThid)       &    SQUEEZE_RIGHT , myThid)
1626  c ToM>>>  c ToM>>>
# Line 1671  Cgf no additional dependency through sno Line 1673  Cgf no additional dependency through sno
1673  Cph)  Cph)
1674  c ToM<<< debug seaice_growth  c ToM<<< debug seaice_growth
1675          WRITE(msgBuf,'(A,7F6.2)')          WRITE(msgBuf,'(A,7F6.2)')
1676    #ifdef SEAICE_ITD
1677       &    ' SEAICE_GROWTH: Heff increments 6, HEFFITD = ',       &    ' SEAICE_GROWTH: Heff increments 6, HEFFITD = ',
1678       &     HEFFITD(20,20,:,bi,bj)       &     HEFFITD(20,20,:,bi,bj)
1679    #else
1680         &    ' SEAICE_GROWTH: Heff increments 6, HEFF = ',
1681         &     HEFF(20,20,bi,bj)
1682    #endif
1683          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
1684       &    SQUEEZE_RIGHT , myThid)       &    SQUEEZE_RIGHT , myThid)
1685  c ToM>>>  c ToM>>>
# Line 1752  c needs to be here to allow use also wit Line 1759  c needs to be here to allow use also wit
1759  #endif /* ALLOW_SITRACER */  #endif /* ALLOW_SITRACER */
1760  c ToM<<< debug seaice_growth  c ToM<<< debug seaice_growth
1761          WRITE(msgBuf,'(A,7F6.2)')          WRITE(msgBuf,'(A,7F6.2)')
1762    #ifdef SEAICE_ITD
1763       &    ' SEAICE_GROWTH: Heff increments 7, HEFFITD = ',       &    ' SEAICE_GROWTH: Heff increments 7, HEFFITD = ',
1764       &     HEFFITD(20,20,:,bi,bj)       &     HEFFITD(20,20,:,bi,bj)
1765    #else
1766         &    ' SEAICE_GROWTH: Heff increments 7, HEFF = ',
1767         &     HEFF(20,20,bi,bj)
1768    #endif
1769          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
1770       &    SQUEEZE_RIGHT , myThid)       &    SQUEEZE_RIGHT , myThid)
1771  c ToM>>>  c ToM>>>
# Line 1799  CADJ STORE hsnow(:,:,bi,bj) = comlev1_bi Line 1811  CADJ STORE hsnow(:,:,bi,bj) = comlev1_bi
1811  #endif /* SEAICE_GROWTH_LEGACY */  #endif /* SEAICE_GROWTH_LEGACY */
1812  c ToM<<< debug seaice_growth  c ToM<<< debug seaice_growth
1813          WRITE(msgBuf,'(A,7F6.2)')          WRITE(msgBuf,'(A,7F6.2)')
1814    #ifdef SEAICE_ITD
1815       &    ' SEAICE_GROWTH: Heff increments 8, HEFFITD = ',       &    ' SEAICE_GROWTH: Heff increments 8, HEFFITD = ',
1816       &     HEFFITD(20,20,:,bi,bj)       &     HEFFITD(20,20,:,bi,bj)
1817    #else
1818         &    ' SEAICE_GROWTH: Heff increments 8, HEFF = ',
1819         &     HEFF(20,20,bi,bj)
1820    #endif
1821          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
1822       &    SQUEEZE_RIGHT , myThid)       &    SQUEEZE_RIGHT , myThid)
1823  c ToM>>>  c ToM>>>
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22