--- trunk/dyn3d/serre.f	2014/03/05 14:57:53	82
+++ trunk/dyn3d/serre.f	2014/09/18 13:36:51	112
@@ -21,4 +21,36 @@
   real:: taux = 3. ! raideur du zoom en X
   real:: tauy = 3. ! raideur du zoom en Y
 
+contains
+
+  subroutine read_serre
+
+    use unit_nml_m, only: unit_nml
+
+    namelist /serre_nml/ clon, clat, grossismx, grossismy, dzoomx, dzoomy, &
+         taux, tauy
+
+    !-------------------------------------------------
+
+    print *, "Enter namelist 'serre_nml'."
+    read(unit=*, nml=serre_nml)
+    write(unit_nml, nml=serre_nml)
+
+    IF (grossismx < 1.) THEN
+       PRINT *, 'Error: grossismx < 1'
+       STOP 1
+    ELSE
+       alphax = 1. - 1. / grossismx
+    ENDIF
+    IF (grossismy < 1.) THEN
+       PRINT *, 'Error: grossismy < 1'
+       STOP 1
+    ELSE
+       alphay = 1. - 1. / grossismy
+    ENDIF
+    PRINT *, 'alphax = ', alphax
+    PRINT *, 'alphay = ', alphay
+
+  end subroutine read_serre
+
 end module serre