Changeset 7652

Timestamp:

2017-02-06T17:14:30+01:00 (8 years ago)

Author:

isabella

Message:

adding output messages to understand what model does

Location:

branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC

Files:

• ASM/bias.F90 (modified) (1 diff)
• DYN/dynhpg.F90 (modified) (2 diffs)

branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90

@@ -317 +317 @@
          IF ( (.NOT. ln_incpc) .AND. ln_incpc_only) &   
             &   CALL ctl_stop (' if you set ln_incpc_only to .true. then you need to set ln_incpc to .true. as well' )
-
+         
+         WRITE(numout,*) '     time step is    = ',rdt,'you choose to write pcbias at nn_bias_itwrt  = ',nn_bias_itwrt,'and end of iau is rday/rdt=',rday/rdt 
       ENDIF
       IF( .NOT. ln_bias ) RETURN

branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90

@@ -99 +99 @@
       !
       IF ( ln_bias .AND. ln_bias_pc_app ) THEN
+         IF(lwp) THEN 
+         WRITE(numout,*) " ! store orig density, use pressure corrected density"
+         ENDIF
          z_rhd_st(:,:,:) = rhd(:,:,:)     ! store orig density 
          rhd(:,:,:)      = rhd_pc(:,:,:)  ! use pressure corrected density
@@ -127 +130 @@
       !
       IF ( ln_bias .AND. ln_bias_pc_app )  THEN
+         IF(lwp) THEN 
+         WRITE(numout,*) " ! restore original density"
+         ENDIF
          rhd(:,:,:) = z_rhd_st(:,:,:)     ! restore original density
          gru(:,:)   = z_gru_st(:,:)