Changeset 7154


Ignore:
Timestamp:
2016-10-28T11:16:33+02:00 (4 years ago)
Author:
isabella
Message:

storing the incremental pressure correction in a different variable so that the output will be easier

Location:
branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC/ASM
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90

    r7125 r7154  
    148148      & tbias_p_out,   &   !: Output temperature bias field for P correction 
    149149      & sbias_p_out,   &   !: Output salinity bias field for P correction 
     150      & tbias_i_out,   &   !: Output temperature bias field for incremental P correction 
     151      & sbias_i_out,   &   !: Output salinity bias field for incremental P correction 
    150152      & tbias_asm_stscale, &   !: Short time scale temperature bias field 
    151153      & sbias_asm_stscale      !: Short time scale salinity bias field 
     
    311313            &      tbias_p(jpi,jpj,jpk), & 
    312314            &      sbias_p(jpi,jpj,jpk), & 
     315            &      tbias_i(jpi,jpj,jpk), & 
     316            &      sbias_i(jpi,jpj,jpk), & 
    313317            &      rhd_pc(jpi,jpj,jpk) , & 
    314318            &      gru_pc(jpi,jpj)     , & 
     
    322326                                      sbias_asm_out(jpi,jpj,jpk), &   
    323327                                      tbias_p_out(jpi,jpj,jpk),   &               
    324                                       sbias_p_out(jpi,jpj,jpk)    ) 
     328                                      sbias_p_out(jpi,jpj,jpk),   & 
     329                                      tbias_i_out(jpi,jpj,jpk),   &   
     330                                      sbias_i_out(jpi,jpj,jpk)   ) 
    325331 
    326332         IF( ln_bias_rlx ) ALLOCATE(  tbias_rlx(jpi,jpj,jpk),     & 
     
    370376         tbias_p(:,:,:) = 0.0_wp 
    371377         sbias_p(:,:,:) = 0.0_wp 
     378         tbias_i(:,:,:) = 0.0_wp 
     379         sbias_i(:,:,:) = 0.0_wp 
    372380         gru_pc(:,:)    = 0.0_wp 
    373381         grv_pc(:,:)    = 0.0_wp 
     
    579587      tbias_p(:,:,:) = 0.0_wp 
    580588      sbias_p(:,:,:) = 0.0_wp 
     589      tbias_i(:,:,:) = 0.0_wp 
     590      sbias_i(:,:,:) = 0.0_wp 
    581591 
    582592      ztscale = 0.1_wp 
     
    633643 
    634644                     DO jk = 1, jpkm1 
    635                         tbias_p(:,:,jk) = tbias_p(:,:,jk) + ( t_bkginc(:,:,jk) * zwgt * ztfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
     645                        tbias_i(:,:,jk) = ( t_bkginc(:,:,jk) * zwgt * ztfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
    636646                        &                - ( tbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 
    637                         sbias_p(:,:,jk) = sbias_p(:,:,jk) + ( s_bkginc(:,:,jk) * zwgt * ztfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
     647                        sbias_i(:,:,jk) = ( s_bkginc(:,:,jk) * zwgt * ztfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
    638648                        &                - ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 
    639649                     ENDDO 
     
    684694                    
    685695                   DO jk = 1, jpkm1 
    686                       tbias_p(:,:,jk) = tbias_p(:,:,jk) + ( t_bkginc(:,:,jk) * zwgt * zfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
     696                      tbias_i(:,:,jk) = ( t_bkginc(:,:,jk) * zwgt * zfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
    687697                      &                - ( tbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 
    688                       sbias_p(:,:,jk) = sbias_p(:,:,jk) + ( s_bkginc(:,:,jk) * zwgt * zfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
     698                      sbias_i(:,:,jk) = ( s_bkginc(:,:,jk) * zwgt * zfrac * (1.0 - fbcoef_stscale(:,:)) )         & 
    689699                      &                - ( sbias_asm_stscale(:,:,jk) * zfrac1 * (1.0 - fbcoef_stscale(:,:)) ) 
    690700                   ENDDO 
     
    819829         tbias_p_out(:,:,:)   = tbias_p(:,:,:) 
    820830         sbias_p_out(:,:,:)   = sbias_p(:,:,:) 
     831         tbias_i_out(:,:,:)   = tbias_i(:,:,:) 
     832         sbias_i_out(:,:,:)   = sbias_i(:,:,:) 
    821833      ENDIF 
    822834 
     
    860872         tsw(:,:,:,jp_tem) = tsa(:,:,:,jp_tem) - tbias_p(:,:,:) 
    861873         tsw(:,:,:,jp_sal) = tsa(:,:,:,jp_sal) - sbias_p(:,:,:) 
     874         IF ( ln_incpc ) THEN 
     875            tsw(:,:,:,jp_sal) = tsa(:,:,:,jp_sal) - sbias_p(:,:,:) - sbias_i(:,:,:) 
     876            tsw(:,:,:,jp_tem) = tsa(:,:,:,jp_tem) - tbias_p(:,:,:) - tbias_i(:,:,:) 
     877         ENDIF 
    862878      ELSE 
    863879         tsw(:,:,:,jp_tem) = tsb(:,:,:,jp_tem) - tbias_p(:,:,:) 
    864880         tsw(:,:,:,jp_sal) = tsb(:,:,:,jp_sal) - sbias_p(:,:,:) 
     881            IF ( ln_incpc ) THEN 
     882               tsw(:,:,:,jp_sal) = tsa(:,:,:,jp_sal) - sbias_p(:,:,:) - sbias_i(:,:,:) 
     883               tsw(:,:,:,jp_tem) = tsa(:,:,:,jp_tem) - tbias_p(:,:,:) - tbias_i(:,:,:) 
     884            ENDIF 
    865885      ENDIF 
    866886 
     
    940960         CALL iom_rstput( nn_bias_itwrt, nn_bias_itwrt, numbias_tot, 'sbias_asm' , sbias_asm_out )    
    941961         CALL iom_rstput( nn_bias_itwrt, nn_bias_itwrt, numbias_tot, 'tbias_p'   , tbias_p_out ) 
    942          CALL iom_rstput( nn_bias_itwrt, nn_bias_itwrt, numbias_tot, 'sbias_p'   , sbias_p_out )          
     962         CALL iom_rstput( nn_bias_itwrt, nn_bias_itwrt, numbias_tot, 'sbias_p'   , sbias_p_out ) 
     963         CALL iom_rstput( nn_bias_itwrt, nn_bias_itwrt, numbias_tot, 'tbias_i'   , tbias_i_out ) 
     964         CALL iom_rstput( nn_bias_itwrt, nn_bias_itwrt, numbias_tot, 'sbias_i'   , sbias_i_out )           
    943965      ENDIF 
    944966 
  • branches/UKMO/dev_r5518_pcbias_ipc/NEMOGCM/NEMO/OPA_SRC/ASM/biaspar.F90

    r7128 r7154  
    2727      & tbias, &       !: Temperature bias field for T correction 
    2828      & tbias_p, &     !:  "           "    "    "   P correction 
     29      & tbias_i, &     !:  "           "    "    "   incremental P correction 
    2930      & sbias, &       !: Salinity bias field    for S correction 
    3031      & sbias_p, &     !:  "           "    "    "   P correction 
     32      & sbias_i, &     !:  "           "    "    "   incremental P correction 
    3133      & rhd_pc         !: Press corrtd density from online to use in dyn_hpg 
    3234 
Note: See TracChangeset for help on using the changeset viewer.