Changeset 8649 for branches/NERC
- Timestamp:
- 2017-10-20T16:46:23+02:00 (7 years ago)
- Location:
- branches/NERC/dev_r5518_GO6_Carb_Debug/NEMOGCM/NEMO/TOP_SRC/MEDUSA
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_GO6_Carb_Debug/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90
r8521 r8649 33 33 !!---------------------------------------------------------------------- 34 34 USE bio_medusa_mod 35 USE dom_oce, ONLY: atfp, atfp1, neuler, rdt, e3t_n, tmask 35 USE dom_oce, ONLY: atfp, atfp1, neuler, rdt, & 36 e3t_0, e3t_n, tmask 36 37 USE in_out_manager, ONLY: lwp, numout 37 38 # if defined key_iomput … … 55 56 USE trcnam_trp, ONLY: ln_trcadv_cen2, ln_trcadv_tvd 56 57 58 !!* Substitution 59 # include "domzgr_substitute.h90" 60 57 61 !! time (integer timestep) 58 62 INTEGER, INTENT( in ) :: kt … … 61 65 INTEGER :: jn 62 66 63 REAL(wp) :: fq0,fq1,fq2,fq3 67 REAL(wp) :: fq0,fq1,fq2,fq3,fq4 64 68 65 69 # if defined key_roam … … 148 152 fq2 = fq0 + fq1 149 153 fq3 = f_benout_n(ji,jj) 150 if (lwp) write (numout,' a,2i3,a,4f15,5)')&154 if (lwp) write (numout,'(a,2i3,a,4f15.5)') & 151 155 'AXY N cons: (i,j)=',ji,jj,', (flx,ben,err,out)=', & 152 156 fq0,fq1,fq2,fq3 … … 162 166 fq2 = fq0 + fq1 163 167 fq3 = f_benout_si(ji,jj) 164 if (lwp) write (numout,' a,2i3,a,4f15,5)')&168 if (lwp) write (numout,'(a,2i3,a,4f15.5)') & 165 169 'AXY Si cons: (i,j)=',ji,jj,', (flx,ben,err,out)=', & 166 170 fq0,fq1,fq2,fq3 … … 177 181 fq3 = fq0 + fq1 178 182 fq4 = f_benout_c(ji,jj) + f_benout_ca(ji,jj) 179 if (lwp) write (numout,' a,2i3,a,5f15,5)')&183 if (lwp) write (numout,'(a,2i3,a,5f15.5)') & 180 184 'AXY C cons: (i,j)=',ji,jj,', (flx,ben,asf,err,out)=', & 181 185 fq0,fq1,fq2,fq3,fq4 … … 191 195 fq2 = fq0 + fq1 192 196 fq3 = 2.0 * f_benout_ca(ji,jj) 193 if (lwp) write (numout,' a,2i3,a,4f15,5)')&197 if (lwp) write (numout,'(a,2i3,a,4f15.5)') & 194 198 'AXY alk cons: (i,j)=',ji,jj,', (flx,ben,err,out)=', & 195 199 fq0,fq1,fq2,fq3 -
branches/NERC/dev_r5518_GO6_Carb_Debug/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90
r8643 r8649 181 181 REAL(wp) :: sumtsn, tsnavg 182 182 INTEGER :: summask 183 CHARACTER(25) :: charout, charout2, charout3 183 184 !! 184 185 !!------------------------------------------------------------------ … … 477 478 !! AXY (28/02/14): check input fields 478 479 if (ztmp(ji,jj) .lt. -3.0 .or. ztmp(ji,jj) .gt. 40.0 ) then 480 write(charout,*) tsn(ji,jj,jk,jp_tem) 481 write(charout2,*) mig(ji), mjg(jj), jk, kt 479 482 Call ctl_warn(' trc_bio_medusa: T WARNING 3D : ', & 480 tsn(ji,jj,jk,jp_tem), 'at I, J, K, kt :',&481 mig(ji), mjg(jj), jk, kt)483 TRIM(charout), 'at I, J, K, kt :', & 484 TRIM(charout2) ) 482 485 # if defined key_debug_medusa 483 486 !! temperature … … 543 546 ( tsnavg .GT. 40.0 ) ) THEN 544 547 IF (ztmp(ji,jj) .LT. -3.0 ) THEN 545 CALL ctl_warn( &546 ' trc_bio_medusa: T SWITCHING 3D, ',&547 tsn(ji,jj,jk,jp_tem), ' -> -3.0 ')548 write(charout,*) tsn(ji,jj,jk,jp_tem) 549 Call ctl_warn(' trc_bio_medusa: T SWITCHING 3D : ', & 550 TRIM(charout), ' -> -3.0 ') 548 551 ztmp(ji,jj) = -3.0 549 552 ENDIF 550 553 IF (ztmp(ji,jj) .GT. 40.0 ) THEN 551 CALL ctl_warn( &552 ' trc_bio_medusa: T SWITCHING 3D, ',&553 tsn(ji,jj,jk,jp_tem), ' -> 40.0 ')554 write(charout,*) tsn(ji,jj,jk,jp_tem) 555 Call ctl_warn(' trc_bio_medusa: T SWITCHING 3D : ', & 556 TRIM(charout), ' -> 40.0 ') 554 557 ztmp(ji,jj) = 40.0 555 558 ENDIF 556 559 ELSE 557 CALL ctl_warn( & 558 ' trc_bio_medusa: T SWITCHING 3D, ', & 559 tsn(ji,jj,jk,jp_tem), ' -> surounding avg : ', tsnavg ) 560 write(charout,*) tsn(ji,jj,jk,jp_tem) 561 write(charout2,*) tsnavg 562 Call ctl_warn(' trc_bio_medusa: T SWITCHING 3D : ', & 563 TRIM(charout), ' -> surounding avg : ', & 564 TRIM(charout2) ) 560 565 ztmp(ji,jj) = tsnavg 561 566 ENDIF … … 563 568 !! end T chack 564 569 if (zsal(ji,jj) .lt. 1.0 .or. zsal(ji,jj) .gt. 47.0 ) then 565 CALL ctl_warn( & 566 ' trc_bio_medusa: S WARNING 2D, ', & 567 tsn(ji,jj,jk,jp_sal), & 568 ' at I, J, K, kt : ', mig(ji), mjg(jj), jk, kt) 570 write(charout,*) tsn(ji,jj,jk,jp_sal) 571 write(charout2,*) mig(ji), mjg(jj), jk, kt 572 Call ctl_warn(' trc_bio_medusa: S WARNING 3D : ', & 573 TRIM(charout), 'at I, J, K, kt :', & 574 TRIM(charout2) ) 569 575 !! Correct out of range values 570 576 IF (zsal(ji,jj) .LT. 1.0 ) THEN 571 CALL ctl_warn( &572 ' trc_bio_medusa: S SWITCHING 3D, ',&573 tsn(ji,jj,jk,jp_sal), ' -> 1.0 ')577 write(charout,*) tsn(ji,jj,jk,jp_sal) 578 Call ctl_warn(' trc_bio_medusa: S SWITCHING 3D : ', & 579 TRIM(charout), ' -> 1.0 ') 574 580 zsal(ji,jj) = 1.0 575 581 ENDIF 576 582 IF (zsal(ji,jj) .GT. 47.0 ) THEN 577 CALL ctl_warn( &578 ' trc_bio_medusa: T SWITCHING 3D, ',&579 tsn(ji,jj,jk,jp_sal), ' -> 47.0 ')583 write(charout,*) tsn(ji,jj,jk,jp_sal) 584 Call ctl_warn(' trc_bio_medusa: S SWITCHING 3D : ', & 585 TRIM(charout), ' -> 47.0 ') 580 586 zsal(ji,jj) = 47.0 581 587 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.