Changeset 13428


Ignore:
Timestamp:
2020-08-21T19:32:10+02:00 (5 months ago)
Author:
smueller
Message:

Adjustments of modifications made in [5168] (see ticket #1360) and bugfix in module asminc

Location:
NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

    r8029 r13428  
    703703        DO jk = 1, jpk 
    704704           DO jj = 1, jpj 
    705               DO ji = 1, jpk 
     705              DO ji = 1, jpi 
    706706                 fzptnz (ji,jj,jk) = ( -0.0575_wp + 1.710523e-3_wp * SQRT( tsn(ji,jj,jk,jp_sal) )                   & 
    707707                                                  - 2.154996e-4_wp *       tsn(ji,jj,jk,jp_sal)   ) * tsn(ji,jj,jk,jp_sal)  & 
     
    735735                  ! Do not apply negative increments if the temperature will fall below freezing 
    736736                  WHERE(t_bkginc(:,:,jk) > 0.0_wp .OR. & 
    737                      &   t_bkg(:,:,jk) + tsa(:,:,jk,jp_tem) + t_bkginc(:,:,jk) * wgtiau(it) > fzptnz(:,:,jk) ) 
     737                     &   tsn(:,:,jk,jp_tem) + tsa(:,:,jk,jp_tem) + t_bkginc(:,:,jk) * wgtiau(it) > fzptnz(:,:,jk) ) 
    738738                     tsa(:,:,jk,jp_tem) = tsa(:,:,jk,jp_tem) + t_bkginc(:,:,jk) * zincwgt 
    739739                  END WHERE 
     
    745745                  ! minimum value salfixmin 
    746746                  WHERE(s_bkginc(:,:,jk) > 0.0_wp .OR. & 
    747                      &  s_bkg(:,:,jk) + tsa(:,:,jk,jp_sal) + s_bkginc(:,:,jk) * wgtiau(it) > salfixmin ) 
     747                     &  tsn(:,:,jk,jp_sal) + tsa(:,:,jk,jp_sal) + s_bkginc(:,:,jk) * wgtiau(it) > salfixmin ) 
    748748                     tsa(:,:,jk,jp_sal) = tsa(:,:,jk,jp_sal) + s_bkginc(:,:,jk) * zincwgt 
    749749                  END WHERE 
     
    775775               ! Do not apply negative increments if the temperature will fall below freezing 
    776776               WHERE(t_bkginc(:,:,:) > 0.0_wp .OR. & 
    777                   &   tsn(:,:,:,jp_tem) + t_bkginc(:,:,:) > fzptnz(:,:,:) ) 
     777                  &   t_bkg(:,:,:) + t_bkginc(:,:,:) > fzptnz(:,:,:) ) 
    778778                  tsn(:,:,:,jp_tem) = t_bkg(:,:,:) + t_bkginc(:,:,:) 
    779779               END WHERE 
     
    785785               ! minimum value salfixmin 
    786786               WHERE(s_bkginc(:,:,:) > 0.0_wp .OR. & 
    787                   &   tsn(:,:,:,jp_sal) + s_bkginc(:,:,:) > salfixmin ) 
     787                  &   s_bkg(:,:,:) + s_bkginc(:,:,:) > salfixmin ) 
    788788                  tsn(:,:,:,jp_sal) = s_bkg(:,:,:) + s_bkginc(:,:,:) 
    789789               END WHERE 
  • NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90

    r5168 r13428  
    111111      !!              Madec, 2008, internal report, IPSL. 
    112112      !!---------------------------------------------------------------------- 
    113       INTEGER ::   istp, jk      ! time step index 
     113      INTEGER ::   istp       ! time step index 
    114114      !!---------------------------------------------------------------------- 
    115115      ! 
     
    156156            CALL Agrif_Step( stp )           ! AGRIF: time stepping 
    157157#else 
    158             !-------------------------------------------------------------! 
    159             ! This trick ensures a minimum consistency between mixing coef. 
    160             ! and other variables 
    161             IF( lk_zdftke .AND. lk_asminc .AND. ln_asmdin .AND. ( istp == nit000) ) THEN 
    162             CALL tke_avn                        ! recompute avt, avm, 
    163                                                 ! avmu, avmv and dissl (approximation) 
    164                DO jk = nit000, nit000 + 2 
    165                   CALL zdf_tke( jk )   ; 
    166                END DO 
    167  
    168             END IF 
    169             !-------------------------------------------------------------! 
    170  
    171158            CALL stp( istp )                 ! standard time stepping 
    172159#endif 
  • NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/step.F90

    r3604 r13428  
    112112      ! 
    113113      !  VERTICAL PHYSICS 
     114 
     115#if defined key_zdftke 
     116      !-------------------------------------------------------------! 
     117      ! This trick ensures a minimum consistency between mixing coef. 
     118      ! and other variables 
     119      IF( lk_zdftke .AND. lk_asminc .AND. ln_asmdin .AND. ( kstp == nit000) ) THEN 
     120         CALL tke_avn                        ! recompute avt, avm, 
     121                                             ! avmu, avmv and dissl (approximation) 
     122         DO jk = nit000, nit000 + 2 
     123            CALL zdf_tke( jk ) 
     124         END DO 
     125          
     126      END IF 
     127      !-------------------------------------------------------------! 
     128#endif 
     129 
    114130                         CALL zdf_bfr( kstp )         ! bottom friction 
    115131 
Note: See TracChangeset for help on using the changeset viewer.