Changeset 9130 for branches/NERC
 Timestamp:
 20171219T14:15:21+01:00 (7 years ago)
 Location:
 branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_8356/NEMOGCM/NEMO
 Files:

 2 edited
Legend:
 Unmodified
 Added
 Removed

branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_8356/NEMOGCM/NEMO/OPA_SRC/stpctl.F90
r9073 r9130 157 157 ztmin = 100.e0 ! arbitrary high minimum value 158 158 DO jj = 2, jpjm1 159 DO ji = 1, jpi159 DO ji = 2, jpim1 160 160 IF( tmask(ji,jj,1) == 1) THEN 161 161 ztmax = MAX(ztmax,tsn(ji,jj,1,jp_tem)) ! find local maximum 
branches/NERC/dev_r5518_GO6_Carb_Fail_from_GO6_8356/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90
r9073 r9130 821 821 !!========================== 822 822 !! DIC Check 823 zmax = 5. e0 ! arbitrary low maximum value823 zmax = 5.0 ! arbitrary low maximum value 824 824 zmin = 4.0E4 ! arbitrary high minimum value 825 825 DO jj = 2, jpjm1 826 DO ji = 1, jpi826 DO ji = 2,jpim1 827 827 IF( tmask(ji,jj,1) == 1) THEN 828 828 zmax = MAX(zmax,zdic(ji,jj)) ! find local maximum … … 834 834 IF( lk_mpp ) CALL mpp_min( zmin ) ! min over the global domain 835 835 ! 836 IF( zmax >4.0E3) THEN ! we've got a problem836 IF( zmax .GT. 4.0E3) THEN ! we've got a problem 837 837 IF (lk_mpp) THEN 838 838 CALL mpp_maxloc ( zdic(:,:),tmask(:,:,1), zmax, ii,ij ) 839 839 ELSE 840 ilocs = MAXLOC( zdic(:,:), mask = tmask(:,:,1) == 1. e0)840 ilocs = MAXLOC( zdic(:,:), mask = tmask(:,:,1) == 1. ) 841 841 ii = ilocs(1) + nimpp  1 842 842 ij = ilocs(2) + njmpp  1 … … 851 851 ENDIF 852 852 ! 853 IF( zmin <=0.0) THEN ! we've got a problem853 IF( zmin .LE. 0.0) THEN ! we've got a problem 854 854 IF (lk_mpp) THEN 855 855 CALL mpp_minloc ( zdic(:,:),tmask(:,:,1), zmin, ii,ij ) 856 856 ELSE 857 ilocs = MINLOC( zdic(:,:), mask = tmask(:,:,1) == 1. e0)857 ilocs = MINLOC( zdic(:,:), mask = tmask(:,:,1) == 1. ) 858 858 ii = ilocs(1) + nimpp  1 859 859 ij = ilocs(2) + njmpp  1 … … 870 870 !!========================== 871 871 !! ALKALINITY Check 872 zmax = 5. e0 ! arbitrary low maximum value872 zmax = 5.0 ! arbitrary low maximum value 873 873 zmin = 4.0E4 ! arbitrary high minimum value 874 874 DO jj = 2, jpjm1 875 DO ji = 1, jpi875 DO ji = 2,jpim1 876 876 IF( tmask(ji,jj,1) == 1) THEN 877 877 zmax = MAX(zmax,zalk(ji,jj)) ! find local maximum … … 883 883 IF( lk_mpp ) CALL mpp_min( zmin ) ! min over the global domain 884 884 ! 885 IF( zmax >4.0E3) THEN ! we've got a problem885 IF( zmax .GT. 4.0E3) THEN ! we've got a problem 886 886 IF (lk_mpp) THEN 887 887 CALL mpp_maxloc ( zalk(:,:),tmask(:,:,1), zmax, ii,ij ) 888 888 ELSE 889 ilocs = MAXLOC( zalk(:,:), mask = tmask(:,:,1) == 1. e0)889 ilocs = MAXLOC( zalk(:,:), mask = tmask(:,:,1) == 1. ) 890 890 ii = ilocs(1) + nimpp  1 891 891 ij = ilocs(2) + njmpp  1 … … 900 900 ENDIF 901 901 ! 902 IF( zmin <=0.0) THEN ! we've got a problem902 IF( zmin .LE. 0.0) THEN ! we've got a problem 903 903 IF (lk_mpp) THEN 904 904 CALL mpp_minloc ( zalk(:,:),tmask(:,:,1), zmin, ii,ij ) 905 905 ELSE 906 ilocs = MINLOC( zalk(:,:), mask = tmask(:,:,1) == 1. e0)906 ilocs = MINLOC( zalk(:,:), mask = tmask(:,:,1) == 1. ) 907 907 ii = ilocs(1) + nimpp  1 908 908 ij = ilocs(2) + njmpp  1
Note: See TracChangeset
for help on using the changeset viewer.