625 |
C-- end of thermodynamic k loop (Nr:1) |
C-- end of thermodynamic k loop (Nr:1) |
626 |
ENDDO |
ENDDO |
627 |
|
|
|
|
|
|
#ifdef ALLOW_AUTODIFF_TAMC |
|
|
C? Patrick? What about this one? |
|
|
cph Keys iikey and idkey dont seem to be needed |
|
|
cph since storing occurs on different tape for each |
|
|
cph impldiff call anyways. |
|
|
cph Thus, common block comlev1_impl isnt needed either. |
|
|
cph Storing below needed in the case useGMREDI. |
|
|
iikey = (ikey-1)*maximpl |
|
|
#endif /* ALLOW_AUTODIFF_TAMC */ |
|
|
|
|
628 |
C-- Implicit diffusion |
C-- Implicit diffusion |
629 |
IF (implicitDiffusion) THEN |
IF (implicitDiffusion) THEN |
630 |
|
|
631 |
IF (tempStepping) THEN |
IF (tempStepping) THEN |
632 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
|
idkey = iikey + 1 |
|
633 |
CADJ STORE gT(:,:,:,bi,bj) = comlev1_bibj , key=ikey, byte=isbyte |
CADJ STORE gT(:,:,:,bi,bj) = comlev1_bibj , key=ikey, byte=isbyte |
634 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
635 |
CALL IMPLDIFF( |
CALL IMPLDIFF( |
641 |
|
|
642 |
IF (saltStepping) THEN |
IF (saltStepping) THEN |
643 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
|
idkey = iikey + 2 |
|
644 |
CADJ STORE gS(:,:,:,bi,bj) = comlev1_bibj , key=ikey, byte=isbyte |
CADJ STORE gS(:,:,:,bi,bj) = comlev1_bibj , key=ikey, byte=isbyte |
645 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
646 |
CALL IMPLDIFF( |
CALL IMPLDIFF( |