release of wrf 3.3

Makefile

@@ -4,6 +4,10 @@ LN      =       ln -s
 MAKE    =       make -i -r
 MV	=	/bin/mv
 RM      =       /bin/rm -f
+CHEM_FILES =	../chem/module_aerosols_sorgam.o \
+		../chem/module_gocart_aerosols.o \
+		../chem/module_mosaic_driver.o \
+		../chem/module_input_tracer.o
 
 deflt :
 		@ echo Please compile the code using ./compile
@@ -74,15 +78,21 @@ all_wrfvar :
 	$(MAKE) MODULE_DIRS="$(DA_WRFVAR_MODULES)" toolsdir
 	if [ $(CRTM) ] ; then \
 	  (cd var/external/crtm; \
-	  . configure/$(SFC).setup; make ) ; \
+	   export ABI_CRTM="${ABI_CRTM}"; . configure/$(SFC_CRTM).setup; $(MAKE) $(J) ) ; \
 	fi
 	if [ $(BUFR) ] ; then \
 	  (cd var/external/bufr;  \
-	  $(MAKE) FC="$(SFC)" CC="$(SCC)" CPP="$(CPP)" CFLAGS="$(CFLAGS)" FFLAGS="$(FCDEBUG) $(FORMAT_FIXED)" RANLIB="$(RANLIB)" AR="$(AR)" ARFLAGS="$(ARFLAGS)" ) ; \
+	  $(MAKE) $(J) FC="$(SFC)" CC="$(SCC)" CPP="$(CPP)" CPPFLAGS="$(CPPFLAGS)" CFLAGS="$(CFLAGS)" FFLAGS="$(FCDEBUG) $(FORMAT_FIXED)" RANLIB="$(RANLIB)" AR="$(AR)" ARFLAGS="$(ARFLAGS)" ) ; \
+	fi
+### Use 'make' to avoid '-i -r' above:
+	if [ $(WAVELET) ] ; then \
+	  ( cd var/external/wavelet; \
+		make CC="$(SC99) -DNOUNDERSCORE" RM="$(RM)" libWavelet.a; \
+		make FC="$(FC)" RM="$(RM)" lib_wavelet.a ) ; \
 	fi
 #	( cd var/build; touch depend.txt; make links; make depend; $(MAKE) $(J) all_wrfvar )
 	( cd var/build; make depend; $(MAKE) $(J) all_wrfvar )
-	( cd var/obsproc; $(MAKE) BUFR_CPP="$(BUFR_CPP)" )
+	( cd var/obsproc; $(MAKE) $(J) BUFR_CPP="$(BUFR_CPP)" )
 	@echo "build started:   $(START_OF_COMPILE)"
 	@echo "build completed:" `date`
 
@@ -109,6 +119,7 @@ em_fire : wrf
 	( cd test/em_fire ; /bin/rm -f README.namelist ; ln -s ../../run/README.namelist . )
 	( cd test/em_fire ; /bin/rm -f gribmap.txt ; ln -s ../../run/gribmap.txt . )
 	( cd test/em_fire ; /bin/rm -f grib2map.tbl ; ln -s ../../run/grib2map.tbl . )
