Changeset 1156

Timestamp:

Jan 8, 2014, 11:49:31 AM (11 years ago)

Author:

emillour

Message:

Mars GCM:

• Update of the read_dust_scenario routine: when input dust scenarios file contain variable "dustop", it is assumed to be visible extinction opacity, and if it is "cdod" (recent change, for the MCDv5.1 dust scenarios), then it is IR absorption opacity (and is multiplied by 2.6 to be converted to visible extinction opacity).

EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/read_dust_scenario.F90 (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -1986 +1986 @@
   module "planetwide_mod", that enable obtaining the min and max of a field
   over the whole planet.
+
+== 08/01/2014 == EM
+- Update of the read_dust_scenario routine: when input dust scenarios file
+  contain variable "dustop", it is assumed to be visible extinction opacity,
+  and if it is "cdod" (recent change, for the MCDv5.1 dust scenarios), then
+  it is IR absorption opacity (and is multiplied by 2.6 to be converted to
+  visible extinction opacity).

trunk/LMDZ.MARS/libf/phymars/read_dust_scenario.F90

@@ -117 +117 @@
    allocate(lat(latlen), lon(lonlen), time(timelen))
 
+   ! "dustop" if loading visible extinction opacity
+   ! "cdod" if loading IR absorption opacity
    ierr=nf90_inq_varid(nid,"dustop",nvarid)
-   ierr=nf90_get_var(nid,nvarid,tautes)
-   IF (ierr .NE. nf90_noerr) THEN
-      PRINT*, "Error: read_dust_scenario <dustop> not found"
-      write(*,*)trim(nf90_strerror(ierr))
-      stop
-   ENDIF
+   if (ierr.eq.nf90_noerr) then
+     ierr=nf90_get_var(nid,nvarid,tautes)
+     IF (ierr .NE. nf90_noerr) THEN
+        PRINT*, "Error: read_dust_scenario <dustop> not found"
+        write(*,*)trim(nf90_strerror(ierr))
+        stop
+     ENDIF
+   else
+     ! did not find "dustop" , look for "cdod"
+     ierr=nf90_inq_varid(nid,"cdod",nvarid)
+     ierr=nf90_get_var(nid,nvarid,tautes)
+     IF (ierr .NE. nf90_noerr) THEN
+        PRINT*, "Error: read_dust_scenario <cdod> not found"
+        write(*,*)trim(nf90_strerror(ierr))
+        stop
+     ENDIF
+     ! and multiply by 2*1.3=2.6 to convert from IR absorption
+     ! to visible extinction opacity
+     tautes(:,:,:)=2.6*tautes(:,:,:)
+   endif
 
    ierr=nf90_inq_varid(nid,"Time",nvarid)