Changeset 1275 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

May 21, 2014, 6:09:22 PM (11 years ago)

Author:

msylvestre

Message:

LMDZ.GENERIC : For 1d model, the planetary radius can now be read in rcm1d.def. It is mandatory to set it when the rings_shadow option (callphys.def) is set to true. The message errors when a parameter is missing in rcm1d.def deals with rcm1d.def and not run.def any more

File:

• trunk/LMDZ.GENERIC/libf/phystd/rcm1d.F (modified) (9 diffs)

trunk/LMDZ.GENERIC/libf/phystd/rcm1d.F

@@ -356 +356 @@
       call getin("g",g)
       IF (g.eq.-99999.) THEN
-          PRINT *,"STOP. I NEED g IN RUN.DEF."
+          PRINT *,"STOP. I NEED g IN RCM1D.DEF."
           STOP
       ELSE
@@ -362 +362 @@
       ENDIF
 
-      !rad = -99999.
-      !PRINT *,'PLANETARY RADIUS in m ?'
-      !call getin("rad",rad)
-      !IF (rad.eq.-99999.) THEN
-      !    PRINT *,"STOP. I NEED rad IN RUN.DEF."
-      !    STOP
-      !ELSE
-      !    PRINT *,"--> rad = ",rad
-      !ENDIF
+      rad = -99999.
+      PRINT *,'PLANETARY RADIUS in m ?'
+      call getin("rad",rad)
+      ! Planetary  radius is needed to compute shadow of the rings
+      IF (rad.eq.-99999. .and. rings_shadow .eq. .true.) THEN
+          PRINT *,"STOP. I NEED rad IN RCM1D.DEF."
+          STOP
+      ELSE
+          PRINT *,"--> rad = ",rad
+      ENDIF
 
       daysec = -99999.
@@ -376 +377 @@
       call getin("daysec",daysec)
       IF (daysec.eq.-99999.) THEN
-          PRINT *,"STOP. I NEED daysec IN RUN.DEF."
+          PRINT *,"STOP. I NEED daysec IN RCM1D.DEF."
           STOP
       ELSE
@@ -389 +390 @@
       call getin("year_day",year_day)
       IF (year_day.eq.-99999.) THEN
-          PRINT *,"STOP. I NEED year_day IN RUN.DEF."
+          PRINT *,"STOP. I NEED year_day IN RCM1D.DEF."
           STOP
       ELSE 
@@ -399 +400 @@
       call getin("periastr",periastr)
       IF (periastr.eq.-99999.) THEN
-          PRINT *,"STOP. I NEED periastr IN RUN.DEF."
+          PRINT *,"STOP. I NEED periastr IN RCM1D.DEF."
           STOP
       ELSE
@@ -409 +410 @@
       call getin("apoastr",apoastr)
       IF (apoastr.eq.-99999.) THEN
-          PRINT *,"STOP. I NEED apoastr IN RUN.DEF."
+          PRINT *,"STOP. I NEED apoastr IN RCM1D.DEF."
           STOP
       ELSE
@@ -419 +420 @@
       call getin("peri_day",peri_day)
       IF (peri_day.eq.-99999.) THEN
-          PRINT *,"STOP. I NEED peri_day IN RUN.DEF."
+          PRINT *,"STOP. I NEED peri_day IN RCM1D.DEF."
           STOP
       ELSE IF (peri_day.gt.year_day) THEN
@@ -432 +433 @@
       call getin("obliquit",obliquit)
       IF (obliquit.eq.-99999.) THEN
-          PRINT *,"STOP. I NEED obliquit IN RUN.DEF."
+          PRINT *,"STOP. I NEED obliquit IN RCM1D.DEF."
           STOP
       ELSE
@@ -442 +443 @@
       call getin("psurf",psurf)
       IF (psurf.eq.-99999.) THEN
-          PRINT *,"STOP. I NEED psurf IN RUN.DEF."
+          PRINT *,"STOP. I NEED psurf IN RCM1D.DEF."
           STOP
       ELSE