Changeset 6 in lmdz_wrf for WRFV3/dyn_em

Timestamp:

Jul 18, 2014, 11:20:11 AM (10 years ago)

Author:

lfita

Message:

Fixing problem with the solar radiation:

1.- Correct values for the julian and the hour (0,1) to the LMDZ physiq

Adding NaN checking for 'PSFC' and 'T_2'

File:

• WRFV3/dyn_em/solve_em.F (modified) (3 diffs)

WRFV3/dyn_em/solve_em.F

@@ -186 +186 @@
    CHARACTER(LEN=256)                         :: mminlu, mminsl
    CHARACTER(LEN=50)                          :: errmsg
+   INTEGER                                    :: hr, minute, sec, ms, julyr, &
+     julday
+   REAL                                       :: gmt
+
 
    errmsg = 'ERROR -- error -- ERROR -- error'
@@ -845 +849 @@
 
        grid%qv_2 = moist(:,:,:,P_QV)
+
+! L. Fita, LMD. July 2014. Getting hour of the day
+       CALL domain_clock_get( grid, current_time=CurrTime )
+       CALL WRFU_TimeGet( CurrTime, YY= julyr, dayOfYear=julday, H=hr, M=minute, S=sec, MS=ms, rc=rc)
+! Julian day hour (0, 1) !!
+       gmt=(hr+real(minute)/60.+real(sec)/3600.+real(ms)/(1000*3600))/24.
+
+! Checking for NaNs (should not be necessary but....)
+   im2 = ims + (ime - ims) / 2
+   jm2 = jms + (jme - jms) / 2
+
+   IF (grid%t_2(im2,iz,jm2) /= grid%t_2(im2,iz,jm2) .OR. ABS(grid%t_2(im2,iz,jm2)) > 10000. ) THEN
+     PRINT errmsg
+     WRITE(wrF_err_message,*)'solve_em: wrong T value=',                 &
+       grid%t_2(im2,iz,jm2),' at: ', im2,', ', iz,', ', jm2,' !!!'
+#ifdef DM_PARALLEL
+     CALL wrf_error_fatal(TRIM(wrf_err_message))
+#else
+     PRINT *,TRIM(wrf_err_message)
+     STOP
+#endif
+   END IF
+! Checking for NaNs (should not be necessary but....)
+   IF (grid%psfc(im2,jm2) /= grid%psfc(im2,jm2) .OR. ABS(grid%psfc(im2,jm2)) > 1000000. ) THEN
+     PRINT errmsg
+     WRITE(wrF_err_message,*)'solve_em: wrong PSFC value=',               &
+       grid%psfc(im2,jm2),' at: ', im2 ,', ', jm2, ' !!!'
+#ifdef DM_PARALLEL
+     CALL wrf_error_fatal(TRIM(wrf_err_message))
+#else
+     PRINT *,TRIM(wrf_err_message)
+     STOP
+#endif
+   END IF
+
        IF (config_flags%lmdz_physics) THEN
 
@@ -865 +904 @@
         &             ,WNUM3DM=num_3d_m, WPARFIRSTSCAL=PARAM_FIRST_SCALAR,           &
         &        WNX=config_flags%e_we, WNY=config_flags%e_sn,                       &
-        &        WNZ=config_flags%e_vert, WJULDAY=grid%julday, WGMT=grid%gmt,        &
+        &        WNZ=config_flags%e_vert, WJULDAY=FLOAT(julday), WGMT=gmt,           &
         &        WTIME_STEP=REAL(config_flags%time_step),                            &
         &        WRF_FULLETA=grid%znw, WRF_HALFETA=grid%znu, WRF_DFULLETA=grid%dnw,  &