+	( cd test/em_fire ; /bin/sh create_links.sh )
 	( cd run ; /bin/rm -f ideal.exe ; ln -s ../main/ideal.exe . )
 	( cd run ; if test -f namelist.input ; then \
 		/bin/cp -f namelist.input namelist.input.backup ; fi ; \
@@ -213,6 +224,20 @@ em_seabreeze2d_x : wrf
 	@echo "build started:   $(START_OF_COMPILE)"
 	@echo "build completed:" `date`
 
+em_tropical_cyclone : wrf
+	@ echo '--------------------------------------'
+	( cd main ; $(MAKE) MODULE_DIRS="$(ALL_MODULES)" SOLVER=em IDEAL_CASE=tropical_cyclone em_ideal )
+	( cd test/em_tropical_cyclone ; /bin/rm -f wrf.exe ; ln -s ../../main/wrf.exe . )
+	( cd test/em_tropical_cyclone ; /bin/rm -f ideal.exe ; ln -s ../../main/ideal.exe . )
+	( cd test/em_tropical_cyclone ; /bin/rm -f README.namelist ; ln -s ../../run/README.namelist . )
+	( cd run ; /bin/rm -f ideal.exe ; ln -s ../main/ideal.exe . )
+	( cd run ; if test -f namelist.input ; then \
+		/bin/cp -f namelist.input namelist.input.backup ; fi ; \
+		/bin/rm -f namelist.input ; ln -s ../test/em_tropical_cyclone/namelist.input . )
+	( cd run ; /bin/rm -f input_sounding ; ln -s ../test/em_tropical_cyclone/input_sounding . )
+	@echo "build started:   $(START_OF_COMPILE)"
+	@echo "build completed:" `date`
+
 em_scm_xy : wrf
 	@ echo '--------------------------------------'
 	( cd main ; $(MAKE) MODULE_DIRS="$(ALL_MODULES)" SOLVER=em IDEAL_CASE=scm_xy em_ideal )
@@ -518,7 +543,12 @@ physics :
 
 em_core :
 	@ echo '--------------------------------------'
-	( cd dyn_em ; $(MAKE) $(J) )
+	if [ $(WRF_CHEM) -eq 0 ] ; then \
+		CF= ; \
+	else \
+		CF=$(CHEM_FILES) ; \
+	fi
+	( cd dyn_em ; $(MAKE) $(J) CF="$(CF)" )
 
 # rule used by configure to test if this will compile with MPI 2 calls MPI_Comm_f2c and _c2f
 mpi2_test :

README

@@ -1,4 +1,4 @@
-WRF Model Version 3.2.1 (August 18, 2010)
+WRF Model Version 3.3 (April 6, 2011)
 http://wrf-model.org/users/users.php
 
 ------------------------
@@ -27,6 +27,18 @@ infringement actions.
 This is the main directory for the WRF Version 3 source code release.
 ======================================
 
+V3.3 Release Notes:
+-------------------
+
+Version 3.3 is released on April 6, 2011 (rev 4896).
+
+- For more information on WRF V3.3 release, visit WRF User's home pages
+  http://www.mmm.ucar.edu/wrf/users/, and
+  http://www.dtcenter.org/wrf-nmm/users/, and read the online User's Guide.
+
+
+======================================
+
 V3.2.1 Release Notes (8/18/10):
 ------------------
 
@@ -127,6 +139,8 @@ WRF update history:
 - V3.1: April 9, 2009
 - V3.1.1: July 31, 2009
 - V3.2: March 31, 2010
+- V3.2.1: August 18, 2010
+- V3.3: April 6, 2011
 
 ======================================
 
@@ -156,6 +170,7 @@ How to compile and run?
       compile em_squall2d_x
       compile em_squall2d_y
       compile em_scm_xy
+      compile em_tropical_cyclone
 
 - If sucessful, this will create either real.exe or ideal.exe and wrf.exe
   in directory main/, and the appropriate executables will be linked into
@@ -209,6 +224,7 @@ What is in WRF V3?
   * hydrostatic option via namelist option
   * positive-definite and monotonic advection for moisture, scalar, tke and chemical tracers
   * global modeling capability on latitude-longitude grid
+  * digital filter initialization
 
 - Two-way nesting:
   * multiple domains and multiple nest levels
@@ -224,17 +240,22 @@ What is in WRF V3?
 
 - Physics options:
   * microphysics (Kessler/ WRF Single Moment 3, 5 and 6 classes / Lin et al./ Mibrandt 2-moment /
-    Eta Ferrier / Thompson / Goddard / 2-moment Morrison / WRF Double Moment 5 and 6 classes )
+    Eta Ferrier / Thompson / Goddard / 2-moment Morrison / WRF Double Moment 5 and 6 classes /
+    SBU-Lin 5-classes )
   * cumulus parameterization (Kain-Fritsch with shallow convection / 
-    Betts-Miller-Janjic / Grell-Devenyi ensemble / New Grell (with shallow convection option) )
-  * planetary boundary layer (Yosei University / Mellor-Yamada-Janjic / ACM2 / QNSE / MYNN / BouLac )
-  * surface layer (similarity theory MM5 / Eta / QNSE / MYNN )
+    Betts-Miller-Janjic / Grell-Devenyi ensemble / New Grell (with shallow convection option) /
+    Tiedtke (with shallow conv and momentum transport) / NSAS (with shallow conv and momentum transport) /
+    SAS (with shallow conv for ARW) / Zhang-McFarlane (with momentum transport) )
+  * UW shallow convection
+  * planetary boundary layer (Yosei University / Mellor-Yamada-Janjic / ACM2 / QNSE / MYNN / BouLac /
+    UW / TEMF )
   * slab soil model (5-layer thermal diffusion / Noah land-surface model (4 levels) /
     RUC LSM (6 levels) / Pleim-Xu (2 levels)
   * Urban canopy model, BEP multi-layer and BEM  (works with Noah LSM, BEP and 
     BEM also requires MYJ or BouLac)
-  * longwave radiation (RRTM / CAM / RRTMG)
-  * shortwave radiation (simple MM5 scheme / Goddard / CAM / RRTMG ), terrain slope and shading effect
+  * longwave radiation (RRTM / CAM / RRTMG / new Goddard )
+  * shortwave radiation (simple MM5 scheme / old Goddard / CAM / RRTMG / new Goddard), 
+    terrain slope and shading effect
   * sub-grid turbulence (constant K diffusion/ 2-D Smagorinsky/ predicted TKE /
     2-D, 6th order diffusion / Nonlinear Backscatter Anisotropic (NBA) sub-grid turbulence stress for LES )
   * Rayleigh damping for w at the upper boundary layer

README.VAR → README.DA

@@ -22,6 +22,30 @@ WRFDA, including infringement actions.
 This is the main directory for the WRFDA Version 3 source code release.
 ======================================
 
+V3.3 Release Notes:
+-----------------------
+
+Version 3.3 is released on April 5, 2011.
+
+- RTTOV interfaces has been upgraded to version 10.0
+  Please download RTTOV v10.0 and install it in appropriate location 
+  if you want to use RTTOV to assimilate radiance observation with bufr format
+
+- crtm_atmosphere is removed from the namelist. 
+  Now the CRTM reference climatology is determined inside WRFDA based on 
+  latitude and month.
+
+- New WRF 4D-Var
+  The WRFPLUS (WRF adjoint and tangent linear model) has been upgraded to V3.3 
+  and it is consistent with the released WRFV version 3.3. However, only dry dynamic 
+  processes are included in WRFPLUS version 3.3 and only serial run is supported
+  at this release.
+  WRFDA is also upgraded to consistent with WRFPLUS version 3.3 and the new 4D-Var 
+  system is a single executable application. Please refer to Users' guide for the
+  details.
+
+- Lots of bug fixes.
+
 V3.2.1 Release Notes:
 -----------------------
 
@@ -143,7 +167,7 @@ How to compile and run?
         setenv NETCDF /usr/local/netcdf-pgi
         setenv BUFR 1 (optional, set to compile WRFDA with NCEP BUFR capability)
         setenv CRTM 1 (optional, set to compile WRFDA with CRTM capability)
-        setenv RTTOV /usr/local/rttov87 (optional, set to compile WRFDA with RTTOV capability)
+        setenv RTTOV /usr/local/rttov10 (optional, set to compile WRFDA with RTTOV capability)
 
 - Type 'compile all_wrfvar' when your are ready:
 

README.windturbine

@@ -0,0 +1,85 @@
+README for wind turbine drag parameterization schemes, added 20101220.  JM
+
+*Compiling the code.* no different from normal
+
+*Running the code.*  There are new variables added to the phys block of
+the namelist.input file.  The string valued variable windturbines_spec
+controls the overall operation of the scheme.  If it's set to "none"
+(quotes included) the wind turbine scheme is off.  If it's set to
+"ideal",  the turbine drag scheme is in idealized mode and other namelist
+variables in the block control the geometry and characteristics of the
+wind farm.  If it's set to the name of a file in the run directory --
+for example, "windspec.in" -- then the scheme is on and the location and
+characteristics of each turbine is set individually from the contents of
+that file.  Which scheme is active will be controlled by another variable
+TBD (not yet implemented, since there's only one scheme right now).
+
+*Idealized configuration.* If windturbines_spec is set to "ideal", the
+geometry and characteristics of the farm are controlled by the following
+additional variables in the phys block of the namelist.input file:
+
+ name        type     sample description
+                      value
+ turbgridid  integer   1     id of the WRF domain containing the farm
+ hubheight   real     100.   height in meters of all turbine hubs in farm
+ diameter    real     100.   diameter in meters of all turbine rotors in farm
+ stdthrcoef  real     0.158  standing thrust coefficient of all turbines in farm
+ cutinspeed  real     3.5    cut in speed (m/s) of all turbines in farm
+ cutoutspeed real     30.    cut out speed (m/s) of all turbines in farm
+ power       real     5.0    nominal power (MW) of all turbines in farm
+ turbpercell real     1.0    number of turbines per grid cell
+ ewfx        integer  10     x-extent of rectangular wind farm in grid cells
+ ewfy        integer  10     y-extent of rectangular wind farm in grid cells
+ pwfx        integer  45     x-coordinate of grid cell in sw corner of farm
+ pwfy        integer  45     y-coordinate of grid cell in sw corner of farm
+
+
+*Real-world configuration.*  The location and characteristics of each
+turbine is specified individually in the file named by the variable
+windturbines_spec.  Each line of the file specifies one turbine.
+The entries of a line are separated by spaces.  The entries are listed
+in order on the line and specify the following, by position on the line:
+
+     1. The WRF grid id into which the turbine is placed [integer]
+     2. True latitude of the wind turbine in degrees   [real]
+     3. True longitude of the wind turbine in degrees  [real]
+     4. Height in meters of the turbine hub [real]
+     5. Diameter in meters of the rotor     [real]
+     6. Standing thrust coefficient         [real]
+     7. Nominal power of turbine (MW)       [real]
+     8. Cutin speed (m/s) [real]
+     9. Cutout speed (m/s) [real]
+
+The location of each turbine is specified using the lat and lon elements
+for the turbine's entry in the windspec file.  That is, the i,j index
+in the grid is computed from the true latitude (entry 2) and the true
+longitude (entry 3).  The coordinate in the j dimension is computed by
+checking true latitude against the LAT_V field in the wrfinput initial
+conditions file.  The i dimension is computed by checking the true
+longitude against LONG_U field in the wrfinput file.
+
+*Pseudo-real configuration.*  This mode is only for testing the real-world
+specification mechanism.  The contents of the windspec.in file in the
+em_seabreeze2d_x directory contains a specification that is identical
+to the ideal specification shown above.  However, instead of true
+latitude and true longitude, elements 2 and 3 of each entry specify
+the i,j grid coordinates of the turbine.  The code knows to treat
+these as grid coordinates instead of latitude and longitude because
+for idealized WRF cases, the LAT and LONG fields in wrfinput files are
+set to all zero.  When the model reads a wrfinput file for an idealized
+case, it checks to see if the lat and lon fields of the case have data
+in them or if they are just zeros.  If they do not have data in them
+(just have zeros) fields 2 and 3 of the windspec file are treated as i,j
+coordinates are used.  A sample windspec file for a pseudo-real run is
+in test/em_seabreeze2d_x/windspec.in).
+
+*Specific note for the Fitch scheme.*  The Fitch scheme is based on
+Blahak et al, 2010, of Wetter Jetzt GbR.  The specific implementation
+comes from Anna Fitch, Alok Gupta, and Idar Barstad at Uni Bergen,
+Norway.  It was added to this release of WRF by Jimy Dudhia, Joe Olson
+(NOAA), Julie Lundquist (U. Colorado/NREL), and John Michalakes (NREL).
+It works with the MYNN PBL and modifies the QKE field representing 2xTKE
+to include the TKE produced by wind turbines. QKE is advected as a part
+of the scalar 4D tracer array in WRF.
+
+

README_test_cases

@@ -1,4 +1,4 @@
-WRFV2 (Weather Research and Forecast) model.   
+WRFV3 (Weather Research and Forecast) model.   
 
 Contents:
 
@@ -108,3 +108,13 @@ The available test cases are
 
 10) Single column model (test/em_scm_xy)
 
+11) 3D tropical cyclone (test/em_tropical_cyclone)
+
+   Idealized tropical cyclone on an f-plane with constant SST in a 
+   specified environment.  Uses capped Newtonian relaxation to mimic 
+   longwave radiation.  The default setup uses the Jordan (1958) mean 
+   hurricane sounding with 28 C SST and no background winds.  The 
+   analytic initial vortex is from Rotunno and Emanuel (1987).  This 
+   case is useful for testing the effects of new model code (e.g., 
+   new physics options) on tropical cyclones in an idealized framework. 
+

Registry/Registry.CONVERT

@@ -470,6 +470,7 @@ package   wsm5scheme    mp_physics==4                -             moist:qv,qc,q
 package   etampnew      mp_physics==5                -             moist:qv,qc
 package   wsm6scheme    mp_physics==6                -             moist:qv,qc,qr,qi,qs,qg
 package   thompson       mp_physics==8               -             moist:qv,qc,qr,qi,qs,qg;scalar:qni
+package   sbu_ylinscheme    mp_physics==13           -             moist:qv,qc,qr,qi,qs
 
 package   rrtmscheme    ra_lw_physics==1             -             -
 package   gfdllwscheme  ra_lw_physics==99            -             moist:qv,qc,qr,qi