/[MITgcm]/MITgcm/pkg/autodiff/autodiff_store.F
ViewVC logotype

Diff of /MITgcm/pkg/autodiff/autodiff_store.F

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

revision 1.7 by heimbach, Thu Feb 12 23:28:20 2009 UTC revision 1.8 by heimbach, Tue Jun 2 16:10:31 2009 UTC
# Line 547  C-      2D arrays Line 547  C-      2D arrays
547            StoreSEAICE(I,J,bi,bj,19) = dwatn(I,J,bi,bj)            StoreSEAICE(I,J,bi,bj,19) = dwatn(I,J,bi,bj)
548            StoreSEAICE(I,J,bi,bj,20) = seaicemasku(I,J,bi,bj)            StoreSEAICE(I,J,bi,bj,20) = seaicemasku(I,J,bi,bj)
549            StoreSEAICE(I,J,bi,bj,21) = seaicemaskv(I,J,bi,bj)            StoreSEAICE(I,J,bi,bj,21) = seaicemaskv(I,J,bi,bj)
550              StoreSEAICE(I,J,bi,bj,22) = stressDivergenceX(I,J,bi,bj)
551              StoreSEAICE(I,J,bi,bj,23) = stressDivergenceY(I,J,bi,bj)
552  #  else  #  else
553            StoreSEAICE(I,J,bi,bj,19) = 0.0                  StoreSEAICE(I,J,bi,bj,19) = 0.0      
554            StoreSEAICE(I,J,bi,bj,20) = 0.0                      StoreSEAICE(I,J,bi,bj,20) = 0.0          
555            StoreSEAICE(I,J,bi,bj,21) = 0.0            StoreSEAICE(I,J,bi,bj,21) = 0.0
556              StoreSEAICE(I,J,bi,bj,22) = 0.0          
557              StoreSEAICE(I,J,bi,bj,23) = 0.0
558  #  endif /* SEAICE_CGRID */  #  endif /* SEAICE_CGRID */
559  #  ifdef SEAICE_ALLOW_EVP  #  ifdef SEAICE_ALLOW_EVP
560            StoreSEAICE(I,J,bi,bj,22) = seaice_sigma1(I,J,bi,bj)            StoreSEAICE(I,J,bi,bj,24) = seaice_sigma1(I,J,bi,bj)
561            StoreSEAICE(I,J,bi,bj,23) = seaice_sigma2(I,J,bi,bj)            StoreSEAICE(I,J,bi,bj,25) = seaice_sigma2(I,J,bi,bj)
562            StoreSEAICE(I,J,bi,bj,24) = seaice_sigma12(I,J,bi,bj)            StoreSEAICE(I,J,bi,bj,26) = seaice_sigma12(I,J,bi,bj)
563  #  else  #  else
564            StoreSEAICE(I,J,bi,bj,22) = 0.0                  StoreSEAICE(I,J,bi,bj,24) = 0.0      
565            StoreSEAICE(I,J,bi,bj,23) = 0.0                      StoreSEAICE(I,J,bi,bj,25) = 0.0          
566            StoreSEAICE(I,J,bi,bj,24) = 0.0            StoreSEAICE(I,J,bi,bj,26) = 0.0
567  #  endif /* SEAICE_ALLOW_EVP */  #  endif /* SEAICE_ALLOW_EVP */
568  # else  # else
569            StoreSEAICE(I,J,bi,bj,11) = 0.0                  StoreSEAICE(I,J,bi,bj,11) = 0.0      
# Line 576  C-      2D arrays Line 580  C-      2D arrays
580            StoreSEAICE(I,J,bi,bj,22) = 0.0                  StoreSEAICE(I,J,bi,bj,22) = 0.0      
581            StoreSEAICE(I,J,bi,bj,23) = 0.0                      StoreSEAICE(I,J,bi,bj,23) = 0.0          
582            StoreSEAICE(I,J,bi,bj,24) = 0.0            StoreSEAICE(I,J,bi,bj,24) = 0.0
583              StoreSEAICE(I,J,bi,bj,25) = 0.0
584              StoreSEAICE(I,J,bi,bj,26) = 0.0
585  # endif /* SEAICE_ALLOW_DYNAMICS */  # endif /* SEAICE_ALLOW_DYNAMICS */
586  # ifdef SEAICE_SALINITY  # ifdef SEAICE_SALINITY
587            StoreSEAICE(I,J,bi,bj,25) = HSALT(I,J,bi,bj)            StoreSEAICE(I,J,bi,bj,27) = HSALT(I,J,bi,bj)
588  # else  # else
589            StoreSEAICE(I,J,bi,bj,25) = 0.0            StoreSEAICE(I,J,bi,bj,27) = 0.0
590  # endif  # endif
591           enddo           enddo
592          enddo          enddo
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22