[Dart-dev] DART/branches Revision: 12252

[dart at ucar.edu]

thoar at ucar.edu
2017-12-29 15:31:10 -0700 (Fri, 29 Dec 2017)
257
Successfully applying the forward operators for P_SFC_LAND, P_SFC_SHIP, various T obs, U, V obs ... have not tested Z_RAOB or WS_SSMI_FF1.
Need to have a geographic subsetting program.
I believe the failing operators are those that are outside the domain.




Modified: DART/branches/coamps/models/coamps_nest/coamps_intrinsic_mod.f90
===================================================================
--- DART/branches/coamps/models/coamps_nest/coamps_intrinsic_mod.f90	2017-12-29 20:56:34 UTC (rev 12251)
+++ DART/branches/coamps/models/coamps_nest/coamps_intrinsic_mod.f90	2017-12-29 22:31:10 UTC (rev 12252)
@@ -704,15 +704,12 @@
 !
 !  compute surface pressure (mb)
 !
-      do i=1,mn   !>@todo ps(m,n) and mn loop ps(mn,1)
+      do i=1,mn
         aoz=sigmwa(1)/(sigmwa(1)-zsfc(i,1))
         ppsfc=p(i,1,kk)&
              -dsigw(kkp1)*g*(th(i,1,kk)-thbm(i,1,kk))&
             /(aoz*cp*thbm(i,1,kk)*thbm(i,1,kk))
-
-        print*,    exbw(i,1,kkp1),ppsfc,   cpor, i, m, n, kk, kkp1
-        ps(i,1) = ((exbw(i,1,kkp1)+ppsfc)**cpor)*1000.0_r8
-        print*,    ps(i,1)
+        ps(i,1) = (exbw(i,1,kkp1)+ppsfc)**cpor * 1000.0_r8
       enddo
 
       return

Modified: DART/branches/coamps/models/coamps_nest/work/input.nml
===================================================================
--- DART/branches/coamps/models/coamps_nest/work/input.nml	2017-12-29 20:56:34 UTC (rev 12251)
+++ DART/branches/coamps/models/coamps_nest/work/input.nml	2017-12-29 22:31:10 UTC (rev 12252)
@@ -158,7 +158,9 @@
                        'U_WSAT_WIND' , 'V_WSAT_WIND',
                        'U_ASCAT_WIND', 'V_ASCAT_WIND',
                        'U_TERRA_WIND', 'V_TERRA_WIND',
-                       'T_RAOB', 'T_AMDAR', 'T_MDCRS'
+                       'T_RAOB', 'T_AMDAR', 'T_MDCRS', 
+                       'WS_SSMI_FF1' , 'T_TOVS_T',
+                       'P_SFC_LAND'  , 'P_SFC_SHIP'
    keep_types        = .true.
    print_only        = .false.
    min_lat           = -90.0
@@ -167,7 +169,7 @@
    max_lon           = 360.0
    /
 
-# 'T_AMDAR', 'T_MDCRS', 'Z_RAOB' , 'P_SFC_LAND'  , 'P_SFC_SHIP', 'WS_SSMI_FF1' , 'T_TOVS_T'
+# 'Z_RAOB', 'WS_SSMI_FF1'
 
 &preprocess_nml
     input_obs_kind_mod_file = '../../../assimilation_code/modules/observations/DEFAULT_obs_kind_mod.F90'
@@ -314,9 +316,9 @@
    Nregions   = 1
    lonlim1    =   0.0,   0.0,   0.0, 235.0
    lonlim2    = 360.0, 360.0, 360.0, 295.0
-   latlim1    =  20.0, -80.0, -20.0,  25.0
-   latlim2    =  80.0, -20.0,  20.0,  55.0
-   reg_names  = 'Northern Hemisphere', 'Southern Hemisphere', 'Tropics', 'North America'
+   latlim1    = -90.0, -80.0, -20.0,  25.0
+   latlim2    =  90.0, -20.0,  20.0,  55.0
+   reg_names  = 'Earth', 'Southern Hemisphere', 'Tropics', 'North America'
    print_mismatched_locs = .false.
    create_rank_histogram = .true.
    outliers_in_histogram = .true.
@@ -358,7 +360,7 @@
    innov_file_name = 'innov.out'
    ngt_file_name   = 'ngt.out'
    obs_seq_in_name = 'obs_seq.out'
-   obs_window      = 240
+   obs_window      = 60
    verbose         = .true.
   /