- Timestamp:
- 2021-07-15T16:05:59+02:00 (3 years ago)
- Location:
- NEMO/branches/UKMO/tools_r4.0-HEAD_dev_MEs/DOMAINcfg/src
- Files:
-
- 2 deleted
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/tools_r4.0-HEAD_dev_MEs/DOMAINcfg/src/zgrmes.f90
r15121 r15124 49 49 ! each vertical sub-zone 50 50 !! * Substitutions 51 #include "domzgr_substitute.h90"52 #include "vectopt_loop_substitute.h90"51 !#include "domzgr_substitute.h90" 52 !#include "vectopt_loop_substitute.h90" 53 53 54 54 CONTAINS … … 906 906 CALL lbc_lnk( e3vw_0, 'V', 1._wp ) 907 907 908 fsdepw(:,:,:) = gdepw_0 (:,:,:)909 fsde3w(:,:,:) = gde3w_0(:,:,:)910 !911 908 where (e3t_0 (:,:,:).eq.0.0) e3t_0(:,:,:) = 1.0 912 909 where (e3u_0 (:,:,:).eq.0.0) e3u_0(:,:,:) = 1.0 … … 917 914 where (e3vw_0 (:,:,:).eq.0.0) e3vw_0(:,:,:) = 1.0 918 915 919 #if defined key_agrif920 ! Ensure meaningful vertical scale factors in ghost lines/columns921 IF( .NOT. Agrif_Root() ) THEN922 !923 IF((nbondi == -1).OR.(nbondi == 2)) THEN924 e3u_0(1,:,:) = e3u_0(2,:,:)925 ENDIF926 !927 IF((nbondi == 1).OR.(nbondi == 2)) THEN928 e3u_0(nlci-1,:,:) = e3u_0(nlci-2,:,:)929 ENDIF930 !931 IF((nbondj == -1).OR.(nbondj == 2)) THEN932 e3v_0(:,1,:) = e3v_0(:,2,:)933 ENDIF934 !935 IF((nbondj == 1).OR.(nbondj == 2)) THEN936 e3v_0(:,nlcj-1,:) = e3v_0(:,nlcj-2,:)937 ENDIF938 !939 ENDIF940 #endif941 942 fsdept(:,:,:) = gdept_0 (:,:,:)943 fsdepw(:,:,:) = gdepw_0 (:,:,:)944 fsde3w(:,:,:) = gde3w_0(:,:,:)945 fse3t (:,:,:) = e3t_0 (:,:,:)946 fse3u (:,:,:) = e3u_0 (:,:,:)947 fse3v (:,:,:) = e3v_0 (:,:,:)948 fse3f (:,:,:) = e3f_0 (:,:,:)949 fse3w (:,:,:) = e3w_0 (:,:,:)950 fse3uw(:,:,:) = e3uw_0 (:,:,:)951 fse3vw(:,:,:) = e3vw_0 (:,:,:)952 953 916 IF( lwp ) WRITE(numout,*) 'Refine mbathy' 954 917 DO jj = 1, jpj 955 918 DO ji = 1, jpi 956 919 DO jk = 1, jpkm1 957 IF( scobot(ji,jj) >= fsdept(ji,jj,jk) ) mbathy(ji,jj) = MAX( 2, jk )920 IF( scobot(ji,jj) >= gdept_0(ji,jj,jk) ) mbathy(ji,jj) = MAX( 2, jk ) 958 921 END DO 959 922 IF( scobot(ji,jj) == 0._wp ) mbathy(ji,jj) = 0 … … 991 954 DO jk = 1, mbathy(ji,jj) 992 955 ! check coordinate is monotonically increasing 993 IF ( fse3w(ji,jj,jk) <= 0._wp .OR. fse3t(ji,jj,jk) <= 0._wp ) THEN956 IF (e3w_0(ji,jj,jk) <= 0._wp .OR. e3t_0(ji,jj,jk) <= 0._wp ) THEN 994 957 WRITE(ctmp1,*) 'ERROR zgr_mes : e3w or e3t =< 0 at point (i,j,k)= ', ji, jj, jk 995 958 WRITE(numout,*) 'ERROR zgr_mes : e3w or e3t =< 0 at point (i,j,k)= ', ji, jj, jk 996 WRITE(numout,*) 'e3w', fse3w(ji,jj,:)997 WRITE(numout,*) 'e3t', fse3t(ji,jj,:)959 WRITE(numout,*) 'e3w',e3w_0(ji,jj,:) 960 WRITE(numout,*) 'e3t',e3t_0(ji,jj,:) 998 961 CALL ctl_stop( ctmp1 ) 999 962 ENDIF 1000 963 ! and check it has never gone negative 1001 IF ( fsdepw(ji,jj,jk) < 0._wp .OR. fsdept(ji,jj,jk) < 0._wp ) THEN964 IF ( gdepw_0(ji,jj,jk) < 0._wp .OR. gdept_0(ji,jj,jk) < 0._wp ) THEN 1002 965 WRITE(ctmp1,*) 'ERROR zgr_mes : gdepw or gdept =< 0 at point (i,j,k)= ', ji, jj, jk 1003 966 WRITE(numout,*) 'ERROR zgr_mes : gdepw or gdept =< 0 at point (i,j,k)= ', ji, jj, jk 1004 WRITE(numout,*) 'gdepw', fsdepw(ji,jj,:)1005 WRITE(numout,*) 'gdept', fsdept(ji,jj,:)967 WRITE(numout,*) 'gdepw',gdepw_0(ji,jj,:) 968 WRITE(numout,*) 'gdept',gdept_0(ji,jj,:) 1006 969 CALL ctl_stop( ctmp1 ) 1007 970 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.